Multi-dimensional simulations of mixing in classical novae

Classical nova explosions are stellar explosions that take place in close binary systems with an energy release only exceeded by gamma-ray bursts and supernova explosions. Matter from the white dwarf flows through the inner lagrangian point and spirals in towards the white dwarf for about 10^4-10^5 years, forming an accretion disk around it. Ultimately, part of this hydrogen-rich matter piles-up on top of the compact object and becomes partially degenerate due to the high densities attained. Consequently, temperature is allowed to rise, but the envelope does not experience any expansion. Actually, this is the key mechanism that controls the subsequent phases and powers a thermonuclear runaway, which is followed by an ejection of part of the accreted envelope. The ejecta are enriched with the products from the nuclear processes, presenting a final metallicity much above solar. This model, introduced in the early 70s, is a solid theory that can account for the gross scenario of nova explosions. Nevertheless, the theory relies on the fact that a mixing episode with matter from the white dwarf core has to take place at the core-envelope interface to successfully account for the high metallicities inferred from observations. During the past 40 years, theoreticians have performed many one-dimensional simulations, which can reproduce the abundances in the ejecta and other important observational properties. However, these calculations performed in spherical symmetry cannot study the mixing process, since they exclude a suite of very important multi-dimensional effects, such as convection. Therefore, multi-dimensional calculations are required to shed light into the mixing episode. In this thesis we have performed two- and three- dimensional simulations of CO novae to study the mixing mechanisms operating at the core-envelope interface, how convection sets in and how the deflagration spreads over the domain, by means of the Eulerian, parallelized, hydrodynamical FLASH code. The two-dimensional results show how convection sets in at the innermost envelope layers, after the appearance of temperature fluctuations that arise from the interface. Convection, in turn, powers the formation of kelvin-Helmholtz instabilities, which efficiently dredge-up 12C from the core and carry it into the envelope, reproducing correctly the high metallicity found in the ejecta. This result solves the controversy generated by the two existing two-dimensional calculations up-to-date. We have also realized a sensitivity study to analyze the impact of some initial parameters, such as the temperature perturbation, resolution of the simulations and the size of the computational domain. The results point out that these parameters have a negligible impact on the degree of mixing and, therefore, the calculations are not affected by numerical artifacts. Although two-dimensional calculations can quantitatively reproduce the mixing episode, they cannot describe correctly the convective pattern due to conservation of vorticity, which translates into recombination of the convective cells. Therefore, we have extended the work to three dimensions and performed the first three-dimesional model of mixing in classical novae up-to-date. These calculations can successfully reproduce the intermittency present in turbulent convection, with an energy cascade into smaller scales which clearly fulfills the Kolmogorov theory, while the thermonuclear runaway continues propagating with almost spherical symmetry. Mixing proceeds through the filamentary structure powered by robust kelvin-Helmholtz instabilitites that arise from the interface, resulting in a CNO enhancement which agrees with observations. This convective profile also generates density contrasts that could be the origin of the inhomogeneous distribution of chemical species.Les explosions de noves tenen lloc en un sistema estel.lar binari, on un dels estels ha arribat a la fi de la seva vida convertit en una nana blanca. En sistemes binaris molt propers, l'estel acompanyant cedeix part del seu gas (material ric en hidrogen), el qual s'arremolina al voltant de la nana blanca durant prop de 10^4 - 10^5 anys. Una fracció d'aquest material acaba apilant-se a la superfície de l'objecte compacte i esdevé parcialment degenerat com a conseqüència de l'elevada densitat. Aquest fet és clau en el procés, ja que permet que la temperatura augmenti sense que es produeixi una expansió de l'embolcall, desencadenant un allau termonuclear i finalment, l'ejecció de matèria. El material ejectat està enriquit amb els isòtops processats en les reaccions nuclears, presentant una metal.licitat molt superior a la solar. Aquest model, presentat a principis dels anys 70, és una teoria sòlida que explica raonablement l'explosió de noves. No obstant, la teoria rau en el fet que s'ha de produir un procés de barreja entre el material de la nana blanca i el material de les capes més internes de l'embolcall per poder explicar l'alta metal.licitat que s'observa en el material ejectat. Durant els últims 40 anys, s'han fet molts estudis en una dimensió que aconsegueixen reproduir correctament les abundàncies del material ejectat i altres importants propietats observacionals, però que no poden explicar com es produeix el procés de barreja, ja que aquests càlculs amb simetria esfèrica exlouen tota una sèrie d'importants fenòmens multidimensionals. Per tant, per estudiar aquests aspectes de la teoria es requereixen estudis multidimensionals. En aquesta tesi hem realitzat simulacions en dues i tres dimensions de noves de CO per estudiar els mecanismes de barreja que es produeixen a la interfície del nucli de la nana blanca i l'embolcall, com s'estableix la convecció i com es propaga el front deflagratiu, mitjançant el codi hidrodinàmic FLASH, que és Eulerià i està paral.lelitzat. Els resultats en dues dimensions mostren com es genera convecció a les capes més internes de l'embolcall, després de la formació de fluctuacions de temperatura a la interfície. La convecció, al seu torn, origina inestabilitats Kelvin-Helmholtz que transporten eficientment 12C del nucli cap a l'embolcall, aconseguint reproduir correctament el grau de metal.licitat observat. Aquest resultat resol la controvèrsia generada pels dos estudis en dues dimensions realitzats fins ara. També hem realitzat un estudi per analitzar l'impacte dels paràmetres inicials tals com la perturbació inicial, la resolució de les simulacions o les dimensions del domini computacional. Els resultats indiquen que cap d'aquests paràmetres influeix en el grau de barreja final i, per tant, que els càlculs no estan condicionats per aspectes numèrics. Finalment, hem presentat el primer model tridimensional de barreja de noves fet fins ara. Aquest càlcul és necessari, ja que les simulacions bidimensionals, tot i que quantitativament reprodueixen la barreja esperada, no poden representar el patró convectiu correctament, degut a la conservació de la vorticitat, fent que les cel.les convectives esdevinguin cada cop més grans. El nostre càlcul aconsegueix reproduir el comportament intermitent de la turbulència, amb una cascada d'energia que flueix cap a escales cada cop més petites, tal i com prediu la teoria de Kolmogorov, alhora que el front convectiu avança pràcticament amb simetria esfèrica. La barreja procedeix a través de l'estructura filamentosa originada per l'aparició de potents inestabilitats Kelvin-Helmholtz a la interfície, obtenint-se una metal.licitat final a l'embolcall que concorda amb els valors observacionals. Aquest patró convectiu també genera contrastos de densitat que podrien ser l'origen de la distribució inhomogènia que presenten les espècies químiques.Postprint (published version

Tests hidrodinàmics amb el codi multidimensional FLASH

El FLASH és un codi estable i flexible, desenvolupat a la Universitat de Chicago, que permet estudiar diferents fenòmens hidrodinàmics. L'objectiu futur és desenvolupar models multidimensionals centrats en el camp de les noves per a caracteritzar algunes fases complexes del fenomen, però degut a la complexitat del codi (més de 500000 línies de programació en C++ i Fortran90), la realització d'una simulació estel.lar necessita tenir un coneixement bàsic del FLASH. Per tant, l'objectiu d'aquest treball és la realització d'alguns tests numèrics per aprendre'n el funcionament. Però a la vegada, l'objectiu és doble, ja que per utilitzar una eina de treball com el FLASH, necessitem testejar el propi codi i comprovar que és capaç de resoldre correctament diferents problemes hidrodinàmics que ja s'han simulat prèviament. Així, s'han treballat els diferents tests i hem observat que el codi FLASH és eficaç per simular-los correctament i, per tant, podem esperar que serà un bon codi per a les futures simulacions estel.lars. D'altra banda, aquest primer contacte amb el codi ens ha permès familiaritzar-nos amb l'estructura, adquirir destresa en el seu ús i començar a entendre'l, tot deixant per a més endavant l'estudi detallat dels diferents models numèrics implementats en el codi i eines auxiliars de representació gràfica

Recent advances in the modelling of classical novae and type I X-ray bursts

Classical nova outbursts and type I X-ray bursts are thermonuclear stellar explosions driven by charged-particle reactions. Extensive numerical simulations of nova explosions have shown that the accreted envelopes attain peak temperatures between 0.1 and 0.4 GK, for about several hundred seconds, and therefore, their ejecta is expected to show signatures of significant nuclear activity. Indeed, it has been claimed that novae play some role in the enrichment of the interstellar medium through a number of intermediate-mass elements. This includes 17O, 15N, and 13C, systematically overproduced in huge amounts with respect to solar abundances, with a lower contribution to a number of species with A<40, such as 7Li, 19F, or 26Al. In this review, we present new 1-D hydrodynamic models of classical nova outbursts, from the onset of accretion up to the explosion and ejection phases. Special emphasis is put on their gross observational properties (including constraints from meteoritic presolar grains and potential gamma-ray signatures) and on their associated nucleosynthesis. Multidimensional models of mixing at the core-envelope interface during outbursts will also be presented. The impact of nuclear uncertainties on the final yields will be also outlined. Detailed analysis of the relevant reactions along the main nuclear path for type I X-ray bursts has only been scarcely addressed, mainly in the context of parameterized one-zone models. Here, we present a detailed study of the nucleosynthesis and nuclear processes powering type I X-ray bursts. The reported bursts have been computed by means of a spherically symmetric (1D), Lagrangian, hydrodynamic code, linked to a nuclear reaction network that contains 325 isotopes (from 1H to 107Te), and 1392 nuclear processes. These evolutionary sequences, followed from the onset of accretion up to the explosion and expansion stages, have been performed for two different metallicities to explore the dependence between the extension of the main nuclear flow and the initial metal content. We carefully analyze the physical parameters that determine the light curve (including recurrence times, ratios between persistent and burst luminosities, or the extent of the envelope expansion). Results are in qualitative agreement with the observed properties of some well-studied bursting sources.Postprint (published version

Autoantibodies against type I IFNs in patients with critical influenza pneumonia

In an international cohort of 279 patients with hypoxemic influenza pneumonia, we identified 13 patients (4.6%) with autoantibodies neutralizing IFN-alpha and/or -omega, which were previously reported to underlie 15% cases of life-threatening COVID-19 pneumonia and one third of severe adverse reactions to live-attenuated yellow fever vaccine. Autoantibodies neutralizing type I interferons (IFNs) can underlie critical COVID-19 pneumonia and yellow fever vaccine disease. We report here on 13 patients harboring autoantibodies neutralizing IFN-alpha 2 alone (five patients) or with IFN-omega (eight patients) from a cohort of 279 patients (4.7%) aged 6-73 yr with critical influenza pneumonia. Nine and four patients had antibodies neutralizing high and low concentrations, respectively, of IFN-alpha 2, and six and two patients had antibodies neutralizing high and low concentrations, respectively, of IFN-omega. The patients' autoantibodies increased influenza A virus replication in both A549 cells and reconstituted human airway epithelia. The prevalence of these antibodies was significantly higher than that in the general population for patients 70 yr of age (3.1 vs. 4.4%, P = 0.68). The risk of critical influenza was highest in patients with antibodies neutralizing high concentrations of both IFN-alpha 2 and IFN-omega (OR = 11.7, P = 1.3 x 10(-5)), especially those <70 yr old (OR = 139.9, P = 3.1 x 10(-10)). We also identified 10 patients in additional influenza patient cohorts. Autoantibodies neutralizing type I IFNs account for similar to 5% of cases of life-threatening influenza pneumonia in patients <70 yr old

Tests hidrodinàmics amb el codi multidimensional FLASH

El FLASH és un codi estable i flexible, desenvolupat a la Universitat de Chicago, que permet estudiar diferents fenòmens hidrodinàmics. L'objectiu futur és desenvolupar models multidimensionals centrats en el camp de les noves per a caracteritzar algunes fases complexes del fenomen, però degut a la complexitat del codi (més de 500000 línies de programació en C++ i Fortran90), la realització d'una simulació estel.lar necessita tenir un coneixement bàsic del FLASH. Per tant, l'objectiu d'aquest treball és la realització d'alguns tests numèrics per aprendre'n el funcionament. Però a la vegada, l'objectiu és doble, ja que per utilitzar una eina de treball com el FLASH, necessitem testejar el propi codi i comprovar que és capaç de resoldre correctament diferents problemes hidrodinàmics que ja s'han simulat prèviament. Així, s'han treballat els diferents tests i hem observat que el codi FLASH és eficaç per simular-los correctament i, per tant, podem esperar que serà un bon codi per a les futures simulacions estel.lars. D'altra banda, aquest primer contacte amb el codi ens ha permès familiaritzar-nos amb l'estructura, adquirir destresa en el seu ús i començar a entendre'l, tot deixant per a més endavant l'estudi detallat dels diferents models numèrics implementats en el codi i eines auxiliars de representació gràfica

Recent advances in the modelling of classical novae and type I X-ray bursts

Classical nova outbursts and type I X-ray bursts are thermonuclear stellar explosions driven by charged-particle reactions. Extensive numerical simulations of nova explosions have shown that the accreted envelopes attain peak temperatures between 0.1 and 0.4 GK, for about several hundred seconds, and therefore, their ejecta is expected to show signatures of significant nuclear activity. Indeed, it has been claimed that novae play some role in the enrichment of the interstellar medium through a number of intermediate-mass elements. This includes 17O, 15N, and 13C, systematically overproduced in huge amounts with respect to solar abundances, with a lower contribution to a number of species with A<40, such as 7Li, 19F, or 26Al. In this review, we present new 1-D hydrodynamic models of classical nova outbursts, from the onset of accretion up to the explosion and ejection phases. Special emphasis is put on their gross observational properties (including constraints from meteoritic presolar grains and potential gamma-ray signatures) and on their associated nucleosynthesis. Multidimensional models of mixing at the core-envelope interface during outbursts will also be presented. The impact of nuclear uncertainties on the final yields will be also outlined. Detailed analysis of the relevant reactions along the main nuclear path for type I X-ray bursts has only been scarcely addressed, mainly in the context of parameterized one-zone models. Here, we present a detailed study of the nucleosynthesis and nuclear processes powering type I X-ray bursts. The reported bursts have been computed by means of a spherically symmetric (1D), Lagrangian, hydrodynamic code, linked to a nuclear reaction network that contains 325 isotopes (from 1H to 107Te), and 1392 nuclear processes. These evolutionary sequences, followed from the onset of accretion up to the explosion and expansion stages, have been performed for two different metallicities to explore the dependence between the extension of the main nuclear flow and the initial metal content. We carefully analyze the physical parameters that determine the light curve (including recurrence times, ratios between persistent and burst luminosities, or the extent of the envelope expansion). Results are in qualitative agreement with the observed properties of some well-studied bursting sources

Towards a theory of mediated communication

SIGLEAvailable from British Library Document Supply Centre- DSC:4335.26205(HPL--97-78) / BLDSC - British Library Document Supply CentreGBUnited Kingdo

Classical and recurrent nova models

Remarkable progress in the understanding of nova outbursts has been achieved through combined efforts in photometry, spectroscopy and numerical simulations. According to the thermonuclear runaway model, novae are powered by thermonuclear explosions in the hydrogen-rich envelopes transferred from a low-mass stellar companion onto a close white dwarf star. Extensive numerical simulations in 1-D have shown that the accreted envelopes attain peak temperatures ranging between 108 and 4 × 108 K, for about several hundred seconds, hence allowing extensive nuclear processing which eventually shows up in the form of nucleosynthetic fingerprints in the ejecta. Indeed, it has been claimed that novae can play a certain role in the enrichment of the interstellar medium through a number of intermediate-mass elements. This includes 17O, 15N, and 13C, systematically overproduced with respect to solar abundances, plus a lower contribution in a number of other species (A < 40), such as 7Li, 19F, or 26Al. At the turn of the XXI Century, classical novae have entered the era of multidimensional models, which provide a new insight into the physical mechanisms that drive mixing at the core-envelope interface. In this review, we will present hydrodynamic models of classical novae, from the onset of accretion up to the explosion and ejection stages, both for classical and recurrent novae, with special emphasis on their gross observational properties and their associated nucleosynthesis. The impact of nuclear uncertainties on the final yields will be discussed. Recent results from 2-D models of mixing during classical nova outbursts will also be presented

Assessment of plasma chitotriosidase activity, CCL18/PARC concentration and NP-C suspicion index in the diagnosis of Niemann-Pick disease type C : A prospective observational study

Niemann-Pick disease type C (NP-C) is a rare, autosomal recessive neurodegenerative disease caused by mutations in either the NPC1 or NPC2 genes. The diagnosis of NP-C remains challenging due to the non-specific, heterogeneous nature of signs/symptoms. This study assessed the utility of plasma chitotriosidase (ChT) and Chemokine (C-C motif) ligand 18 (CCL18)/pulmonary and activation-regulated chemokine (PARC) in conjunction with the NP-C suspicion index (NP-C SI) for guiding confirmatory laboratory testing in patients with suspected NP-C. In a prospective observational cohort study, incorporating a retrospective determination of NP-C SI scores, two different diagnostic approaches were applied in two separate groups of unrelated patients from 51 Spanish medical centers (n = 118 in both groups). From Jan 2010 to Apr 2012 (Period 1), patients with ≥2 clinical signs/symptoms of NP-C were considered 'suspected NP-C' cases, and NPC1/NPC2 sequencing, plasma chitotriosidase (ChT), CCL18/PARC and sphingomyelinase levels were assessed. Based on findings in Period 1, plasma ChT and CCL18/PARC, and NP-C SI prediction scores were determined in a second group of patients between May 2012 and Apr 2014 (Period 2), and NPC1 and NPC2 were sequenced only in those with elevated ChT and/or elevated CCL18/PARC and/or NP-C SI ≥70. Filipin staining and 7-ketocholesterol (7-KC) measurements were performed in all patients with NP-C gene mutations, where possible. In total across Periods 1 and 2, 10/236 (4%) patients had a confirmed diagnosis o NP-C based on gene sequencing (5/118 [4.2%] in each Period): all of these patients had two causal NPC1 mutations. Single mutant NPC1 alleles were detected in 8/236 (3%) patients, overall. Positive filipin staining results comprised three classical and five variant biochemical phenotypes. No NPC2 mutations were detected. All patients with NPC1 mutations had high ChT activity, high CCL18/PARC concentrations and/or NP-C SI scores ≥70. Plasma 7-KC was higher than control cut-off values in all patients with two NPC1 mutations, and in the majority of patients with single mutations. Family studies identified three further NP-C patients. This approach may be very useful for laboratories that do not have mass spectrometry facilities and therefore, they cannot use other NP-C biomarkers for diagnosis

Inborn errors of OAS-RNase L in SARS-CoV-2-related multisystem inflammatory syndrome in children

Multisystem inflammatory syndrome in children (MIS-C) is a rare and severe condition that follows benign COVID-19. We report autosomal recessive deficiencies of OAS1, OAS2, or RNASEL in five unrelated children with MIS-C. The cytosolic double-stranded RNA (dsRNA)-sensing OAS1 and OAS2 generate 2'-5'-linked oligoadenylates (2-5A) that activate the single-stranded RNA-degrading ribonuclease L (RNase L). Monocytic cell lines and primary myeloid cells with OAS1, OAS2, or RNase L deficiencies produce excessive amounts of inflammatory cytokines upon dsRNA or severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) stimulation. Exogenous 2-5A suppresses cytokine production in OAS1-deficient but not RNase L-deficient cells. Cytokine production in RNase L-deficient cells is impaired by MDA5 or RIG-I deficiency and abolished by mitochondrial antiviral-signaling protein (MAVS) deficiency. Recessive OAS-RNase L deficiencies in these patients unleash the production of SARS-CoV-2-triggered, MAVS-mediated inflammatory cytokines by mononuclear phagocytes, thereby underlying MIS-C