Self-Organising Maps and Correlation Analysis as a Tool to Explore Patterns in Excitation-Emission Matrix Data Sets and to Discriminate Dissolved Organic Matter Fluorescence Components

Dissolved organic matter (DOM) is a complex mixture of organic compounds, ubiquitous in marine and freshwater systems. Fluorescence spectroscopy, by means of Excitation-Emission Matrices (EEM), has become an indispensable tool to study DOM sources, transport and fate in aquatic ecosystems. However the statistical treatment of large and heterogeneous EEM data sets still represents an important challenge for biogeochemists. Recently, Self-Organising Maps (SOM) has been proposed as a tool to explore patterns in large EEM data sets. SOM is a pattern recognition method which clusterizes and reduces the dimensionality of input EEMs without relying on any assumption about the data structure. In this paper, we show how SOM, coupled with a correlation analysis of the component planes, can be used both to explore patterns among samples, as well as to identify individual fluorescence components. We analysed a large and heterogeneous EEM data set, including samples from a river catchment collected under a range of hydrological conditions, along a 60-km downstream gradient, and under the influence of different degrees of anthropogenic impact. According to our results, chemical industry effluents appeared to have unique and distinctive spectral characteristics. On the other hand, river samples collected under flash flood conditions showed homogeneous EEM shapes. The correlation analysis of the component planes suggested the presence of four fluorescence components, consistent with DOM components previously described in the literature. A remarkable strength of this methodology was that outlier samples appeared naturally integrated in the analysis. We conclude that SOM coupled with a correlation analysis procedure is a promising tool for studying large and heterogeneous EEM data sets.

A Be-type star with a black-hole companion

Stellar-mass black holes have all been discovered through X-ray emission, which arises from the accretion of gas from their binary companions (this gas is either stripped from low-mass stars or supplied as winds from massive ones). Binary evolution models also predict the existence of black holes accreting from the equatorial envelope of rapidly spinning Be-type stars (stars of the Be type are hot blue irregular variables showing characteristic spectral emission lines of hydrogen). Of the ~80 Be X-ray binaries known in the Galaxy, however, only pulsating neutron stars have been found as companions. A black hole was formally allowed as a solution for the companion to the Be star MWC 656 (also known as HD 215227), although that was based on a single radial velocity curve of the Be star, a mistaken spectral classification and rough estimates of the inclination angle. Here we report observations of an accretion disk line mirroring the orbit of the Be star. This, together with an improved radial velocity curve of the Be star through fitting sharp Fe II profiles from the equatorial disk, and a refined Be classification (to that of a B1.5-B2 III star), reveals a black hole of 3.8 to 6.9 solar masses orbiting MWC 656, the candidate counterpart of the gamma-ray source AGL J2241+4454. The black hole is X-ray quiescent and fed by a radiatively inefficient accretion flow giving a luminosity less than 1.6 x 10-7 times the Eddington luminosity. This implies that Be binaries with black-hole companions are difficult to detect by conventional X-ray surveys.

Technical Note: Dissolved organic matter fluorescence a finite mixture approach to deconvolve excitation-emission matrices

The analysis of the shape of excitation-emission matrices (EEMs) is a relevant tool for exploring the origin, transport and fate of dissolved organic matter (DOM) in aquatic ecosystems. Within this context, the decomposition of EEMs is acquiring a notable relevance. A simple mathematical algorithm that automatically deconvolves individual EEMs is described, creating new possibilities for the comparison of DOM fluorescence properties and EEMs that are very different from each other. A mixture model approach is adopted to decompose complex surfaces into sub-peaks. The laplacian operator and the Nelder-Mead optimisation algorithm are implemented to individuate and automatically locate potential peaks in the EEM landscape. The EEMs of a simple artificial mixture of fluorophores and DOM samples collected in a Mediterranean river are used to describe the model application and to illustrate a strategy that optimises the search for the optimal output.

Dynamic and quasistatic trajectories in quasifission reactions and particle emission

We show that the quasifission paths predicted by the one-body dissipation dynamics, in the slowest phase of a binary reaction, follow a quasistatic path, which represents a sequence of states of thermal equilibrium at a fixed value of the deformation coordinate. This establishes the use of the statistical particle-evaporation model in the case of dynamical time-evolving systems. Pre- and post-scission multiplicities of neutrons and total multiplicities of protons and α particles in fission reactions of 63Cu+92Mo, 60Ni+100Mo, 63Cu+100Mo at 10 MeV/u and 20Ne+144,148,154Sm at 20 MeV/u are reproduced reasonably well with statistical model calculations performed along dynamic trajectories whose slow stage (from the most compact configuration up to the point where the neck starts to develop) lasts some 35×10−21 s.

Neutron skin thickness in the droplet model with surface width dependence: Indications of softness of the nuclear symmetry energy

We analyze the neutron skin thickness in finite nuclei with the droplet model and effective nuclear interactions. The ratio of the bulk symmetry energy J to the so-called surface stiffness coefficient Q has in the droplet model a prominent role in driving the size of neutron skins. We present a correlation between the density derivative of the nuclear symmetry energy at saturation and the J/Q ratio. We emphasize the role of the surface widths of the neutron and proton density profiles in the calculation of the neutron skin thickness when one uses realistic mean-field effective interactions. Next, taking as experimental baseline the neutron skin sizes measured in 26 antiprotonic atoms along the mass table, we explore constraints arising from neutron skins on the value of the J/Q ratio. The results favor a relatively soft symmetry energy at subsaturation densities. Our predictions are compared with the recent constraints derived from other experimental observables. Though the various extractions predict different ranges of values, one finds a narrow window L∼45-75 MeV for the coefficient L that characterizes the density derivative of the symmetry energy that is compatible with all the different empirical indications.

Radio Imaging of the Very-High-Energy γ-Ray Emission Region in the Central Engine of a Radio Galaxy

The accretion of matter onto a massive black hole is believed to feed the relativistic plasma jets found in many active galactic nuclei (AGN). Although some AGN accelerate particles to energies exceeding 1012 electron volts and are bright sources of very-high-energy (VHE) γ-ray emission, it is not yet known where the VHE emission originates. Here we report on radio and VHE observations of the radio galaxy Messier 87, revealing a period of extremely strong VHE γ-ray flares accompanied by a strong increase of the radio flux from its nucleus. These results imply that charged particles are accelerated to very high energies in the immediate vicinity of the black hole.

Gamma-ray emission from microquasars: a numerical model for LSI+61º303

We explore the possible association between the microquasar LSI +61°303 and the EGRET source 2CG 135+01/3EG J0241+6103 by studying, with a detailed numerical model, whether this system can produce the emission and the variability detected by EGRET (>100 MeV) through inverse Compton (IC) scattering. Our numerical approach considers a population of relativistic electrons entrained in a cylindrical inhomogeneous jet, interacting with both the radiation and the magnetic fields, taking into account the Thomson and Klein-Nishina regimes of interaction. Our results reproduce the observed spectral characteristics and variability at γ-rays, thus strengthening the identification of LSI +61°303 as a high-energy γ-ray source.

Three-step pathway engineering results in more incidence rate and higher emission of nerolidol and improved attraction of Diadegma semiclausum

The concentration and ratio of terpenoids in the headspace volatile blend of plants have a fundamental role in the communication of plants and insects. The sesquiterpene (E)-nerolidol is one of the important volatiles with effect on beneficial carnivores for biologic pest management in the field. To optimize de novo biosynthesis and reliable and uniform emission of (E)-nerolidol, we engineered different steps of the (E)-nerolidol biosynthesis pathway in Arabidopsis thaliana. Introduction of a mitochondrial nerolidol synthase gene mediates de novo emission of (E)-nerolidol and linalool. Co-expression of the mitochondrial FPS1 and cytosolic HMGR1 increased the number of emitting transgenic plants (incidence rate) and the emission rate of both volatiles. No association between the emission rate of transgenic volatiles and their growth inhibitory effect could be established. (E)-Nerolidol was to a large extent metabolized to non-volatile conjugates.

Incompressibility of neutron-rich matter

The saturation properties of neutron-rich matter are investigated in a relativistic mean-field formalism using two accurately calibrated models: NL3 and FSUGold. The saturation properties density, binding energy per nucleon, and incompressibility coefficient are calculated as a function of the neutron-proton asymmetry α≡(N-Z)/A to all orders in α. Good agreement (at the 10% level or better) is found between these numerical calculations and analytic expansions that are given in terms of a handful of bulk parameters determined at saturation density. Using insights developed from the analytic approach and a general expression for the incompressibility coefficient of infinite neutron-rich matter, i.e., K0(α)=K0+Kτα2+ , we construct a hybrid model with values for K0 and Kτ as suggested by recent experimental findings. Whereas the hybrid model provides a better description of the measured distribution of isoscalar monopole strength in the Sn isotopes relative to both NL3 and FSUGold, it significantly underestimates the distribution of strength in 208Pb. Thus, we conclude that the incompressibility coefficient of neutron-rich matter remains an important open problem.

Nuclear Symmetry Energy Probed by Neutron Skin Thickness of Nuclei

We describe a relation between the symmetry energy coefficients csym(ρ) of nuclear matter and asym(A) of finite nuclei that accommodates other correlations of nuclear properties with the low-density behavior of csym(ρ). Here, we take advantage of this relation to explore the prospects for constraining csym(ρ) of systematic measurements of neutron skin sizes across the mass table, using as example present data from antiprotonic atoms. The found constraints from neutron skins are in harmony with the recent determinations from reactions and giant resonances.

Self-Organising Maps and Correlation Analysis as a Tool to Explore Patterns in Excitation-Emission Matrix Data Sets and to Discriminate Dissolved Organic Matter Fluorescence Components

Dissolved organic matter (DOM) is a complex mixture of organic compounds, ubiquitous in marine and freshwater systems. Fluorescence spectroscopy, by means of Excitation-Emission Matrices (EEM), has become an indispensable tool to study DOM sources, transport and fate in aquatic ecosystems. However the statistical treatment of large and heterogeneous EEM data sets still represents an important challenge for biogeochemists. Recently, Self-Organising Maps (SOM) has been proposed as a tool to explore patterns in large EEM data sets. SOM is a pattern recognition method which clusterizes and reduces the dimensionality of input EEMs without relying on any assumption about the data structure. In this paper, we show how SOM, coupled with a correlation analysis of the component planes, can be used both to explore patterns among samples, as well as to identify individual fluorescence components. We analysed a large and heterogeneous EEM data set, including samples from a river catchment collected under a range of hydrological conditions, along a 60-km downstream gradient, and under the influence of different degrees of anthropogenic impact. According to our results, chemical industry effluents appeared to have unique and distinctive spectral characteristics. On the other hand, river samples collected under flash flood conditions showed homogeneous EEM shapes. The correlation analysis of the component planes suggested the presence of four fluorescence components, consistent with DOM components previously described in the literature. A remarkable strength of this methodology was that outlier samples appeared naturally integrated in the analysis. We conclude that SOM coupled with a correlation analysis procedure is a promising tool for studying large and heterogeneous EEM data sets.

Glassy magnetic phase driven by short range charge and magnetic ordering in nanocrystalline La1/3Sr2/3FeO3-δ: Magnetization, Mössbauer, and polarized neutron studies

The charge ordered La1/3Sr2/3FeO3−δ (LSFO) in bulk and nanocrystalline forms are investigated using ac and dc magnetization, M¨ossbauer, and polarized neutron studies. A complex scenario of short-range charge and magnetic ordering is realized from the polarized neutron studies in nanocrystalline specimen. This short-range ordering does not involve any change in spin state and modification in the charge disproportion between Fe3+ and Fe5+ compared to bulk counterpart as evident in the M¨ossbauer results. The refinement of magnetic diffraction peaks provides magnetic moments of Fe3+ and Fe5+ are about 3.15 μB and 1.57 μB for bulk, and 2.7 μB and 0.53 μB for nanocrystalline specimen, respectively. The destabilization of charge ordering leads to magnetic phase separation, giving rise to the robust exchange bias (EB) effect. Strikingly, EB field at 5 K attains a value as high as 4.4 kOe for average size ∼70 nm, which is zero for the bulk counterpart. A strong frequency dependence of ac susceptibility reveals cluster-glass-like transition around ∼65 K, below which EB appears. Overall results propose that finite-size effect directs the complex glassy magnetic behavior driven by unconventional short-range charge and magnetic ordering, and magnetic phase separation appears in nanocrystalline LSFO.

Analysis of bulk and surface contributions in the neutron skin of nuclei

The neutron skin thickness of nuclei is a sensitive probe of the nuclear symmetry energy and has multiple implications for nuclear and astrophysical studies. However, precision measurements of this observable are difficult to obtain. The analysis of the experimental data may imply some assumptions about the bulk or surface nature of the formation of the neutron skin. Here we study the bulk or surface character of neutron skins of nuclei following from calculations with Gogny, Skyrme, and covariant nuclear mean-field interactions. These interactions are successful in describing nuclear charge radii and binding energies but predict different values for neutron skins. We perform the study by fitting two-parameter Fermi distributions to the calculated self-consistent neutron and proton densities. We note that the equivalent sharp radius is a more suitable reference quantity than the half-density radius parameter of the Fermi distributions to discern between the bulk and surface contributions in neutron skins. We present calculations for nuclei in the stability valley and for the isotopic chains of Sn and Pb.

Gamma-ray emission from massive young stellar objects

Context.Massive stars form in dense and massive molecular cores. The exact formation mechanism is unclear, but it is possible that some massive stars are formed by processes similar to those that produce the low-mass stars, with accretion/ejection phenomena occurring at some point of the evolution of the protostar. This picture seems to be supported by the detection of a collimated stellar wind emanating from the massive protostar IRAS 16547-4247. A triple radio source is associated with the protostar: a compact core and two radio lobes. The emission of the southern lobe is clearly non-thermal. Such emission is interpreted as synchrotron radiation produced by relativistic electrons locally accelerated at the termination point of a thermal jet. Since the ambient medium is determined by the properties of the molecular cloud in which the whole system is embedded, we can expect high densities of particles and infrared photons. Because of the confirmed presence of relativistic electrons, inverse Compton and relativistic Bremsstrahlung interactions are unavoidable. Aims.We aim to make quantitative predictions of the spectral energy distribution of the non-thermal spots generated by massive young stellar objects, with emphasis on the particular case of IRAS 16547-4247. Methods.We study the high-energy emission generated by the relativistic electrons which produce the non-thermal radio source in IRAS 16547-4247. We also study the result of proton acceleration at the terminal shock of the thermal jet and make estimates of the secondary gamma rays and electron-positron pairs produced by pion decay. Results.We present spectral energy distributions for the southern lobe of IRAS 16547-4247, for a variety of conditions. We show that high-energy emission might be detectable from this object in the gamma-ray domain. The source may also be detectable in X-rays through long exposures with current X-ray instruments. Conclusions.Gamma-ray telescopes such as GLAST, and even ground-based Cherenkov arrays of new generation can be used to study non-thermal processes occurring during the formation of massive stars.

Classical line shapes based on analytical solutions of bimolecular trajectories in collision induced emission. II. Reactive collisions

The classical theory of collision induced emission (CIE) from pairs of dissimilar rare gas atoms was developed in Paper I [D. Reguera and G. Birnbaum, J. Chem. Phys. 125, 184304 (2006)] from a knowledge of the straight line collision trajectory and the assumption that the magnitude of the dipole could be represented by an exponential function of the inter-nuclear distance. This theory is extended here to deal with other functional forms of the induced dipole as revealed by ab initio calculations. Accurate analytical expression for the CIE can be obtained by least square fitting of the ab initio values of the dipole as a function of inter-atomic separation using a sum of exponentials and then proceeding as in Paper I. However, we also show how the multi-exponential fit can be replaced by a simpler fit using only two analytic functions. Our analysis is applied to the polar molecules HF and HBr. Unlike the rare gas atoms considered previously, these atomic pairs form stable bound diatomic molecules. We show that, interestingly, the spectra of these reactive molecules are characterized by the presence of multiple peaks. We also discuss the CIE arising from half collisions in excited electronic states, which in principle could be probed in photo-dissociation experiments.