Effects of high-energy astrophysical events on habitable planets

With the possible exception of meteor impacts, high-energy astrophysical events such as supernovae, gamma-ray bursts (GRB) and flares are usually not taken into account for biological and evolutionary studies due to their low rates of occurrence. We show that a class of these events may occur at distances and time scales in which their biological effects are non-negligible, maybe more frequent than the impacts of large asteroids. We review the effects of four transient astrophysical sources of ionizing radiation on biospheres - stellar flares, giant flares from soft gamma repeaters (SGR), supernovae and GRB. The main damaging features of them are briefly discussed and illustrated. We point out some open problems and ongoing work. Received 28 February 2012, accepted 6 July 2012, first published online 10 August 2012

Magnetars

Recent x-ray observations have shown that a substantial fraction of newly born neutron stars have magnetic fields of several 1014 G. They reveal themselves as soft gamma repeaters and anomalous x-ray pulsars and may account for the missing radio pulsars in young supernova remnants.

Role of Landau quantization on the neutron-drip transition in magnetar crusts

The role of a strong magnetic field on the neutron-drip transition in the crust of a magnetar is studied. The composition of the crust and the neutron-drip threshold are determined numerically for different magnetic field strengths using the experimental atomic mass measurements from the 2012 Atomic Mass Evaluation complemented with theoretical masses calculated from the Brussels-Montreal Hartree-Fock-Bogoliubov nuclear mass model HFB-24. The equilibrium nucleus at the neutron-drip point is found to be independent of the magnetic field strength. As demonstrated analytically, the neutron-drip density and pressure increase almost linearly with the magnetic field strength in the strongly quantizing regime for which electrons lie in the lowest Landau level. For weaker magnetic fields, the neutron-drip density exhibits typical quantum oscillations. In this case, the neutron-drip density can be either increased by about 14% or decreased by 25% depending on the magnetic field strength. These variations are shown to be almost universal, independently of the nuclear mass model employed. These results may have important implications for the physical interpretation of timing irregularities and quasiperiodic oscillations detected in soft gamma-ray repeaters and anomalous x-ray pulsars, as well as for the cooling of strongly magnetized neutron stars.

Is SGR 0418+5729 indeed a waning magnetar?

SGR 0418+5729 is a transient soft gamma-ray repeater which underwent a major outburst in 2009 June, during which the emission of short bursts was observed. Its properties appeared quite typical of other sources of the same class until long-term X-ray monitoring failed to detect any period derivative. The present upper limit on P implies that the surface dipole field is Bp lsim 7.5 × 1012 G, well below those measured in other soft gamma-ray repeaters (SGRs) and in the Anomalous X-ray Pulsars (AXPs), a group of similar sources. Both SGRs and AXPs are currently believed to be powered by ultra-magnetized neutron stars (magnetars, Bp ≈ 1014-1015 G). SGR 0418+5729 hardly seems to fit in such a picture. We show that the magneto-rotational properties of SGR 0418+5729 can be reproduced if this is an aged magnetar, ≈1 Myr old, which experienced substantial field decay. The large initial toroidal component of the internal field required to match the observed properties of SGR 0418+5729 ensures that crustal fractures, and hence bursting activity, can still occur at the present time. The thermal spectrum observed during the outburst decay is compatible with the predictions of a resonant Compton scattering model (as in other SGRs/AXPs) if the field is low and the magnetospheric twist is moderate.

Quiescent state and outburst evolution of SGR 0501+4516

We report on the quiescent state of the soft gamma repeater SGR 0501+4516 observed by XMM–Newton on 2009 August 30. The source exhibits an absorbed flux ∼75 times lower than that measured at the peak of the 2008 outburst, and a rather soft spectrum, with the same value of the blackbody temperature observed with ROSAT back in 1992. This new observation is put into the context of all existing X-ray data since its discovery in 2008 August, allowing us to complete the study of the timing and spectral evolution of the source from outburst until its quiescent state. The set of deep XMM–Newton observations performed during the few years time-scale of its outburst allows us to monitor the spectral characteristics of this magnetar as a function of its rotational period, and their evolution along these years. After the first ∼10 d, the initially hot and bright surface spot progressively cooled down during the decay. We discuss the behaviour of this magnetar in the context of its simulated secular evolution, inferring a plausible dipolar field at birth of 3 × 1014 G, and a current (magnetothermal) age of ∼10 kyr.

A VLT/FLAMES survey for massive binaries in Westerlund 1. IV. Wd1-5 – binary product and a pre-supernova companion for the magnetar CXOU J1647-45?

Context. The first soft gamma-ray repeater was discovered over three decades ago, and was subsequently identified as a magnetar, a class of highly magnetised neutron star. It has been hypothesised that these stars power some of the brightest supernovae known, and that they may form the central engines of some long duration gamma-ray bursts. However there is currently no consenus on the formation channel(s) of these objects. Aims. The presence of a magnetar in the starburst cluster Westerlund 1 implies a progenitor with a mass ≥40 M⊙, which favours its formation in a binary that was disrupted at supernova. To test this hypothesis we conducted a search for the putative pre-SN companion. Methods. This was accomplished via a radial velocity survey to identify high-velocity runaways, with subsequent non-LTE model atmosphere analysis of the resultant candidate, Wd1-5. Results. Wd1-5 closely resembles the primaries in the short-period binaries, Wd1-13 and 44, suggesting a similar evolutionary history, although it currently appears single. It is overluminous for its spectroscopic mass and we find evidence of He- and N-enrichement, O-depletion, and critically C-enrichment, a combination of properties that is difficult to explain under single star evolutionary paradigms. We infer a pre-SN history for Wd1-5 which supposes an initial close binary comprising two stars of comparable (~ 41 M⊙ + 35 M⊙) masses. Efficient mass transfer from the initially more massive component leads to the mass-gainer evolving more rapidly, initiating luminous blue variable/common envelope evolution. Reverse, wind-driven mass transfer during its subsequent WC Wolf-Rayet phase leads to the carbon pollution of Wd1-5, before a type Ibc supernova disrupts the binary system. Under the assumption of a physical association between Wd1-5 and J1647-45, the secondary is identified as the magnetar progenitor; its common envelope evolutionary phase prevents spin-down of its core prior to SN and the seed magnetic field for the magnetar forms either in this phase or during the earlier episode of mass transfer in which it was spun-up. Conclusions. Our results suggest that binarity is a key ingredient in the formation of at least a subset of magnetars by preventing spin-down via core-coupling and potentially generating a seed magnetic field. The apparent formation of a magnetar in a Type Ibc supernova is consistent with recent suggestions that superluminous Type Ibc supernovae are powered by the rapid spin-down of these objects.

Instabilities of soft elastic microtubes filled with viscous fluids: Pearls, wrinkles, and sausage strings

A linear stability analysis is presented to study the self-organized instabilities of a highly compliant elastic cylindrical shell filled with a viscous liquid and submerged in another viscous medium. The prototype closely mimics many components of micro-or nanofluidic devices and biological processes such as the budding of a string of pearls inside cells and sausage-string formation of blood vessels. The cylindrical shell is considered to be a soft linear elastic solid with small storage modulus. When the destabilizing capillary force derived from the cross-sectional curvature overcomes the stabilizing elastic and in-plane capillary forces, the microtube can spontaneously self-organize into one of several possible configurations; namely, pearling, in which the viscous fluid in the core of the elastic shell breaks up into droplets; sausage strings, in which the outer interface of the mircrotube deforms more than the inner interface; and wrinkles, in which both interfaces of the thin-walled mircrotube deform in phase with small amplitudes. This study identifies the conditions for the existence of these modes and demonstrates that the ratios of the interfacial tensions at the interfaces, the viscosities, and the thickness of the microtube play crucial roles in the mode selection and the relative amplitudes of deformations at the two interfaces. The analysis also shows asymptotically that an elastic fiber submerged in a viscous liquid is unstable for Y = gamma/(G(e)R) > 6 and an elastic microchannel filled with a viscous liquid should rupture to form spherical cavities (pearling) for Y > 2, where gamma, G(e), and R are the surface tension, elastic shear modulus, and radius, respectively, of the fiber or microchannel.

Bearing capacity of piles in soft clay underlaid by cohesive frictional soil

By using the axisymmetric quasi-lower bound finite-element limit analysis, the bearing capacity factors N-c(p) and N-gamma q(p) have been computed for axially loaded piles, with the shaft embedded in a fully cohesive soil medium and the tip placed over cohesive frictional soil strata. The results were obtained for various combinations of L/D, phi(l), and c(l)/c(u); the subscripts l and u refer to lower and upper soil strata, respectively. The factors N-c(p) and N-gamma q(p) increase continuously with increases in L/D and phi(l); the rate of increase of N-c(p) and N-gamma q(p) with L/D, however, decreases with an increase in L/D. For c(l)/c(u) > 100, the factor N-c(p) hardly depends on L/D.

Gamma Knife Radiosurgery as a Therapeutic Strategy for Intracranial Sarcomatous Metastases

Purpose: To determine the indication and outcomes for Gamma Knife stereotactic radiosurgery (GKSRS) in the care of patients with intracranial sarcomatous metastases. Methods and Materials: Data from 21 patients who underwent radiosurgery for 60 sarcomatous intracranial metastases (54 parenchymal and 6 dural-based) were studied. Nine patients had radiosurgery for solitary tumors and 12 for multiple tumors. The primary pathology was metastatic leiomyosarcoma (4 patients), osteosarcoma (3 patients), soft-tissue sarcoma (5 patients), chondrosarcoma (2 patients), alveolar soft part sarcoma (2 patients), and rhabdomyosarcoma, Ewing's sarcoma, liposarcoma, neurofibrosarcoma, and synovial sarcoma (1 patient each). Twenty patients received multimodality management for their primary tumor, and 1 patient had no evidence of systemic disease. The mean tumor volume was 6.2 cm 3 (range, 0.07-40.9 cm 3), and a median margin dose of 16 Gy was administered. Three patients had progressive intracranial disease despite fractionated whole-brain radiotherapy before SRS. Results: A local tumor control rate of 88% was achieved (including patients receiving boost, up-front, and salvage SRS). New remote brain metastases developed in 7 patients (33%). The median survival after diagnosis of intracranial metastasis was 16 months, and the 1-year survival rate was 61%. Conclusions: Gamma Knife radiosurgery was a well-tolerated and initially effective therapy in the management of patients with sarcomatous intracranial metastases. However, many patients, including those who also received fractionated whole-brain radiotherapy, developed progressive new brain disease. © 2010 Elsevier Inc. All rights reserved.

Electron impact excitation of Mg VIII: Collision strengths, transition probabilities and theoretical EUV and soft X-ray line intensities for Mg VIII

Context: Mg VIII emission lines are observed in a range of astronomical objects such as the Sun, other cool stars and in the coronal line region of Seyfert galaxies. Under coronal conditions Mg VIII emits strongly in the extreme ultraviolet (EUV) and soft X-ray spectral regions which makes it an ideal ion for plasma diagnostics.

Aims. Two theoretical atomic models, consisting of 125 fine structure levels, are developed for the Mg VIII ion. The 125 levels arise from the 2s(2)2p, 2s(2)p2, 2p(3), 2s(2)3s, 2s(2)3p, 2s(2)3d, 2s2p3s, 2s2p3p, 2s2p3d, 2p(2)3s, 2p(2)3p and 2p(2)3d configurations. Electron impact excitation collision strengths and radiative transition probabilities are calculated for both Mg VIII models, compared with existing data, and the best model selected to generate a set of theoretical emission line intensities. The EUV lines, covering 312-790 angstrom, are compared with existing solar spectra (SERTS-89 and SUMER), while the soft X-ray transitions (69-97 angstrom) are examined for potential density diagnostic line ratios and also compared with the limited available solar and stellar observational data.

Methods. The R-matrix codes Breit-Pauli RMATRXI and RMATRXII are utilised, along with the PSTGF code, to calculate the collision strengths for two Mg VIII models. Collision strengths are averaged over a Maxwellian distribution to produce the corresponding effective collision strengths for use in astrophysical applications. Transition probabilities are also calculated using the CIV3 atomic structure code. The best data are then incorporated into the modelling code CLOUDY and line intensities generated for a range of electron temperatures and densities appropriate to solar and stellar coronal plasmas.

Results. The present effective collision strengths are compared with two previous calculations. Good levels of agreement are found with the most recent, but there are large differences with the other for forbidden transitions. The resulting line intensities compare favourably with the observed values from the SERTS-89 and SUMER spectra. Theoretical soft X-ray emission lines are presented and several density diagnostic line ratios examined, which are in reasonable agreement with the limited observational data available.

Nanostructure formation in poly(gamma-benzyl-L-glutamate)-poly(ethylene glycol)-poly(gamma-benzyl-L-glutamate) triblock copolymers in the solid state

The morphology in the solid state of a series of triblock copolymers comprising a poly(ethylene glycol) (PEG) midblock and symmetric poly(gamma-benzyl-L-glutamate) (PBLG) end blocks has been studied using X-ray scattering and microscopy techniques. Transmission electron microscopy (TEM) on samples selectively stained with uranyl acetate provided clear assignment of morphologies for as-cast and annealed samples. The thickness of both PEG and PBLG domains was in good agreement with calculations based on the conformations of the respective chains, allowing for the crystal or amorphous state of PEG and the a-helical or P-sheet structure of the PBLG. Atomic force microscopy provided complementary information on surface morphology for several samples that was in good agreement with the structure observed by TEM. A morphology diagram was constructed. Cylindrical structures were observed for ordered samples with low f(PBLG), whereas at higher f(PLBG) there was evidence for broken lamellar and "hockey puck" nanostructures. Regular lamellae were observed for intermediate compositions.

DECIPHERING THE QUARK-GLUON STRUCTURE OF THE PHOTON IN E-GAMMA COLLISIONS

We investigate the capability of an egamma collider to unravel the hadronic content of the photon. The experimental problem for probing the gluonic structure of the photon is that small-x triggers overwhelmingly select soft photons rather than soft gluons in hard photons. We show that the problem can be finessed in experiments where laser back-scattering is used to prepare a source of very hard photons. We illustrate their power for studying the parton distributions of the photon and, specifically, for separating the quark and gluon components in events where dijets, jet-gamma pairs, and heavy quark pairs are produced.

Modeling and simulations of some anisotropic soft-matter systems: from biaxial to chiral materials

We have modeled various soft-matter systems with molecular dynamics (MD) simulations. The first topic concerns liquid crystal (LC) biaxial nematic (Nb) phases, that can be possibly used in fast displays. We have investigated the phase organization of biaxial Gay-Berne (GB) mesogens, considering the effects of the orientation, strength and position of a molecular dipole. We have observed that for systems with a central dipole, nematic biaxial phases disappear when increasing dipole strength, while for systems characterized by an offset dipole, the Nb phase is stabilized at very low temperatures. In a second project, in view of their increasing importance as nanomaterials in LC phases, we are developing a DNA coarse-grained (CG) model, in which sugar and phosphate groups are represented with Lennard-Jones spheres, while bases with GB ellipsoids. We have obtained shape, position and orientation parameters for each bead, to best reproduce the atomistic structure of a B-DNA helix. Starting from atomistic simulations results, we have completed a first parametrization of the force field terms, accounting for bonded (bonds, angles and dihedrals) and non-bonded interactions (H-bond and stacking). We are currently validating the model, by investigating stability and melting temperature of various sequences. Finally, in a third project, we aim to explain the mechanism of enantiomeric discrimination due to the presence of a chiral helix of poly(gamma-benzyl L-glutamate) (PBLG), in solution of dimethylformamide (DMF), interacting with chiral or pro-chiral molecules (in our case heptyl butyrate, HEP), after tuning properly an atomistic force field (AMBER). We have observed that DMF and HEP molecules solvate uniformly the PBLG helix, but the pro-chiral solute is on average found closer to the helix with respect to the DMF. The solvent presents a faster isotropic diffusion, twice as HEP, also indicating a stronger interaction of the solute with the helix.

Factorization and non-local 1/m b corrections in the decay B X s gamma

In this thesis, a systematic analysis of the bar B to X_sgamma photon spectrum in the endpoint region is presented. The endpoint region refers to a kinematic configuration of the final state, in which the photon has a large energy m_b-2E_gamma = O(Lambda_QCD), while the jet has a large energy but small invariant mass. Using methods of soft-collinear effective theory and heavy-quark effective theory, it is shown that the spectrum can be factorized into hard, jet, and soft functions, each encoding the dynamics at a certain scale. The relevant scales in the endpoint region are the heavy-quark mass m_b, the hadronic energy scale Lambda_QCD and an intermediate scale sqrt{Lambda_QCD m_b} associated with the invariant mass of the jet. It is found that the factorization formula contains two different types of contributions, distinguishable by the space-time structure of the underlying diagrams. On the one hand, there are the direct photon contributions which correspond to diagrams with the photon emitted directly from the weak vertex. The resolved photon contributions on the other hand arise at O(1/m_b) whenever the photon couples to light partons. In this work, these contributions will be explicitly defined in terms of convolutions of jet functions with subleading shape functions. While the direct photon contributions can be expressed in terms of a local operator product expansion, when the photon spectrum is integrated over a range larger than the endpoint region, the resolved photon contributions always remain non-local. Thus, they are responsible for a non-perturbative uncertainty on the partonic predictions. In this thesis, the effect of these uncertainties is estimated in two different phenomenological contexts. First, the hadronic uncertainties in the bar B to X_sgamma branching fraction, defined with a cut E_gamma > 1.6 GeV are discussed. It is found, that the resolved photon contributions give rise to an irreducible theory uncertainty of approximately 5 %. As a second application of the formalism, the influence of the long-distance effects on the direct CP asymmetry will be considered. It will be shown that these effects are dominant in the Standard Model and that a range of -0.6 < A_CP^SM < 2.8 % is possible for the asymmetry, if resolved photon contributions are taken into account.

NLL QCD contribution of the electromagnetic dipole operator to B -> Xs gamma gamma with a massive strange quark

We calculate the O(αs) corrections to the double differential decay width dΓ77/(ds1ds2) for the process B¯→Xsγγ, originating from diagrams involving the electromagnetic dipole operator O7. The kinematical variables s1 and s2 are defined as si=(pb−qi)2/m2b, where pb, q1, q2 are the momenta of the b quark and two photons. We introduce a nonzero mass ms for the strange quark to regulate configurations where the gluon or one of the photons become collinear with the strange quark and retain terms which are logarithmic in ms, while discarding terms which go to zero in the limit ms→0. When combining virtual and bremsstrahlung corrections, the infrared and collinear singularities induced by soft and/or collinear gluons drop out. By our cuts the photons do not become soft, but one of them can become collinear with the strange quark. This implies that in the final result a single logarithm of ms survives. In principle, the configurations with collinear photon emission could be treated using fragmentation functions. In a related work we find that similar results can be obtained when simply interpreting ms appearing in the final result as a constituent mass. We do so in the present paper and vary ms between 400 and 600 MeV in the numerics. This work extends a previous paper by us, where only the leading power terms with respect to the (normalized) hadronic mass s3=(pb−q1−q2)2/m2b were taken into account in the underlying triple differential decay width dΓ77/(ds1ds2ds3).