Fiducial Case Research Articles

Statistical gravitational waveform models: What to simulate next?

Models of gravitational waveforms play a critical role in detecting and characterizing the gravitational waves (GWs) from compact binary coalescences. Waveforms from numerical relativity (NR), while highly accurate, are too computationally expensive to produce to be directly used with Bayesian parameter estimation tools like Markov-chain-Monte-Carlo and nested sampling. We propose a Gaussian process regression (GPR) method to generate accurate reduced-order-model waveforms based only on existing accurate (e.g. NR) simulations. Using a training set of simulated waveforms, our GPR approach produces interpolated waveforms along with uncertainties across the parameter space. As a proof of concept, we use a training set of IMRPhenomD waveforms to build a GPR model in the 2-d parameter space of mass ratio $q$ and equal-and-aligned spin $\chi_1=\chi_2$. Using a regular, equally-spaced grid of 120 IMRPhenomD training waveforms in $q\in[1,3]$ and $\chi_1 \in [-0.5,0.5]$, the GPR mean approximates IMRPhenomD in this space to mismatches below $4.3\times 10^{-5}$. Our approach can alternatively use training waveforms directly from numerical relativity. Beyond interpolation of waveforms, we also present a greedy algorithm that utilizes the errors provided by our GPR model to optimize the placement of future simulations. In a fiducial test case we find that using the greedy algorithm to iteratively add simulations achieves GPR errors that are $\sim 1$ order of magnitude lower than the errors from using Latin-hypercube or square training grids.

Fates of the dense cores formed by fragmentation of filaments: do they fragment again or not?

Fragmentation of filaments into dense cores is thought to be an important step in forming stars. The bar-mode instability of spherically collapsing cores found in previous linear analysis invokes a possibility of re-fragmentation of the cores due to their ellipsoidal (prolate or oblate) deformation. To investigate this possibility, here we perform three-dimensional self-gravitational hydrodynamics simulations that follow all the way from filament fragmentation to subsequent core collapse. We assume the gas is polytropic with index \gamma, which determines the stability of the bar-mode. For the case that the fragmentation of isolated hydrostatic filaments is triggered by the most unstable fragmentation mode, we find the bar mode grows as collapse proceeds if \gamma < 1.1, in agreement with the linear analysis. However, it takes more than ten orders-of-magnitude increase in the central density for the distortion to become non-linear. In addition to this fiducial case, we also study non-fiducial ones such as the fragmentation is triggered by a fragmentation mode with a longer wavelength and it occurs during radial collapse of filaments and find the distortion rapidly grows. In most of astrophysical applications, the effective polytropic index of collapsing gas exceeds 1.1 before ten orders-of-magnitude increase in the central density. Thus, supposing the fiducial case of filament fragmentation, re-fragmentation of dense cores would not be likely and their final mass would be determined when the filaments fragment.

Star Cluster Formation from Turbulent Clumps. I. The Fast Formation Limit

Abstract We investigate the formation and early evolution of star clusters, assuming that they form from a turbulent starless clump of a given mass bounded inside a parent self-gravitating molecular cloud characterized by a particular mass surface density. As a first step, we assume instantaneous star cluster formation and gas expulsion. We draw our initial conditions from observed properties of starless clumps. We follow the early evolution of the clusters up to 20 Myr, investigating the effects of different star formation efficiencies, primordial binary fractions and eccentricities, and primordial mass segregation levels. We investigate clumps with initial masses of embedded in ambient cloud environments with mass surface densities and . We show that these models of fast star cluster formation result, in the fiducial case, in clusters that expand rapidly, even considering only the bound members. Clusters formed from higher environments tend to expand more quickly and thus are soon larger than clusters born from lower conditions. To form a young cluster of a given age, stellar mass, and mass surface density, these models need to assume a parent molecular clump that is many times denser, which is unrealistic compared to observed systems. We also show that, in these models, the initial binary properties are only slightly modified by interactions, meaning that the binary properties, e.g., at 20 Myr, are very similar to those at birth. With this study, we set up the foundation for future work, where we will investigate more realistic models of star formation compared to this instantaneous, baseline case.

Molecular clouds photoevaporation and FIR line emission

With the aim of improving predictions on far infrared (FIR) line emission from Giant Molecular Clouds (GMC), we study the effects of photoevaporation (PE) produced by external far-ultraviolet (FUV) and ionizing (extreme-ultraviolet, EUV) radiation on GMC structure. We consider three different GMCs with mass in the range $M_{\rm GMC} = 10^{3-6}\,\rm{M_{\odot}}$. Our model includes: (i) an observationally-based inhomogeneous GMC density field, and (ii) its time evolution during the PE process. In the fiducial case ($M_{\rm GMC}\approx10^5 M_{\odot}$), the photoevaporation time ($t_{pe}$) increases from 1 Myr to 30 Myr for gas metallicity $Z=0.05-1\,\rm Z_{\odot}$, respectively. Next, we compute the time-dependent luminosity of key FIR lines tracing the neutral and ionized gas layers of the GMCs, ([CII] at $158\,\rm{\mu m}$, [OIII] at $88\,\rm \mu m$) as a function of $G_0$, and $Z$ until complete photoevaporation at $t_{pe}$. We find that the specific [CII] luminosity is almost independent on the GMC model within the survival time of the cloud. Stronger FUV fluxes produce higher [CII] and [OIII] luminosities, however lasting for progressively shorter times. At $Z=Z_{\odot}$ the [CII] emission is maximized ($L_{\rm CII}\approx 10^4\,\rm{L_{\odot}}$ for the fiducial model) for $t<1\,\rm{Myr}$ and $\log G_0\geq 3$. Noticeably, and consistently with the recent detection by Inoue et al. (2016) of a galaxy at redshift $z\approx 7.2$, for $Z\leq 0.2\,\rm{Z_{\odot}}$ the [OIII] line might outshine [CII] emission by up to $\approx 1000$ times. We conclude that the [OIII] line is a key diagnostic of low metallicity ISM, especially in galaxies with very young stellar populations.

PROSPECTS FOR CHEMICALLY TAGGING STARS IN THE GALAXY

It is now well-established that the elemental abundance patterns of stars hold key clues not only to their formation, but also to the assembly histories of galaxies. One of the most exciting possibilities is the use of stellar abundance patterns as "chemical tags" to identify stars that were born in the same molecular cloud. In this paper, we assess the prospects of chemical tagging as a function of several key underlying parameters. We show that in the fiducial case of 104 distinct cells in chemical space and stars in the survey, one can expect to detect groups that are overdensities in the chemical space. However, we find that even very large overdensities in chemical space do not guarantee that the overdensity is due to a single set of stars from a common birth cloud. In fact, for our fiducial model parameters, the typical overdensity is comprised of stars from a wide range of clusters with the most dominant cluster contributing only 25% of the stars. The most important factors limiting the identification of disrupted clusters via chemical tagging are the number of chemical cells in the chemical space and the survey sampling rate of the underlying stellar population. Both of these factors can be improved through strategic observational plans. While recovering individual clusters through chemical tagging may prove challenging, we show, in agreement with previous work, that different CMFs imprint different degrees of clumpiness in chemical space. These differences provide the opportunity to statistically reconstruct the slope and high-mass cutoff of CMF and its evolution through cosmic time.

THE ZURICH ENVIRONMENTAL STUDY (ZENS) OF GALAXIES IN GROUPS ALONG THE COSMIC WEB. V. PROPERTIES AND FREQUENCY OF MERGING SATELLITES AND CENTRALS IN DIFFERENT ENVIRONMENTS

We use the Zurich ENvironmental Study (ZENS) database to investigate the environmental dependence of the merger fraction $\Gamma$ and merging galaxy properties in a sample of ~1300 group galaxies with $M>10^{9.2}M_\odot$ and 0.05<z<0.0585. In all galaxy mass bins investigated in our study, we find that $\Gamma$ decreases by a factor of ~2-3 in groups with halo masses $M_{HALO}>10^{13.5} M_\odot$ relative to less massive systems, indicating a suppression of merger activity in large potential wells. In the fiducial case of relaxed groups only, we measure a variation $\Delta\Gamma/\Delta \log (M_{HALO}) \sim - 0.07$ dex$^{-1}$, which is almost independent of galaxy mass and merger stage. At galaxy masses $>10^{10.2} M_\odot$, most mergers are dry accretions of quenched satellites onto quenched centrals, leading to a strong increase of $\Gamma$ with decreasing group-centric distance at these mass scales.Both satellite and central galaxies in these high mass mergers do not differ in color and structural properties from a control sample of nonmerging galaxies of equal mass and rank. At galaxy masses $<10^{10.2} M_\odot$, where we mostly probe satellite-satellite pairs and mergers between star-forming systems, close pairs (projected distance $<10-20$ kpc) show instead $\sim2\times$ enhanced (specific) star formation rates and $\sim1.5\times$ larger sizes than similar mass, nonmerging satellites. The increase in both size and SFR leads to similar surface star-formation densities in the merging and control-sample satellite populations.

Sensitivity of CTA to dark matter signals from the Galactic Center

The Galactic Center is one of the most promising targets for indirect detection of dark matter with gamma rays. We investigate the sensitivity of the upcoming Cherenkov Telescope Array (CTA) to dark matter annihilation and decay in the Galactic Center. As the inner density profile of the Milky Way's dark matter halo is uncertain, we study the impact of the slope of the Galactic density profile, inwards of the Sun, on the prospects for detecting a dark matter signal with CTA. Adopting the Ring Method to define the signal and background regions in an ON-OFF analysis approach, we find that the sensitivity achieved by CTA to annihilation signals is strongly dependent on the inner profile slope, whereas the dependence is more mild in the case of dark matter decay. Surprisingly, we find that the optimal choice of signal and background regions is virtually independent of the assumed density profile. For the fiducial case of a Navarro-Frenk-White profile,we find that CTA will be able to probe annihilation cross-sections well below the canonical thermal relic value for dark matter masses from a few tens of GeV up to ∼ 5 TeV for annihilation to τ+τ−, and will achieve only a slightly weaker sensitivity for annihilation to bb̄ or μ+μ−. CTA will improve significantly on current sensitivity to annihilation signals for dark matter masses above ∼ 100 GeV, covering parameter space that is complementary to that probed by searches with the Fermi Large Area Telescope. The interpretation of apparent excesses in the measured cosmic-ray electron and positron spectra as signals of dark matter decay will also be testable with CTA observations of the Galactic Center. We demonstrate that both for annihilation and for decay, including spectral information for hard channels (such as μ+μ− and τ+τ−) leads to enhanced sensitivity for dark matter masses above mDM ∼ 200 GeV.

THE STRUCTURE OF SURFACE H2O LAYERS OF ICE-COVERED PLANETS WITH HIGH-PRESSURE ICE

Many extrasolar (bound) terrestrial planets and free-floating (unbound) planets have been discovered. The existence of bound and unbound terrestrial planets with liquid water is an important question, and of particular importance is the question of their habitability. Even for a globally ice-covered planet, geothermal heat from the planetary interior may melt the interior ice, creating an internal ocean covered by an ice shell. In this paper, we discuss the conditions that terrestrial planets must satisfy for such an internal ocean to exist on the timescale of planetary evolution. The question is addressed in terms of planetary mass, distance from a central star, water abundance, and abundance of radiogenic heat sources. In addition, we investigate the structures of the surface-H2O layers of ice-covered planets by considering the effects of ice under high pressure (high-pressure ice). As a fiducial case, 1M$\oplus$ planet at 1 AU from its central star and with 0.6 to 25 times the H2O mass of Earth could have an internal ocean. We find that high-pressure ice layers may appear between the internal ocean and the rock portion on a planet with an H2O mass over 25 times that of Earth. The planetary mass and abundance of surface water strongly restrict the conditions under which an extrasolar terrestrial planet may have an internal ocean with no high-pressure ice under the ocean. Such high-pressure-ice layers underlying the internal ocean are likely to affect the habitability of the planet.

THERMAL AND CHEMICAL EVOLUTION OF COLLAPSING FILAMENTS

Intergalactic filaments form the foundation of the cosmic web that connect galaxies together, and provide an important reservoir of gas for galaxy growth and accretion. Here we present very high resolution two-dimensional simulations of the thermal and chemical evolution of such filaments, making use of a 32 species chemistry network that tracks the evolution of key molecules formed from hydrogen, oxygen, and carbon. We study the evolution of filaments over a wide range of parameters including the initial density, initial temperature, strength of the dissociating UV background, and metallicity. In low-redshift, $Z \approx 0.1 Z_\odot $ filaments, the evolution is determined completely by the initial cooling time. If this is sufficiently short, the center of the filament always collapses to form dense, cold core containing a substantial fraction of molecules. In high-redshift, $Z=10^{-3} Z_\odot$ filaments, the collapse proceeds much more slowly. This is due mostly to the lower initial temperatures, which leads to a much more modest increase in density before the atomic cooling limit is reached, making subsequent molecular cooling much less efficient. Finally, we study how the gravitational potential from a nearby dwarf galaxy affects the collapse of the filament and compare this to NGC 5253, a nearby starbusting dwarf galaxy thought to be fueled by the accretion of filament gas. In contrast to our fiducial case, a substantial density peak forms at the center of the potential. This peak evolves faster than the rest of the filament due to the increased rate at which chemical species form and cooling occur. We find that we achieve similar accretion rates as NGC 5253 but our two-dimensional simulations do not recover the formation of the giant molecular clouds that are seen in radio observations.

PROTOSTELLAR FEEDBACK AND FINAL MASS OF THE SECOND-GENERATION PRIMORDIAL STARS

The first stars in the universe ionized the ambient primordial gas through various feedback processes. "Second-generation" primordial stars potentially form from this disturbed gas after its recombination. In this Letter, we study the late formation stage of such second-generation stars, where a large amount of gas accretes onto the protostar and the final stellar mass is determined when the accretion terminates. We directly compute the complex interplay between the accretion flow and stellar ultraviolet (UV) radiation, performing radiation-hydrodynamic simulations coupled with stellar evolution calculations. Because of more efficient H2 and HD cooling in the pre-stellar stage, the accretion rates onto the star are ten times lower than in the case of the formation of the first stars. The lower accretion rates and envelope density result in the occurrence of an expanding bipolar HII region at a lower protostellar mass M_* \simeq 10Msun, which blows out the circumstellar material, thereby quenching the mass supply from the envelope to the accretion disk. At the same time the disk loses mass due to photoevaporation by the growing star. In our fiducial case the stellar UV feedback terminates mass accretion onto the star at M_* \simeq 17Msun. Although the derived masses of the second-generation primordial stars are systematically lower than those of the first generation, the difference is within a factor of only a few. Our results suggest a new scenario, whereby the majority of the primordial stars are born as massive stars with tens of solar masses, regardless of their generations.

VISCOUS-LIKE INTERACTION OF THE SOLAR WIND WITH THE PLASMA TAIL OF COMET SWIFT-TUTTLE

We compare the results of the numerical simulation of the viscous-like interaction of the solar wind with the plasma tail of a comet, with velocities of H2O+ ions in the tail of comet Swift-Tuttle determined by means of spectroscopic, ground based observations. Our aim is to constrain the value of the basic parameters in the viscous-like interaction model: the effective Reynolds number of the flow and the interspecies coupling timescale. We find that in our simulations the flow rapidly evolves from an arbitrary initial condition to a quasi-steady state for which there is a good agreement between the simulated tailward velocity of H2O+ ions and the kinematics derived from the observations. The fiducial case of our model, characterized by a low effective Reynolds number (Re \approx 20 and selected on the basis of a comparison to in situ measurements of the plasma flow at comet Halley, yields an excellent fit to the observed kinematics. Given the agreement between model and observations, with no ad hoc assumptions, we believe that this result suggests that viscous-like momentum transport may play an important role in the interaction of the solar wind and the cometary plasma environment.

Impact of early dark energy on the Planck SZ cluster sample

Context. One science goal of the upcoming Planck mission is to perform a full-sky cluster survey based on the Sunyaev-Zel'dovich (SZ) effect, which leads to the question of how such a survey would be affected by cosmological models with a different history of structure formation than LCDM. One class of these models are early dark energy (EDE) cosmologies, where the dark energy contribution does not vanish at early times. Aims. Since structures grow slower in the presence of EDE, one expects an increase in the number of galaxy clusters compared to LCDM at intermediate and high redshifts, which could explain the reported excess of the angular CMB power spectrum on cluster scales via an enhanced SZ contribution. We study the impact of EDE on Planck's expected cluster sample. Methods. To obtain realistic simulations, we constructed full-sky SZ maps for EDE and LCDM cosmologies, taking angular cluster correlation into account. Using these maps, we simulated Planck observations with and without Galactic foregrounds and fed the results into our filter pipeline based on the spherical multi-frequency matched filters. Results. For the case of EDE cosmologies, we clearly find an increase in the detected number of clusters compared to the fiducial LCDM case. This shows that the spherical multi-frequency matched filter is sensitive enough to find deviations from the LCDM sample, being caused by EDE. In addition we find an interesting effect of EDE on the completeness of the cluster sample, such that EDE helps to obtain cleaner samples.

The Formation of the First Stars. II. Radiative Feedback Processes and Implications for the Initial Mass Function

We consider the radiative feedback processes that operate during the formation of the first stars, including the photodissociation of H_2, Ly-alpha radiation pressure, formation and expansion of an HII region, and disk photoevaporation. These processes may inhibit continued accretion once the stellar mass has reached a critical value, and we evaluate this mass separately for each process. Photodissociation of H_2 in the local dark matter minihalo occurs relatively early in the growth of the protostar, but we argue this does not affect subsequent accretion since by this time the depth of the potential is large enough for accretion to be mediated by atomic cooling. However, neighboring starless minihalos can be affected. Ionization creates an HII region in the infalling envelope above and below the accretion disk. Ly-alpha radiation pressure acting at the boundary of the HII region is effective at reversing infall from narrow polar directions when the star reaches ~20-30Msun, but cannot prevent infall from other directions. Expansion of the HII region beyond the gravitational escape radius for ionized gas occurs at masses ~50-100Msun, depending on the accretion rate and angular momentum of the inflow. However, again, accretion from the equatorial regions can continue since the neutral accretion disk has a finite thickness and shields a substantial fraction of the accretion envelope from direct ionizing flux. At higher stellar masses, ~140Msun in the fiducial case, the combination of declining accretion rates and increasing photoevaporation-driven mass loss from the disk act to effectively halt the increase in the protostellar mass. We identify this process as the mechanism that terminates the growth of Population III stars... (abridged)

Solar winds driven by nonlinear low‐frequency Alfvén waves from the photosphere: Parametric study for fast/slow winds and disappearance of solar winds

We investigate how properties of the corona and solar wind in open coronal holes depend on properties of magnetic fields and their footpoint motions at the surface. We perform one‐dimensional magnetohydrodynamical (MHD) simulations for the heating and the acceleration in coronal holes by low‐frequency Alfvén waves from the photosphere to 0.3 or 0.1 AU. We impose low‐frequency (≲0.05 Hz) transverse fluctuations of the field lines at the photosphere with various amplitude, spectrum, and polarization in the open flux tubes with different photospheric field strength, Br,0, and superradial expansion of the cross section, fmax. We find that transonic solar winds are universal consequences. The atmosphere is also stably heated up to ≳106 K by the dissipation of the Alfvén waves through compressive‐wave generation and wave reflection in the cases of the sufficient wave input with photospheric amplitude, 〈dv⊥,0〉 ≳ 0.7 km s−1. The density, and accordingly the mass flux, of solar winds show a quite sensitive dependence on 〈dv⊥,0〉 because of an unstable aspect of the heating by the nonlinear Alfvén waves. A case with 〈dv⊥,0〉 = 0.4 km s−1 gives ≃50 times smaller mass flux than the fiducial case for the fast wind with 〈dv⊥,0〉 = 0.7 km s−1; solar wind virtually disappears only if 〈dv⊥,0〉 becomes ≃1/2. We also find that the solar wind speed has a positive correlation with Br,0/fmax, which is consistent with recent observations by Kojima et al. On the basis of these findings, we show that both fast and slow solar winds can be explained by the single process, the dissipation of the low‐frequency Alfvén waves, with different sets of 〈dv⊥,0〉 and Br,0/fmax. Our simulations naturally explain the observed (1) anticorrelation of the solar wind speed and the coronal temperature and (2) larger amplitude of Alfvénic fluctuations in the fast wind. In Appendix A, we also explain our implementation of the outgoing boundary condition of the MHD waves with some numerical tests.

The Formation of the First Stars. I. Mass Infall Rates, Accretion Disk Structure, and Protostellar Evolution

We present a theoretical model for primordial star formation. First we describe the structure of the initial gas cores as virialized, quasi-hydrostatic objects in accord with recent high-resolution numerical studies. The accretion rate can then be related to characteristic densities and temperatures that are set by the cooling properties of molecular hydrogen. We allow for rotation of the gas core, assuming angular momentum conservation inside the sonic point of the flow. In the typical case, most mass then reaches the star via an accretion disk. The structure of the inner region of this disk is described with the standard theory of viscous disks, but with allowance for the substantial energies absorbed in ionizing and dissociating the gas. The size of the protostar and its luminosity depend on the accretion rate, the energetics of the accreting gas, and the ability of the radiation to escape from the stellar accretion shock. We combine these models for the infall rate, inner disk structure, and protostellar evolution to predict the radiation field that is the basis for radiative feedback processes acting against infall (second paper in the series). For realistic initial angular momenta, the photosphere of the protostar is much smaller and hotter than in the spherical case, leading to stronger radiative feedback at earlier stages in the evolution. In particular, once the star is older than its Kelvin-Helmholtz time, contraction toward the main sequence causes a rapid increase in ionizing and far-ultraviolet luminosity at masses ~30 M? in the fiducial case. Since the cores out of which the first stars formed were much more massive than 30 M? and since feedback is dynamically unimportant at lower masses, we conclude that the first stars should have had masses 30 M?.

Orbital dynamics of three-dimensional bars - I. The backbone of three-dimensional bars. A fiducial case

In this series of papers we investigate the orbital structure of three-dimensional (3D) models representing barred galaxies. In the present introductory paper we use a fiducial case to describe all families of periodic orbits that may play a role in the morphology of three-dimensional bars. We show that, in a 3D bar, the backbone of the orbital structure is not just the x1 family, as in two-dimensional (2D) models, but a tree of 2D and 3D families bifurcating from x1. Besides the main tree we have also found another group of families of lesser importance around the radial 3:1 resonance. The families of this group bifurcate from x1 and influence the dynamics of the system only locally. We also find that 3D orbits elongated along the bar minor axis can be formed by bifurcations of the planar x2 family. They can support 3D bar-like structures along the minor axis of the main bar. Banana-like orbits around the stable Lagrangian points build a forest of 2D and 3D families as well. The importance of the 3D x1-tree families at the outer parts of the bar depends critically on whether they are introduced in the system as bifurcations in z or in z˙.

Dynamical Chaos in the Wisdom-Holman Integrator: Origins and Solutions

We examine the nonlinear stability of the Wisdom-Holman (WH) symplectic mapping applied to the integration of perturbed, highly eccentric (e ≳ 0.9) two-body orbits. We find that the method is unstable and introduces artificial chaos into the computed trajectories for this class of problems, unless the step size chosen is small enough that periapse is always resolved, in which case the method is generically stable. This "radial orbit instability" persists even for weakly perturbed systems. Using the Stark problem as a fiducial test case, we investigate the dynamical origin of this instability and argue that the numerical chaos results from the overlap of step-size resonances; interestingly, for the Stark problem many of these resonances appear to be absolutely stable.