Galaxies in box: A simulated view of the interstellar medium

We review progress in the development of physically realistic three dimensional simulated models of the galaxy.We consider the scales from star forming molecular clouds to the full spiral disc. Models are computed using hydrodynamic (HD) or magnetohydrodynamic (MHD) equations and may include cosmic ray or tracer particles. The range of dynamical scales between the full galaxy structure and the turbulent scales of supernova (SN) explosions and even cloud collapse to form stars, make it impossible with current computing tools and resources to resolve all of these in one model. We therefore consider a hierarchy of models and how they can be related to enhance our understanding of the complete galaxy.Comment: Chapter in Large Scale Magnetic Fields in the Univers

Simulating Turbulence Using the Astrophysical Discontinuous Galerkin Code TENET

In astrophysics, the two main methods traditionally in use for solving the Euler equations of ideal fluid dynamics are smoothed particle hydrodynamics and finite volume discretization on a stationary mesh. However, the goal to efficiently make use of future exascale machines with their ever higher degree of parallel concurrency motivates the search for more efficient and more accurate techniques for computing hydrodynamics. Discontinuous Galerkin (DG) methods represent a promising class of methods in this regard, as they can be straightforwardly extended to arbitrarily high order while requiring only small stencils. Especially for applications involving comparatively smooth problems, higher-order approaches promise significant gains in computational speed for reaching a desired target accuracy. Here, we introduce our new astrophysical DG code TENET designed for applications in cosmology, and discuss our first results for 3D simulations of subsonic turbulence. We show that our new DG implementation provides accurate results for subsonic turbulence, at considerably reduced computational cost compared with traditional finite volume methods. In particular, we find that DG needs about 1.8 times fewer degrees of freedom to achieve the same accuracy and at the same time is more than 1.5 times faster, confirming its substantial promise for astrophysical applications.Comment: 21 pages, 7 figures, to appear in Proceedings of the SPPEXA symposium, Lecture Notes in Computational Science and Engineering (LNCSE), Springe

The Science of Sungrazers, Sunskirters, and Other Near-Sun Comets

This review addresses our current understanding of comets that venture close to the Sun, and are hence exposed to much more extreme conditions than comets that are typically studied from Earth. The extreme solar heating and plasma environments that these objects encounter change many aspects of their behaviour, thus yielding valuable information on both the comets themselves that complements other data we have on primitive solar system bodies, as well as on the near-solar environment which they traverse. We propose clear definitions for these comets: We use the term near-Sun comets to encompass all objects that pass sunward of the perihelion distance of planet Mercury (0.307 AU). Sunskirters are defined as objects that pass within 33 solar radii of the Sun’s centre, equal to half of Mercury’s perihelion distance, and the commonly-used phrase sungrazers to be objects that reach perihelion within 3.45 solar radii, i.e. the fluid Roche limit. Finally, comets with orbits that intersect the solar photosphere are termed sundivers. We summarize past studies of these objects, as well as the instruments and facilities used to study them, including space-based platforms that have led to a recent revolution in the quantity and quality of relevant observations. Relevant comet populations are described, including the Kreutz, Marsden, Kracht, and Meyer groups, near-Sun asteroids, and a brief discussion of their origins. The importance of light curves and the clues they provide on cometary composition are emphasized, together with what information has been gleaned about nucleus parameters, including the sizes and masses of objects and their families, and their tensile strengths. The physical processes occurring at these objects are considered in some detail, including the disruption of nuclei, sublimation, and ionisation, and we consider the mass, momentum, and energy loss of comets in the corona and those that venture to lower altitudes. The different components of comae and tails are described, including dust, neutral and ionised gases, their chemical reactions, and their contributions to the near-Sun environment. Comet-solar wind interactions are discussed, including the use of comets as probes of solar wind and coronal conditions in their vicinities. We address the relevance of work on comets near the Sun to similar objects orbiting other stars, and conclude with a discussion of future directions for the field and the planned ground- and space-based facilities that will allow us to address those science topics

Physical Processes in Star Formation

© 2020 Springer-Verlag. The final publication is available at Springer via https://doi.org/10.1007/s11214-020-00693-8.Star formation is a complex multi-scale phenomenon that is of significant importance for astrophysics in general. Stars and star formation are key pillars in observational astronomy from local star forming regions in the Milky Way up to high-redshift galaxies. From a theoretical perspective, star formation and feedback processes (radiation, winds, and supernovae) play a pivotal role in advancing our understanding of the physical processes at work, both individually and of their interactions. In this review we will give an overview of the main processes that are important for the understanding of star formation. We start with an observationally motivated view on star formation from a global perspective and outline the general paradigm of the life-cycle of molecular clouds, in which star formation is the key process to close the cycle. After that we focus on the thermal and chemical aspects in star forming regions, discuss turbulence and magnetic fields as well as gravitational forces. Finally, we review the most important stellar feedback mechanisms.Peer reviewedFinal Accepted Versio

Dense Stellar Populations: Initial Conditions

This chapter is based on four lectures given at the Cambridge N-body school "Cambody". The material covered includes the IMF, the 6D structure of dense clusters, residual gas expulsion and the initial binary population. It is aimed at those needing to initialise stellar populations for a variety of purposes (N-body experiments, stellar population synthesis).Comment: 85 pages. To appear in The Cambridge N-body Lectures, Sverre Aarseth, Christopher Tout, Rosemary Mardling (eds), Lecture Notes in Physics Series, Springer Verla

The Physics of Star Cluster Formation and Evolution

© 2020 Springer-Verlag. The final publication is available at Springer via https://doi.org/10.1007/s11214-020-00689-4.Star clusters form in dense, hierarchically collapsing gas clouds. Bulk kinetic energy is transformed to turbulence with stars forming from cores fed by filaments. In the most compact regions, stellar feedback is least effective in removing the gas and stars may form very efficiently. These are also the regions where, in high-mass clusters, ejecta from some kind of high-mass stars are effectively captured during the formation phase of some of the low mass stars and effectively channeled into the latter to form multiple populations. Star formation epochs in star clusters are generally set by gas flows that determine the abundance of gas in the cluster. We argue that there is likely only one star formation epoch after which clusters remain essentially clear of gas by cluster winds. Collisional dynamics is important in this phase leading to core collapse, expansion and eventual dispersion of every cluster. We review recent developments in the field with a focus on theoretical work.Peer reviewe

Dust in Supernovae and Supernova Remnants II: Processing and survival

Observations have recently shown that supernovae are efficient dust factories, as predicted for a long time by theoretical models. The rapid evolution of their stellar progenitors combined with their efficiency in precipitating refractory elements from the gas phase into dust grains make supernovae the major potential suppliers of dust in the early Universe, where more conventional sources like Asymptotic Giant Branch (AGB) stars did not have time to evolve. However, dust yields inferred from observations of young supernovae or derived from models do not reflect the net amount of supernova-condensed dust able to be expelled from the remnants and reach the interstellar medium. The cavity where the dust is formed and initially resides is crossed by the high velocity reverse shock which is generated by the pressure of the circumstellar material shocked by the expanding supernova blast wave. Depending on grain composition and initial size, processing by the reverse shock may lead to substantial dust erosion and even complete destruction. The goal of this review is to present the state of the art about processing and survival of dust inside supernova remnants, in terms of theoretical modelling and comparison to observations

A MODEST review

We present an account of the state of the art in the fields explored by the research community invested in 'Modeling and Observing DEnse STellar systems'. For this purpose, we take as a basis the activities of the MODEST-17 conference, which was held at Charles University, Prague, in September 2017. Reviewed topics include recent advances in fundamental stellar dynamics, numerical methods for the solution of the gravitational N-body problem, formation and evolution of young and old star clusters and galactic nuclei, their elusive stellar populations, planetary systems, and exotic compact objects, with timely attention to black holes of different classes of mass and their role as sources of gravitational waves. Such a breadth of topics reflects the growing role played by collisional stellar dynamics in numerous areas of modern astrophysics. Indeed, in the next decade, many revolutionary instruments will enable the derivation of positions and velocities of individual stars in the Milky Way and its satellites and will detect signals from a range of astrophysical sources in different portions of the electromagnetic and gravitational spectrum, with an unprecedented sensitivity. On the one hand, this wealth of data will allow us to address a number of long-standing open questions in star cluster studies; on the other hand, many unexpected properties of these systems will come to light, stimulating further progress of our understanding of their formation and evolution.Comment: 42 pages; accepted for publication in 'Computational Astrophysics and Cosmology'. We are much grateful to the organisers of the MODEST-17 conference (Charles University, Prague, September 2017). We acknowledge the input provided by all MODEST-17 participants, and, more generally, by the members of the MODEST communit

Whole-genome sequencing reveals host factors underlying critical COVID-19

Critical COVID-19 is caused by immune-mediated inflammatory lung injury. Host genetic variation influences the development of illness requiring critical care1 or hospitalization2,3,4 after infection with SARS-CoV-2. The GenOMICC (Genetics of Mortality in Critical Care) study enables the comparison of genomes from individuals who are critically ill with those of population controls to find underlying disease mechanisms. Here we use whole-genome sequencing in 7,491 critically ill individuals compared with 48,400 controls to discover and replicate 23 independent variants that significantly predispose to critical COVID-19. We identify 16 new independent associations, including variants within genes that are involved in interferon signalling (IL10RB and PLSCR1), leucocyte differentiation (BCL11A) and blood-type antigen secretor status (FUT2). Using transcriptome-wide association and colocalization to infer the effect of gene expression on disease severity, we find evidence that implicates multiple genes—including reduced expression of a membrane flippase (ATP11A), and increased expression of a mucin (MUC1)—in critical disease. Mendelian randomization provides evidence in support of causal roles for myeloid cell adhesion molecules (SELE, ICAM5 and CD209) and the coagulation factor F8, all of which are potentially druggable targets. Our results are broadly consistent with a multi-component model of COVID-19 pathophysiology, in which at least two distinct mechanisms can predispose to life-threatening disease: failure to control viral replication; or an enhanced tendency towards pulmonary inflammation and intravascular coagulation. We show that comparison between cases of critical illness and population controls is highly efficient for the detection of therapeutically relevant mechanisms of disease