Galactic fountain flows in high-resolution hydrodynamical simulations of galaxy disks

Feedback from the most massive components of a young stellar cluster deeply affects the surrounding ISM driving an expanding over-pressured hot gas cavity in it. In spiral galaxies these structures may have sufficient energy to break the disk and eject large amount of material into the halo. The cycling of this gas, which eventually will fall back onto the disk, is known as galactic fountains. We aim at better understanding the dynamics of such fountain flow in a Galactic context, frame the problem in a more dynamic environment possibly learning about its connection and regulation to the local driving mechanism and understand its role as a metal diffusion channel. The interaction of the fountain with a hot corona is hereby analyzed, trying to understand the properties and evolution of the extraplanar material. We perform high resolution hydrodynamical simulations with the moving-mesh code AREPO to model the multi-phase ISM of a Milky Way type galaxy. A non-equilibrium chemical network is included to self consistently follow the evolution of the main coolants of the ISM. Spiral arm perturbations in the potential are considered so that large molecular gas structures are able to dynamically form here, self shielded from the interstellar radiation field. We model the effect of SN feedback from a new-born stellar cluster inside such a giant molecular cloud, as the driving force of the fountain. Passive Lagrangian tracer particles are used in conjunction to the SN energy deposition to model and study diffusion of freshly synthesized metals. We find that both interactions with hot coronal gas and local ISM properties and motions are equally important in shaping the fountain. We notice a bimodal morphology where most of the ejected gas is in a cold $10^4$ K clumpy state while the majority of the affected volume is occupied by a hot diffuse medium. While only about 20\% of the produced metals stay local, most of them quickly diffuse through this hot regime to great scales.

Optical evidence for the unification of active galactic nuclei and quasi-stellar objects.

There is a variety of optical evidence for some unification of different types of active galactic nuclei and quasi-stellar objects (QSOs). The case is very strong for the unification of at least some Seyfert galaxies, where polarization data show that the type assigned to the Seyfert galaxy must depend on viewing direction. It has been proposed that Fanaroff-Riley type 2 (FR2) radio galaxies are quasars seen in a direction from which the quasar is obscured, and there is some limited direct evidence for this picture. The broad absorption line QSOs may be normal QSOs seen from a special direction. Some of the sources observed to have high luminosities in the far infrared could be obscured QSOs and active nuclei. Mergers and interactions are likely to play an important role in nuclear activity, and active galaxies and QSOs could change their apparent types through these encounters followed by subsequent evolution.

A VLT/FLAMES survey for massive binaries in Westerlund 1. III. The WC9d binary W239 and implications for massive stellar evolution

Context. There is growing evidence that a treatment of binarity amongst OB stars is essential for a full theory of stellar evolution. However the binary properties of massive stars – frequency, mass ratio & orbital separation – are still poorly constrained. Aims. In order to address this shortcoming we have undertaken a multiepoch spectroscopic study of the stellar population of the young massive cluster Westerlund 1. In this paper we present an investigation into the nature of the dusty Wolf-Rayet star and candidate binary W239. Methods. To accomplish this we have utilised our spectroscopic data in conjunction with multi-year optical and near-IR photometric observations in order to search for binary signatures. Comparison of these data to synthetic non-LTE model atmosphere spectra were used to derive the fundamental properties of the WC9 primary. Results. We found W239 to have an orbital period of only ~5.05 days, making it one of the most compact WC binaries yet identified. Analysis of the long term near-IR lightcurve reveals a significant flare between 2004-6. We interpret this as evidence for a third massive stellar component in the system in a long period (>6 yr), eccentric orbit, with dust production occuring at periastron leading to the flare. The presence of a near-IR excess characteristic of hot (~1300 K) dust at every epoch is consistent with the expectation that the subset of persistent dust forming WC stars are short (<1 yr) period binaries, although confirmation will require further observations. Non-LTE model atmosphere analysis of the spectrum reveals the physical properties of the WC9 component to be fully consistent with other Galactic examples. Conclusions. The simultaneous presence of both short period Wolf-Rayet binaries and cool hypergiants within Wd 1 provides compelling evidence for a bifurcation in the post-Main Sequence evolution of massive stars due to binarity. Short period O+OB binaries will evolve directly to the Wolf-Rayet phase, either due to an episode of binary mediated mass loss – likely via case A mass transfer or a contact configuration – or via chemically homogenous evolution. Conversely, long period binaries and single stars will instead undergo a red loop across the HR diagram via a cool hypergiant phase. Future analysis of the full spectroscopic dataset for Wd 1 will constrain the proportion of massive stars experiencing each pathway; hence quantifying the importance of binarity in massive stellar evolution up to and beyond supernova and the resultant production of relativistic remnants.

On the nature of the galactic early-B hypergiants

Aims. Despite their importance to a number of astrophysical fields, the lifecycles of very massive stars are still poorly defined. In order to address this shortcoming, we present a detailed quantitative study of the physical properties of four early-B hypergiants (BHGs) of spectral type B1-4 Ia+; Cyg OB2 #12, ζ1 Sco, HD 190603 and BP Cru. These are combined with an analysis of their long-term spectroscopic and photometric behaviour in order to determine their evolutionary status. Methods. Quantitative analysis of UV–radio photometric and spectroscopic datasets was undertaken with a non-LTE model atmosphere code in order to derive physical parameters for comparison with apparently closely related objects, such as B supergiants (BSGs) and luminous blue variables (LBVs), and theoretical evolutionary predictions. Results. The long-term photospheric and spectroscopic datasets compiled for the early-B HGs revealed that they are remarkably stable over long periods ( ≥ 40 yrs), with the possible exception of ζ1 Sco prior to the 20th century; in contrast to the typical excursions that characterise LBVs. Quantitative analysis of ζ1 Sco, HD 190603 and BP Cru yielded physical properties intermediate between BSGs and LBVs; we therefore suggest that BHGs are the immediate descendants and progenitors (respectively) of such stars, for initial masses in the range ~30−60 M⊙. Comparison of the properties of ζ1 Sco with the stellar population of its host cluster/association NGC 6231/Sco OB1 provides further support for such an evolutionary scenario. In contrast, while the wind properties of Cyg OB2 #12 are consistent with this hypothesis, the combination of extreme luminosity and spectroscopic mass (~110 M⊙) and comparatively low temperature means it cannot be accommodated in such a scheme. Likewise, despite its co-location with several LBVs above the Humphreys-Davidson (HD) limit, the lack of long term variability and its unevolved chemistry apparently excludes such an identification. Since such massive stars are not expected to evolve to such cool temperatures, instead traversing an O4-6Ia → O4-6Ia+ → WN7-9ha pathway, the properties of Cyg OB2 #12 are therefore difficult to understand under current evolutionary paradigms. Finally, we note that as with AG Car in its cool phase, despite exceeding the HD limit, the properties of Cyg OB2 #12 imply that it lies below the Eddington limit – thus we conclude that the HD limit does not define a region of the HR diagram inherently inimical to the presence of massive stars.

IRAS 18357-0604 – an analogue of the galactic yellow hypergiant IRC +10420?

Context. Yellow hypergiants represent a short-lived evolutionary episode experienced by massive stars as they transit to and from a red supergiant phase. As such, their properties provide a critical test of stellar evolutionary theory, while recent observations unexpectedly suggest that a subset may explode as Type II supernovae. Aims. The galactic yellow hypergiant IRC +10420 is a cornerstone system for understanding this phase since it is the strongest post-RSG candidate known, has demonstrated real-time evolution across the Hertzsprung-Russell diagram and been subject to extensive mass loss. In this paper we report on the discovery of a twin of IRC +10420 - IRAS 18357-0604. Methods. Optical and near-IR spectroscopy are used to investigate the physical properties of IRAS 18357-0604 and also provide an estimate of its systemic velocity, while near- to mid-IR photometry probes the nature of its circumstellar environment. Results. These observations reveal pronounced spectral similarities between IRAS 18357-0604 and IRC +10420, suggesting comparable temperatures and wind geometries. IR photometric data reveals a similarly dusty circumstellar environment, although historical mass loss appears to have been heavier in IRC +10420. The systemic velocity implies a distance compatible with the red supergiant-dominated complex at the base of the Scutum Crux arm; the resultant luminosity determination is consistent with a physical association but suggests a lower initial mass than inferred for IRC +10420 (≲20 M⊙ versus ~40 M⊙). Evolutionary predictions for the physical properties of supernova progenitors derived from ~18–20 M⊙ stars – or ~12–15 M⊙ stars that have experienced enhanced mass loss as red supergiants – compare favourably with those of IRAS 18357-0604, which in turn appears to be similar to the the progenitor of SN2011dh; it may therefore provide an important insight into the nature of the apparently H-depleted yellow hypergiant progenitors of some Type IIb SNe.

The X-ray outburst of the Galactic Centre magnetar SGR J1745−2900 during the first 1.5 year

In 2013 April a new magnetar, SGR 1745−2900, was discovered as it entered an outburst, at only 2.4 arcsec angular distance from the supermassive black hole at the centre of the Milky Way, Sagittarius A*. SGR 1745−2900 has a surface dipolar magnetic field of ∼2 × 1014 G, and it is the neutron star closest to a black hole ever observed. The new source was detected both in the radio and X-ray bands, with a peak X-ray luminosity LX ∼ 5 × 1035 erg s−1. Here we report on the long-term Chandra (25 observations) and XMM–Newton (eight observations) X-ray monitoring campaign of SGR 1745−2900 from the onset of the outburst in 2013 April until 2014 September. This unprecedented data set allows us to refine the timing properties of the source, as well as to study the outburst spectral evolution as a function of time and rotational phase. Our timing analysis confirms the increase in the spin period derivative by a factor of ∼2 around 2013 June, and reveals that a further increase occurred between 2013 October 30 and 2014 February 21. We find that the period derivative changed from 6.6 × 10−12 to 3.3 × 10−11 s s−1 in 1.5 yr. On the other hand, this magnetar shows a slow flux decay compared to other magnetars and a rather inefficient surface cooling. In particular, starquake-induced crustal cooling models alone have difficulty in explaining the high luminosity of the source for the first ∼200 d of its outburst, and additional heating of the star surface from currents flowing in a twisted magnetic bundle is probably playing an important role in the outburst evolution.

On the formation and evolution of the first Be star in a black hole binary MWC 656

We find that the formation of MWC 656 (the first Be binary containing a black hole) involves a common envelope phase and a supernova explosion. This result supports the idea that a rapidly rotating Be star can emerge out of a common envelope phase, which is very intriguing because this evolutionary stage is thought to be too fast to lead to significant accretion and spin up of the B star. We predict ∼10–100 of B-BH binaries to currently reside in the Galactic disc, among which around 1/3 contain a Be star, but there is only a small chance to observe a system with parameters resembling MWC 656. If MWC 656 is representative of intrinsic Galactic Be-BH binary population, it may indicate that standard evolutionary theory needs to be revised. This would pose another evolutionary problem in understanding black hole (BH) binaries, with BH X-ray novae formation issue being the prime example. Future evolution of MWC 656 with an ∼5 M⊙ BH and with an ∼13 M⊙ main-sequence companion on an ∼60 d orbit may lead to the formation of a coalescing BH–NS (neutron star) system. The estimated Advanced LIGO/Virgo detection rate of such systems is up to ∼0.2 yr−1. This empirical estimate is a lower limit as it is obtained with only one particular evolutionary scenario, the MWC 656 binary. This is only a third such estimate available (after Cyg X-1 and Cyg X-3), and it lends additional support to the existence of so far undetected BH–NS binaries.

Constraining the GRB-Magnetar Model by Means of the Galactic Pulsar Population

A large fraction of Gamma-ray bursts (GRBs) displays an X-ray plateau phase within <105 s from the prompt emission, proposed to be powered by the spin-down energy of a rapidly spinning newly born magnetar. In this work we use the properties of the Galactic neutron star population to constrain the GRB-magnetar scenario. We re-analyze the X-ray plateaus of all Swift GRBs with known redshift, between 2005 January and 2014 August. From the derived initial magnetic field distribution for the possible magnetars left behind by the GRBs, we study the evolution and properties of a simulated GRB-magnetar population using numerical simulations of magnetic field evolution, coupled with Monte Carlo simulations of Pulsar Population Synthesis in our Galaxy. We find that if the GRB X-ray plateaus are powered by the rotational energy of a newly formed magnetar, the current observational properties of the Galactic magnetar population are not compatible with being formed within the GRB scenario (regardless of the GRB type or rate at z = 0). Direct consequences would be that we should allow the existence of magnetars and "super-magnetars" having different progenitors, and that Type Ib/c SNe related to Long GRBs form systematically neutron stars with higher initial magnetic fields. We put an upper limit of ≤16 "super-magnetars" formed by a GRB in our Galaxy in the past Myr (at 99% c.l.). This limit is somewhat smaller than what is roughly expected from Long GRB rates, although the very large uncertainties do not allow us to draw strong conclusion in this respect.

The differential evolution of hot- and cold-dust in galaxies in the last 7.5 GYR

Tese de mestrado em Física, apresentada à Universidade de Lisboa, através da Faculdade de Ciências, 2016

Examining the growth and evolution of magnetic fields in clusters of galaxies: a numerical perspective

Magnetic fields are ubiquitous in galaxy cluster atmospheres and have a variety of astrophysical and cosmological consequences. Magnetic fields can contribute to the pressure support of clusters, affect thermal conduction, and modify the evolution of bubbles driven by active galactic nuclei. However, we currently do not fully understand the origin and evolution of these fields throughout cosmic time. Furthermore, we do not have a general understanding of the relationship between magnetic field strength and topology and other cluster properties, such as mass and X-ray luminosity. We can now begin to answer some of these questions using large-scale cosmological magnetohydrodynamic (MHD) simulations of the formation of galaxy clusters including the seeding and growth of magnetic fields. Using large-scale cosmological simulations with the FLASH code combined with a simplified model of the acceleration of cosmic rays responsible for the generation of radio halos, we find that the galaxy cluster frequency distribution and expected number counts of radio halos from upcoming low-frequency sur- veys are strongly dependent on the strength of magnetic fields. Thus, a more complete understanding of the origin and evolution of magnetic fields is necessary to understand and constrain models of diffuse synchrotron emission from clusters. One favored model for generating magnetic fields is through the amplification of weak seed fields in active galactic nuclei (AGN) accretion disks and their subsequent injection into cluster atmospheres via AGN-driven jets and bubbles. However, current large-scale cosmological simulations cannot directly include the physical processes associated with the accretion and feedback processes of AGN or the seeding and merging of the associated SMBHs. Thus, we must include these effects as subgrid models. In order to carefully study the growth of magnetic fields in clusters via AGN-driven outflows, we present a systematic study of SMBH and AGN subgrid models. Using dark-matter only cosmological simulations, we find that many important quantities, such as the relationship between SMBH mass and galactic bulge velocity dispersion and the merger rate of black holes, are highly sensitive to the subgrid model assumptions of SMBHs. In addition, using MHD calculations of an isolated cluster, we find that magnetic field strengths, extent, topology, and relationship to other gas quantities such as temperature and density are also highly dependent on the chosen model of accretion and feedback. We use these systematic studies of SMBHs and AGN inform and constrain our choice of subgrid models, and we use those results to outline a fully cosmological MHD simulation to study the injection and growth of magnetic fields in clusters of galaxies. This simulation will be the first to study the birth and evolution of magnetic fields using a fully closed accretion-feedback cycle, with as few assumptions as possible and a clearer understanding of the effects of the various parameter choices.

Sub-kiloparsec alma imaging of compact star-forming galaxies at z ~ 2.5: revealing the formation of dense galactic cores in the progenitors of compact quiescent galaxies

We present spatially resolved Atacama Large Millimeter/submillimeter Array (ALMA) 870 μm dust continuum maps of six massive, compact, dusty star-forming galaxies at z ~ 2.5. These galaxies are selected for their small rest-frame optical sizes (r_e,F160W ~ 1.6 kpc) and high stellar mass densities that suggest that they are direct progenitors of compact quiescent galaxies at z ~ 2. The deep observations yield high far-infrared (FIR) luminosities of L_IR = 10^12.3-12.8 L_⨀ and star formation rates (SFRs) of SFR = 200–700 M_⊙ yr^−1, consistent with those of typical star-forming "main sequence" galaxies. The high spatial resolution (FWHM ~ 0 12–0 18) ALMA and Hubble Space Telescope photometry are combined to construct deconvolved, mean radial profiles of their stellar mass and (UV+IR) SFR. We find that the dusty, nuclear IR–SFR overwhelmingly dominates the bolometric SFR up to r ~ 5 kpc, by a factor of over 100× from the unobscured UV–SFR. Furthermore, the effective radius of the mean SFR profile (r_e,SFR ~ 1 kpc) is ~30% smaller than that of the stellar mass profile. The implied structural evolution, if such nuclear starburst last for the estimated gas depletion time of Δt = ±100 Myr, is a 4×increase of the stellar mass density within the central 1 kpc and a 1.6× decrease of the half-mass–radius. This structural evolution fully supports dissipation-driven, formation scenarios in which strong nuclear starbursts transform larger, star-forming progenitors into compact quiescent galaxies.

New indicators of star cluster dynamical evolution

This thesis presents a study of globular clusters (GCs), based on analysis of Monte Carlo simulations of globular clusters (GCs) with the aim to define new empirical parameters measurable from observations and able to trace the different phases of their dynamical evolution history. During their long term dynamical evolution, due to mass segregation and and dynamical friction, massive stars transfer kinetic energy to lower-mass objects, causing them to sink toward the cluster center. This continuous transfer of kinetic energy from the core to the outskirts triggers the runaway contraction of the core, known as "core collapse" (CC), followed by episodes of expansion and contraction called gravothermal oscillations. Clearly, such an internal dynamical evolution corresponds to significant variations also of the structure of the system. Determining the dynamical age of a cluster can be challenging as it depends on various internal and external properties. The traditional classification of GCs as CC or post-CC systems relies on detecting a steep power-law cusp in the central density profile, which may not always be reliable due to post-CC oscillations or other processes. In this thesis, based on the normalized cumulative radial distribution (nCRD) within a fraction of the half-mass radius is analyzed, and three diagnostics (A5, P5, and S2.5) are defined. These diagnostics show sensitivity to dynamical evolution and can distinguish pre-CC clusters from post-CC clusters.The analysis performed using multiple simulations with different initial conditions, including varying binary fractions and the presence of dark remnants showed the time variations of the diagnostics follow distinct patterns depending on the binary fraction and the retention or ejection of black holes. This analysis is extended to a larger set of simulations matching the observed properties of Galactic GCs, and the parameters show a potential to distinguish the dynamical stages of the observed clusters as well.

Prognostic Value Of Dna And Mrna E6/e7 Of Human Papillomavirus In The Evolution Of Cervical Intraepithelial Neoplasia Grade 2.

This study aimed at evaluating whether human papillomavirus (HPV) groups and E6/E7 mRNA of HPV 16, 18, 31, 33, and 45 are prognostic of cervical intraepithelial neoplasia (CIN) 2 outcome in women with a cervical smear showing a low-grade squamous intraepithelial lesion (LSIL). This cohort study included women with biopsy-confirmed CIN 2 who were followed up for 12 months, with cervical smear and colposcopy performed every three months. Women with a negative or low-risk HPV status showed 100% CIN 2 regression. The CIN 2 regression rates at the 12-month follow-up were 69.4% for women with alpha-9 HPV versus 91.7% for other HPV species or HPV-negative status (P < 0.05). For women with HPV 16, the CIN 2 regression rate at the 12-month follow-up was 61.4% versus 89.5% for other HPV types or HPV-negative status (P < 0.05). The CIN 2 regression rate was 68.3% for women who tested positive for HPV E6/E7 mRNA versus 82.0% for the negative results, but this difference was not statistically significant. The expectant management for women with biopsy-confirmed CIN 2 and previous cytological tests showing LSIL exhibited a very high rate of spontaneous regression. HPV 16 is associated with a higher CIN 2 progression rate than other HPV infections. HPV E6/E7 mRNA is not a prognostic marker of the CIN 2 clinical outcome, although this analysis cannot be considered conclusive. Given the small sample size, this study could be considered a pilot for future larger studies on the role of predictive markers of CIN 2 evolution.

Effects Of Therapeutic Approach On The Neonatal Evolution Of Very Low Birth Weight Infants With Patent Ductus Arteriosus.

To analyze the effects of treatment approach on the outcomes of newborns (birth weight [BW] < 1,000 g) with patent ductus arteriosus (PDA), from the Brazilian Neonatal Research Network (BNRN) on: death, bronchopulmonary dysplasia (BPD), severe intraventricular hemorrhage (IVH III/IV), retinopathy of prematurity requiring surgical (ROPsur), necrotizing enterocolitis requiring surgery (NECsur), and death/BPD. This was a multicentric, cohort study, retrospective data collection, including newborns (BW < 1000 g) with gestational age (GA) < 33 weeks and echocardiographic diagnosis of PDA, from 16 neonatal units of the BNRN from January 1, 2010 to Dec 31, 2011. Newborns who died or were transferred until the third day of life, and those with presence of congenital malformation or infection were excluded. Groups: G1 - conservative approach (without treatment), G2 - pharmacologic (indomethacin or ibuprofen), G3 - surgical ligation (independent of previous treatment). Factors analyzed: antenatal corticosteroid, cesarean section, BW, GA, 5 min. Apgar score < 4, male gender, Score for Neonatal Acute Physiology Perinatal Extension (SNAPPE II), respiratory distress syndrome (RDS), late sepsis (LS), mechanical ventilation (MV), surfactant (< 2 h of life), and time of MV. death, O2 dependence at 36 weeks (BPD36wks), IVH III/IV, ROPsur, NECsur, and death/BPD36wks. Student's t-test, chi-squared test, or Fisher's exact test; Odds ratio (95% CI); logistic binary regression and backward stepwise multiple regression. Software: MedCalc (Medical Calculator) software, version 12.1.4.0. p-values < 0.05 were considered statistically significant. 1,097 newborns were selected and 494 newborns were included: G1 - 187 (37.8%), G2 - 205 (41.5%), and G3 - 102 (20.6%). The highest mortality was observed in G1 (51.3%) and the lowest in G3 (14.7%). The highest frequencies of BPD36wks (70.6%) and ROPsur were observed in G3 (23.5%). The lowest occurrence of death/BPD36wks occurred in G2 (58.0%). Pharmacological (OR 0.29; 95% CI: 0.14-0.62) and conservative (OR 0.34; 95% CI: 0.14-0.79) treatments were protective for the outcome death/BPD36wks. The conservative approach of PDA was associated to high mortality, the surgical approach to the occurrence of BPD36wks and ROPsur, and the pharmacological treatment was protective for the outcome death/BPD36wks.

Evolution Of The Deaths Registry System In Brazil: Associations With Changes In The Mortality Profile, Under-registration Of Death Counts, And Ill-defined Causes Of Death.

This paper examines the spatial pattern of ill-defined causes of death across Brazilian regions, and its relationship with the evolution of completeness of the deaths registry and changes in the mortality age profile. We make use of the Brazilian Health Informatics Department mortality database and population censuses from 1980 to 2010. We applied demographic methods to evaluate the quality of mortality data for 137 small areas and correct for under-registration of death counts when necessary. The second part of the analysis uses linear regression models to investigate the relationship between, on the one hand, changes in death counts coverage and age profile of mortality, and on the other, changes in the reporting of ill-defined causes of death. The completeness of death counts coverage increases from about 80% in 1980-1991 to over 95% in 2000-2010 at the same time the percentage of ill-defined causes of deaths reduced about 53% in the country. The analysis suggests that the government's efforts to improve data quality are proving successful, and they will allow for a better understanding of the dynamics of health and the mortality transition.