Skip to main content
Open Access Publications from the University of California

UCSC is one of the world's leading centers for both observational and theoretical research in astronomy and astrophysics. The department was recently ranked first in the country in research impact, based on citation studies. Faculty and students in the department and our affiliated research centers are building and using first-rank telescopes and instrumentation—on Earth and in space—extending humanity’s vision to planets orbiting nearby stars and the first stirrings of the Universe.

The department includes 24 faculty members, whose research interests range from our solar system and the Milky Way to the most distant galaxies in the Universe and the most fundamental questions of cosmology.

UCSC is a leader in astrophysics education, and we attract some the best graduate students in the country, enrolling approximately 40 students working towards the Ph.D. degree.

Currently this page is for hosting only ISIMA (International Summer Institute for Modeling in Astrophysics) conference proceedings.

Cover page of Singular tidal modes and the regularization of the tidal singularity

Singular tidal modes and the regularization of the tidal singularity


Following the work on ray orbits in spatially hyperbolic systems by Mass and Lam (1995) and Rieutord and Valdettaro (1997) we seek to examine the behaviour of the shear layer emitted at the critical latitude in 3D in a spherical shell fllled with rotating fluid. We compare the (previously known) 3D and the 2D solutions for a sphere in an infinite domain to find the major difference being a logarithmic singularity on the rotation axis formed by a cone of shear converging to an apex. We then consider the "split disc" arrangement first considered by Walton to examine this singularity in more detail. We also consider the behaviour of the Moore and Saffman shear layers under the influence of a large-scale forcing; our motivation is primarily the dissipation of tidal energiesin astrophysical binary systems.

Cover page of Taming jets in magnetised fluids

Taming jets in magnetised fluids


The effects of a uniform horizontal magnetic field on jets dynamics in 2D Boussinesq turbulence, i.e. Howard-Krishnamurti problem are studied with a numerical simulation. For a fixed fluid and magnetic diffusivity, it is shown that as the imposed field strength becomes larger jets start behaving in a more organized way, i.e. achieve stationary state and are finally quenched. The time evolution of total stress, Reynolds stress, Maxwell stress is examined and all the stresses are shown to vanish when jets are quenched. The quenching of jets is confirmed for different values of magnetic diffusivity, albeit the required field strength increases. It is also shown that the inclusion of overstable modes reinforces jets where Maxwell stress overcomes Reynolds stress. For a larger imposed field jets are shown to quench. A possible mechanism for the transition to the reinforcement of jets by Maxwell stress is discussed based on the transition in the most unstable mode in the underlying turbulence.

Cover page of The orbital decay of a retrograde planet in a protoplanetary disk

The orbital decay of a retrograde planet in a protoplanetary disk


Motivated by recent observations of retrograde planets, we investigated the orbital decay of a retrograde planet embedded in a protoplanetary disk. We treated both gravitational and hydrodynamic drag, and found the migration time scale ranges from 103 to 105 years for planet masses between 10-3 to 101 Jupiter masses. We also found that a highly inclined orbit can increase this time scale by a factor of 10 and that due to inclination damping, the final inclination is unlikely to be greater than 50 degrees.

Cover page of Radiative Rayleigh-Taylor instabilities

Radiative Rayleigh-Taylor instabilities


This project investigates the role of radiation in Rayleigh-Taylor instabilities by performing linear stability analyses of a plane parallel background equilibrium, with a semi-infinite medium 1 overlying a semi-infinite medium 2, in a gravitational field g and a radiation flux F normal to the discontinuity.

Cover page of Geostrophic turbulence with a magnetic field

Geostrophic turbulence with a magnetic field


The project is an extension of the work on f-plane magnetohydrodynamic (MHD) turbulence and its consequences on momentum transport. A somewhat detailed overview is given, with the physical mechanisms explained. The quasi-geostrophic equations, so well known in the Geophysical Fluid Dynamics (GFD) community, is derived with the Lorentz force present. The two-layer model is proposed as a simplified model for our studies. Progress with magnetically influenced barotropic and baroclinic instabilities are given, and some proposed future work concludes the document.

Cover page of MHD jet propagation in the case of DG Tau

MHD jet propagation in the case of DG Tau


The aim of the work is to perform numerical simulations of the propagation of stellar jets with consistent nozzle conditions obtained from launching simulations. This novel approach provides a global picture of the jet from its launching to its interaction with the ambient medium. The flow parameters observed at a distance of a few AU from the protostellar jet DG Tau were used to constrain the global inflow conditions whereas the actual profiles of different quantities are obtained from steady-state launching simulations. A new simulation was run on time and length scales typical of stellar jets. We also investigated the effects of cooling in these jets. We find evidence of density knots in our adiabatic simulations whereas simulations with cooling have much fewer and weaker knots.

Cover page of Stoked Dynamos

Stoked Dynamos


In this project we address the question of whether a flow that is not a dynamo can be made to exhibit dynamo-like properties by feeding it with a small amount of magnetic field. This may be pertinent to the solar dynamo and the processes that sustain it. We present a 3-D fully nonlinear magnetohydrodynamic simulation of the dynamo properties of a time-dependent ABC flow and discuss a method for leaking magnetic field into the computational domain. Our results suggest that sufficient magnetic feeding significantly boosts the magnetic energy of nondynamo flows and can maintain a mangetic field for long times.

Cover page of Production of Elephant Trunks in HII Regions by Radiation-Magnetohydrodynamic Instabilities

Production of Elephant Trunks in HII Regions by Radiation-Magnetohydrodynamic Instabilities


Recent SPH and grid code simulations showed, that ionizing radiation can amplify overdensities in turbulent molecular clouds and produce molecular pillars. The relevance of magnetic fields for the structure and stability of molecular clouds is still under discussion. We investigate whether an ionization front hitting a medium with small distortions of the magnetic field can produce the observed pillar-like structures in star forming regions (e.g. Eagle Nebula). Numerical MHD simulations with the Athena 2.0 grid code with ionizing radiation were performed. It turns out that the ionizing radiation drives a shock wave into the cold magnetized cloud and amplifies overdensities seeded by Alfven waves. Alfven waves can be seeds for molecular pillars. However, the magnetic field in structures created by Alfven waves makes these regions hostile to star formation.

Cover page of The role of radiation pressure in the dynamics of HII regions at z>1

The role of radiation pressure in the dynamics of HII regions at z>1


Observations of starburst galaxy at high redshift hint that the ionization parameter at z ~ 2 is higher than in the local universe. Following Krumholz & Matzner (2009), a physical explanation of a higher ionization parameter can be a radiation pressure-dominated HII region population. We wrote a population synthesis code to generate a family of HII regions and let them evolve following a solution that acconts for both radiation pressure- and gas pressure-dominated evolution. We suppose that the galaxy is spatially unresolved, and that the star formation rate and the ambient density are constant. We find that the ionization parameter increases for Lyman-Break galaxies.

Cover page of Sweet-Parker Reconnection with Anomalous Resistivity — A Toy Model

Sweet-Parker Reconnection with Anomalous Resistivity — A Toy Model


Magnetic reconnection is a common phenomenon in astrophysical contexts. The conventional Sweet-Parker model describes magnetic reconnection due resistivity. However, microscopic resistivity appears too small to reproduce the observed rate of reconnection. In this report, we describe the basic idea of anomalous resistivity in non-relativistic collisionless ion-electron plasma. We build a one-dimensional model along the direction of current in the current sheet. When the ion temperature is much less than the electron temperature, ion-acoustic instability develops when current density is sufficiently large so that the electron drift speed exceeds a few times the sound speed. The instability generates ion-acoustic waves, which are damped by non-linear wave-particle interaction. Anomalous resistivity arises due to the momentum exchange between waves and particles. The calculated anomalous resistivity strongly depends on the current density in the current sheet, and is typically much larger than the microscopic resistivity. However, matching the anomalous resistivity to the Sweet-Parker model, the resulting reconnection rate still falls off the observed rate by a large factor.