Aggregate impacts in Saturn's rings

Ring particle aggregates are formed in the outer parts of Saturn's main rings. We study how collisions between aggregates can lead to destruction or coalescence of these aggregates, with local N-body simulations taking into account the dissipative impacts and gravitational forces between particles. Impacts of aggregates with different mass ratios are studied, as well as aggregates that consist of particles with different physical properties. We find that the outcome of the collision is very sensitive to the shape of the aggregate, in the sense that more elongated aggregates are more prone to be destroyed. We were interested in testing the accretion criterion Barbara and Esposito [Barbara, J.M., Esposito, L.W., 2002. Icarus 160, 161-171] used in their F ring simulations, according to which accretion requires that the masses of the colliding bodies differ at least by a factor of 100. We confirm that such a critical mass ratio exists. In particular, simulations indicate that the exact critical mass ratio depends on the internal density and elasticity of particles, and the distance from the planet. The zone of transition, defined by the distance where individual particles or small aggregates first start to stick on the larger aggregate, and by the distance where two similar sized aggregates on the average eventually coalesce is only about 5000 km wide, if fixed particle properties are used. The rotational state of the aggregates that form via aggregate collision rapidly reaches synchronous rotation, similarly to the aggregates that form via gradual growth.     

Spectral fit of Cygnus X-1 in high energy— a self-consistent study

We fit the spectra of Cyg X-1 using two component advective flows with Keplerian accretion disks on the equatorial plane surrounded by sub-Keplerian disks when standing shocks are present. The soft photons generated by the bremsstrahlung and synchrotron processes in the sub-Keplerian flow, as well as the multi-colour black body emission from the Keplerian disk are Comptonized by the thermal and non-thermal electrons. By varying Keplerian and sub-Keplerian rates we are able to reproduce the observed soft and hard states as far as X-ray region is concerned and ‘low γ-ray intensity’ and ‘high γ-ray intensity’ states as far as the soft γ-ray region is concerned. We also find two pivotal points where the spectra intersect as is observed in Cyg X-1.      

Black hole spin dependence of general relativistic multi-transonic accretion close to the horizon

We introduce a novel formalism to investigate the role of the spin angular momentum of astrophysical black holes in influencing the behavior of low angular momentum general relativistic accretion. We propose a metric independent analysis of axisymmetric general relativistic flow, and consequently formulate the space and time dependent equations describing the general relativistic hydrodynamic accretion flow in the Kerr metric. The associated stationary critical solutions for such flow equations are provided and the stability of the stationary transonic configuration is examined using an elegant linear perturbation technique. We examine the properties of infalling material for both prograde and retrograde accretion as a function of the Kerr parameter at extremely close proximity to the event horizon. Our formalism can be used to identify a new spectral signature of black hole spin, and has the potential of performing the black hole shadow imaging corresponding to the low angular momentum accretion flow.     

Particle stirring in turbulent gas disks: Including orbital oscillations

We describe the diffusion and random velocities of solid particles due to stochastic forcing by turbulent gas. We include the orbital dynamics of Keplerian disks, both in-plane epicycles and vertical oscillations. We obtain a new result for the diffusion of solids. The Schmidt number (ratio of gas to particle diffusivity) is Sc≈1+(Ωt_stop)², in terms of the particle stopping time t_stop and the orbital frequency Ω. The standard result, Sc=1+t_stop/t_eddy, in terms of the eddy turnover time, t_eddy, is shown to be incorrect. The main difference is that Sc rises quadratically, not linearly, with stopping time. Consequently, particles larger than 10 cm in protoplanetary disks will suffer less radial diffusion and will settle closer to the midplane. Such a layer of boulders would be more prone to gravitational collapse. Our predictions of RMS speeds, vertical scale height and diffusion coefficients will help interpret numerical simulations. We confirm previous results for the vertical stirring of particles (scale heights and random velocities), and add a correction for arbitrary ratios of eddy to orbital times. The particle layer becomes thinner for t_eddy>1/Ω with the strength of turbulent diffusion held fixed. We use two analytic techniques—the Hinze–Tchen formalism and the Fokker–Planck equation with velocity diffusion—with identical results when the regimes of validity overlap. We include simple physical arguments for the scaling of our results.     

Fission track and fault kinematics analyses for new insight into the Late Cenozoic tectonic regime changes in West-Central Sulawesi (Indonesia)

A 3-D density model for the Cretan and Libyan Seas and Crete was developed by gravity modelling constrained by five 2-D seismic lines. Velocity values of these cross-sections were used to obtain the initial densities using the Nafe–Drake and Birch empirical functions for the sediments, the crust and the upper mantle. The crust outside the Cretan Arc is 18 to 24 km thick, including 10 to 14 km thick sediments. The crust below central Crete at its thickest section, has values between 32 and 34 km, consisting of continental crust of the Aegean microplate, which is thickened by the subducted oceanic plate below the Cretan Arc. The oceanic lithosphere is decoupled from the continental along a NW–SE striking front between eastern Crete and the Island of Kythera south of Peloponnese. It plunges steeply below the southern Aegean Sea and is probably associated with the present volcanic activity of the southern Aegean Sea in agreement with published seismological observations of intermediate seismicity. Low density and velocity upper mantle below the Cretan Sea with ρ  3.25 × 10³ kg/m³ and V_p velocity of compressional waves around 7.7 km/s, which are also in agreement with observed high heat flow density values, point out at the mobilization of the upper mantle material here. Outside the Hellenic Arc the upper mantle density and velocity are ρ ≥ 3.32 × 10³ kg/m³ and V_p = 8.0 km/s, respectively. The crust below the Cretan Sea is thin continental of 15 to 20 km thickness, including 3 to 4 km of sediments. Thick accumulations of sediments, located to the SSW and SSE of Crete, are separated by a block of continental crust extended for more than 100 km south of Central Crete. These deep sedimentary basins are located on the oceanic crust backstopped by the continental crust of the Aegean microplate. The stretched continental margin of Africa, north of Cyrenaica, and the abruptly terminated continental Aegean microplate south of Crete are separated by oceanic lithosphere of only 60 to 80 km width at their closest proximity. To the east and west, the areas are floored by oceanic lithosphere, which rapidly widens towards the Herodotus Abyssal plain and the deep Ionian Basin of the central Mediterranean Sea. Crustal shortening between the continental margins of the Aegean microplate and Cyrenaica of North Africa influence the deformation of the sediments of the Mediterranean Ridge that has been divided in an internal and external zone. The continental margin of Cyrenaica extends for more than 80 km to the north of the African coast in form of a huge ramp, while that of the Aegean microplate is abruptly truncated by very steep fractures towards the Mediterranean Ridge. Changes in the deformation style of the sediments express differences of the tectonic processes that control them. That is, subduction to the northeast and crustal subsidence to the south of Crete. Strike-slip movement between Crete and Libya is required by seismological observations.     

Models of Jupiter's growth incorporating thermal and hydrodynamic constraints

We model the growth of Jupiter via core nucleated accretion, applying constraints from hydrodynamical processes that result from the disk-planet interaction. We compute the planet's internal structure using a well tested planetary formation code that is based upon a Henyey-type stellar evolution code. The planet's interactions with the protoplanetary disk are calculated using 3-D hydrodynamic simulations. Previous models of Jupiter's growth have taken the radius of the planet to be approximately one Hill sphere radius, RH. However, 3-D hydrodynamic simulations show that only gas within ∼0.25RH remains bound to the planet, with the more distant gas eventually participating in the shear flow of the protoplanetary disk. Therefore in our new simulations, the planet's outer boundary is placed at the location where gas has the thermal energy to reach the portion of the flow not bound to the planet. We find that the smaller radius increases the time required for planetary growth by ∼5%. Thermal pressure limits the rate at which a planet less than a few dozen times as massive as Earth can accumulate gas from the protoplanetary disk, whereas hydrodynamics regulates the growth rate for more massive planets. Within a moderately viscous disk, the accretion rate peaks when the planet's mass is about equal to the mass of Saturn. In a less viscous disk hydrodynamical limits to accretion are smaller, and the accretion rate peaks at lower mass. Observations suggest that the typical lifetime of massive disks around young stellar objects is ∼3 Myr. To account for the dissipation of such disks, we perform some of our simulations of Jupiter's growth within a disk whose surface gas density decreases on this timescale. In all of the cases that we simulate, the planet's effective radiating temperature rises to well above 1000 K soon after hydrodynamic limits begin to control the rate of gas accretion and the planet's distended envelope begins to contract. According to our simulations, proto-Jupiter's distended and thermally-supported envelope was too small to capture the planet's current retinue of irregular satellites as advocated by Pollack et al. [Pollack, J.B., Burns, J.A., Tauber, M.E., 1979. Icarus 37, 587-611].     

Structure of self-gravity wakes in Saturn's A ring as measured by Cassini CIRS

The CIRS infrared spectrometer onboard the Cassini spacecraft has scanned Saturn's A ring azimuthally from several viewing angles since its orbit insertion in 2004. A quadrupolar asymmetry has been detected in this ring at spacecraft elevations ranging between 16° to 37°. Its fractional amplitude decreases from 22% to 8% from 20° to 37° elevations. The patterns observed in two almost complete azimuthal scans at elevations 20° and 36° strongly favor the self-gravity wakes as the origin of the asymmetry. The elliptical, infinite cylinder model of Hedman et al. [Hedman, M.M., Nicholson, P.D., Salo, H., Wallis, B.D., Buratti, B.J., Baines, K.H., Brown, R.H., Clark, R.N., 2007. Astron. J. 133, 2624-2629] can reproduce the CIRS observations well. Such wakes are found to have an average height-to-spacing ratio H/λ=0.1607±0.0002, a width-over-spacing W/λ=0.3833±0.0008. Gaps between wakes, which are filled with particles, have an optical depth τG=0.1231±0.0005. The wakes mean pitch angle ΦW is 70.70°±0.07°, relative to the radial direction. The comparison of ground-based visible data with CIRS observations constrains the A ring to be a monolayer. For a surface mass density of 40 g cm⁻² [Tiscarino, M.S., Burns, J.A., Nicholson, P.D., Hedman, M.M., Porco, C.C., 2007. Icarus 189, 14-34], the expected spacing of wakes is λ≈60 m. Their height and width would then be H≈10 m and W≈24 m, values that match the maximum size of particles in this ring as determined from ground-based stellar occultations [French, R.G., Nicholson, P.D., 2000. Icarus 145, 502-523].     

Saturn in hot water: Viscous evolution of the Enceladus torus

The detection of outgassing water vapor from Enceladus is one of the great breakthroughs of the Cassini mission. The fate of this water once ionized has been widely studied; here we investigate the effects of purely neutral-neutral interactions within the Enceladus torus. We find that, thanks in part to the polar nature of the water molecule, a cold (∼180 K) neutral torus would undergo rapid viscous heating and spread to the extent of the observed hydroxyl cloud, before plasma effects become important. We investigate the physics behind the spreading of the torus, paying particular attention to the competition between heating and rotational line cooling. A steady-state torus model is constructed, and it is demonstrated that the torus will be observable in the millimeter band with the upcoming Herschel satellite. The relative strength of rotational lines could be used to distinguish between physical models for the neutral cloud.     

The molecular hydrogen explorer H2EX

The Molecular Hydrogen Explorer, H2EX, was proposed in response to the ESA 2015 - 2025 Cosmic Vision Call as a medium class space mission with NASA and CSA participations. The mission, conceived to understand the formation of galaxies, stars and planets from molecular hydrogen, is designed to observe the first rotational lines of the H₂ molecule (28.2, 17.0, 12.3 and 9.7 μm) over a wide field, and at high spectral resolution. H2EX can provide an inventory of warm (≥ 100 K) molecular gas in a broad variety of objects, including nearby young star clusters, galactic molecular clouds, active galactic nuclei, local and distant galaxies. The rich array of molecular, atomic and ionic lines, as well as solid state features available in the 8 to 29 μm spectral range brings additional science dimensions to H2EX. We present the optical and mechanical design of the H2EX payload based on an innovative Imaging Fourier Transform Spectrometer fed by a 1.2 m telescope. The 20’×20’ field of view is imaged on two 1024×1024 Si:As detectors. The maximum resolution of 0.032 cm^− 1 (full width at half maximum) means a velocity resolution of 10 km s^− 1 for the 0 – 0 S(3) line at 9.7 μm. This instrument offers the large field of view necessary to survey extended emission in the Galaxy and local Universe galaxies as well as to perform unbiased extragalactic and circumstellar disks surveys. The high spectral resolution makes H2EX uniquely suited to study the dynamics of H₂ in all these environments. The mission plan is made of seven wide-field spectro-imaging legacy programs, from the cosmic web to galactic young star clusters, within a nominal two years mission. The payload has been designed to re-use the Planck platform and passive cooling design.

Clumps in the outer disk by disk instability: Why they are initially gas giants and the legacy of disruption

We explore the initial conditions for fragments in the extended regions of gravitationally unstable disks. We combine analytic estimates for the fragmentation of spiral arms with 3D SPH simulations to show that initial fragment masses are in the gas giant regime. These initial fragments will have substantial angular momentum, and should form disks with radii of a few AU. We show that clumps will survive for multiple orbits before they undergo a second, rapid collapse due to H₂ dissociation and that it is possible to destroy bound clumps by transporting them into the inner disk. The consequences of disrupted clumps for planet formation, dust processing, and disk evolution are discussed. We argue that it is possible to produce Earth-mass cores in the outer disk during the earliest phases of disk evolution.