Predicting spectral features in galaxy spectra from broad-band photometry

We explore the prospects of predicting emission-line features present in galaxy spectra given broad-band photometry alone. There is a general consent that colours, and spectral features, most notably the 4000 angstrom break, can predict many properties of galaxies, including star formation rates and hence they could infer some of the line properties. We argue that these techniques have great prospects in helping us understand line emission in extragalactic objects and might speed up future galaxy redshift surveys if they are to target emission-line objects only. We use two independent methods, Artificial Neural Networks (based on the ANNz code) and Locally Weighted Regression (LWR), to retrieve correlations present in the colour N-dimensional space and to predict the equivalent widths present in the corresponding spectra. We also investigate how well it is possible to separate galaxies with and without lines from broad-band photometry only. We find, unsurprisingly, that recombination lines can be well predicted by galaxy colours. However, among collisional lines some can and some cannot be predicted well from galaxy colours alone, without any further redshift information. We also use our techniques to estimate how much information contained in spectral diagnostic diagrams can be recovered from broad-band photometry alone. We find that it is possible to classify active galactic nuclei and star formation objects relatively well using colours only. We suggest that this technique could be used to considerably improve redshift surveys such as the upcoming Fibre Multi Object Spectrograph (FMOS) survey and the planned Wide Field Multi Object Spectrograph (WFMOS) survey.

Alternative diagnostic diagrams and the `forgotten` population of weak line galaxies in the SDSS

A numerous population of weak line galaxies (WLGs) is often left out of statistical studies on emission-line galaxies (ELGs) due to the absence of an adequate classification scheme, since classical diagnostic diagrams, such as [O iii]/H beta versus [N ii]/H alpha (the BPT diagram), require the measurement of at least four emission lines. This paper aims to remedy this situation by transposing the usual divisory lines between star-forming (SF) galaxies and active galactic nuclei (AGN) hosts and between Seyferts and LINERs to diagrams that are more economical in terms of line quality requirements. By doing this, we rescue from the classification limbo a substantial number of sources and modify the global census of ELGs. More specifically, (1) we use the Sloan Digital Sky Survey Data Release 7 to constitute a suitable sample of 280 000 ELGs, one-third of which are WLGs. (2) Galaxies with strong emission lines are classified using the widely applied criteria of Kewley et al., Kauffmann et al. and Stasinska et al. to distinguish SF galaxies and AGN hosts and Kewley et al. to distinguish Seyferts from LINERs. (3) We transpose these classification schemes to alternative diagrams keeping [N ii]/H alpha as a horizontal axis, but replacing H beta by a stronger line (H alpha or [O ii]), or substituting the ionization-level sensitive [O iii]/H beta ratio with the equivalent width of H alpha (W(H alpha)). Optimized equations for the transposed divisory lines are provided. (4) We show that nothing significant is lost in the translation, but that the new diagrams allow one to classify up to 50 per cent more ELGs. (5) Introducing WLGs in the census of galaxies in the local Universe increases the proportion of metal-rich SF galaxies and especially LINERs. In the course of this analysis, we were led to make the following points. (i) The Kewley et al. BPT line for galaxy classification is generally ill-used. (ii) Replacing [O iii]/H beta by W(H alpha) in the classification introduces a change in the philosophy of the distinction between LINERs and Seyferts, but not in its results. Because the W(H alpha) versus [N ii]/H alpha diagram can be applied to the largest sample of ELGs without loss of discriminating power between Seyferts and LINERs, we recommend its use in further studies. (iii) The dichotomy between Seyferts and LINERs is washed out by WLGs in the BPT plane, but it subsists in other diagnostic diagrams. This suggests that the right wing in the BPT diagram is indeed populated by at least two classes, tentatively identified with bona fide AGN and `retired` galaxies that have stopped forming stars and are ionized by their old stellar populations.

The evolution of the mass-metallicity relation in SDSS galaxies uncovered by astropaleontology

We have obtained the mass-metallicity (M-Z) relation at different lookback times for the same set of galaxies from the Sloan Digital Sky Survey, using the stellar metallicities estimated with our spectral synthesis code STARLIGHT. We have found that this relation steepens and spans a wider range in both mass and metallicity at higher redshifts. We have modelled the time evolution of stellar metallicity with a closed-box chemical evolution model, for galaxies of different types and masses. Our results suggest that the M-Z relation for galaxies with present-day stellar masses down to 10(10) M(circle dot) is mainly driven by the history of star formation and not by inflows or outflows.

Analysis of drop formation at conical tips : 1. Theory

The phenomenon of drop formation at conical tips under near zero flow conditions has been investigated using a theoretical approach. The analysis permits the prediction of drop profile and drop volume, until the onset of instability. A semiempirical approach based on the similarity of drop shapes has been adopted to predict the detaching drop volumes at conical tips. The effects of base diameter of the cone, cone angle, interfacial tension, and the densities of the drop and the surrounding fluid on the maximum and detached drop volumes are predicted.

Conformational analysis of D-pantothenic acid

The conformational analysis of d-pantothenic acid using classical semiempirical methods has been carried out. The pantothenic acid molecule can exist in the neutral form (I) or in the ionised form (II) with a deprotonated negatively charged carboxyl group. The neutral molecule as well as the anion is highly flexible and has an ensemble of several allowed conformations rather than one or two unique conformations. The distribution of allowed conformations indicate that the β-alanine as well as the pantoic acid part of the molecule prefers partially folded conformations. The conformation of the former is greatly affected by the ionisation state of the carboxyl group whereas that of the latter is not. Possibility of intramolecular hydrogen bonding in different allowed conformations has also been explored. A bifurcated hydrogen bond involving a carboxyl (or carboxylate) oxygen atom and a hydroxyl oxygen atom, as acceptors, and the amide nitrogen atom as the donor occurs frequently in both I and II. Amongst the two crystal structures containing pantothenic acid reported so far, the conformation of the molecule in l-lysine d-pantothenate lies in the allowed region and is stabilised by a bifurcated intramolecular hydrogen bond, whereas that in the calcium bromide salt falls in a disallowed region, presumably due to the requirement of tridentate metal coordination.

Conformational analysis of 2,2 '-bithiophene under interaction of external electric field defined by point charges

Conformational analysis of 2,2'-bithiophene (BT) under the influence of an electric field (EF) constructed by point charges has been performed by using semi-empirical Austin Model 1 (AM1) and Parametric model number 3 (PM3) calculations. When the EF perpendicular to the molecular conjugation chain is applied, both AM1 and PM3 calculations show an energy increase of the anti-conformation. AM1 predicts that the global minimum shifts to syn-conformation when the EF strength is larger than a critical value. and PM predicts that the local minimum in anti-conformation vanishes. This kind of EF effect has been ascribed to the EF and dipole moment interaction.

Chemical bond analysis of nonlinearity of urea crystal

A novel and quantitative study on structure-property relationships has been carried out in urea crystal, based on the dielectric theory of complex crystals and the modified Levine bond charge model, mainly from the chemical bond viewpoint. For the first time, it was treated like this, and the bond parameters and linear and nonlinear characteristics of constituent chemical bonds were presented quantitatively. The theoretical result agrees satisfactorily with the experimental datum and can reasonably explain the nonlinear origin of urea crystal, that is, the C-N bond in the conjugated system of bonds O double left arrow C<--N-H. At the same time, the novel method should be a useful tool toward the future development of the search for new nonlinear optical (NLO) materials in the organic crystal field.

Analysis of the electronic structure, hyperfine structure, and volume isotope shifts in the low lying states of BA I and BA II

Relativistic multi-configuration Dirac Fock (MCDF) wavefunctions coupled to good angular momentum J have been calculated for low lying states of Ba I and Ba II. These wavefunctions are compared with semiempirical ones derived from experimental atomic energy levels. It is found that significantly better agreement is obtained when close configurations are included in the MCDF wavefunctions. Calculations of the electronic part of the field isotope shift lead to very good agreement with electronic factors derived from experimental data. Furthermore, the slopes of the lines in a King plot analysis of many of the optical lines are predicted accurately by these calculations. However, the MCDF wavefunctions seem not to be of sufficient accuracy to give agreement with the experimental magnetic dipole and electric quadrupole hyperfine structure constants.

Theoretical Analysis of Aggregation in Block-Copolymer Films: The Optical Signature

We focus this work on the theoretical investigation of the block-copolymer poly [oxyoctyleneoxy-(2,6-dimethoxy-1,4phenylene-1,2-ethinylene-phenanthrene-2,4diyl) named as LaPPS19, recently proposed for optoelectronic applications. We used for that a variety of methods, from molecular mechanics to quantum semiempirical techniques (AMI, ZINDO/S-CIS). Our results show that as expected isolated LaPPS19 chains present relevant electron localization over the phenanthrene group. We found, however, that LaPPS19 could assemble in a pi-stacked form, leading to impressive interchain interaction; the stacking induces electronic delocalization between neighbor chains and introduces new states below the phenanthrene-related absorption; these results allowed us to associate the red-shift of the absorption edge, seen in the experimental results, to spontaneous pi-stack aggregation of the chains. (C) 2009 Wiley Periodicals, Inc. Int J Quantum Chem 110: 885-892, 2010

New insights on surface-enhanced Raman scattering based on controlled aggregation and spectroscopic studies, DFT calculations and symmetry analysis for 3,6-bi-2-pyridyl-1,2,4,5-tetrazine adsorbed onto citrate-stabilized gold nanoparticles

Asystematic study on the surface-enhanced Raman scattering (SERS) for 3,6-bi-2-pyridyl-1,2,4,5-tetrazine (bptz) adsorbed onto citrate-modified gold nanoparticles (cit-AuNps) was carried out based on electronic and vibrational spectroscopy and density functional methods. The citrate/bptz exchange was carefully controlled by the stepwise addition of bptz to the cit-AuNps, inducing flocculation and leading to the rise of a characteristic plasmon coupling band in the visible region. Such stepwise procedure led to a uniform decrease of the citrate SERS signals and to the rise of characteristic peaks of bptz, consistent with surface binding via the N heterocyclic atoms. In contrast, single addition of a large amount of bptz promoted complete aggregation of the nanoparticles, leading to a strong enhancement of the SERS signals. In this case, from the distinct Raman profiles involved, the formation of a new SERS environment became apparent, conjugating the influence of the local hot spots and charge-transfer (CT) effects. The most strongly enhanced vibrations belong to a(1) and b(2) representations, and were interpreted in terms of the electromagnetic and the CT mechanisms: the latter involving significant contribution of vibronic coupling in the system. Copyright (C) 2010 John Wiley & Sons, Ltd.

Computation of (3)J(HH) coupling constants with a combination of density functional theory and semiempirical calculations. Application to complex molecules

This work evaluates the efficiency of economic levels of theory for the prediction of (3)J(HH) spin-spin coupling constants, to be used when robust electronic structure methods are prohibitive. To that purpose, DFT methods like mPW1PW91. B3LYP and PBEPBE were used to obtain coupling constants for a test set whose coupling constants are well known. Satisfactory results were obtained in most of cases, with the mPW1PW91/6-31G(d,p)//B3LYP/6-31G(d,p) leading the set. In a second step. B3LYP was replaced by the semiempirical methods PM6 and RM1 in the geometry optimizations. Coupling constants calculated with these latter structures were at least as good as the ones obtained by pure DFT methods. This is a promising result, because some of the main objectives of computational chemistry - low computational cost and time, allied to high performance and precision - were attained together. (C) 2012 Elsevier B.V. All rights reserved.

Structural Analysis of alpha-Fetoprotein (AFP)-like Peptides with Anti-Breast-Cancer Properties

The abundance of alpha-fetoprotein (AFP), a natural protein produced by the fetal yolk sac during pregnancy, correlates with lower incidence of estrogen receptor positive (ER+) breast cancer. The pharmacophore region of AFP has been narrowed down to a four amino acid (AA) region in the third domain of the 591 AA peptide. Our computational study focuses on a 4-mer segment consisting of the amino acids threonine-proline-valine-asparagine (TPVN). We have run replica exchange molecular dynamics (REMD) simulations and used 120 configurational snapshots from the total trajectory as starting configurations for quantum chemical calculations. We optimized structures using semiempirical (PM3, PM6, PM6-D2, PM6-H2, PM6-DH+, PM6-DH2) and density functional methods (TPSS, PBE0, M06-2X). By comparing the accuracy of these methods against RI-MP2 benchmarks, we devised a protocol for calculating the lowest energy conformers of these peptides accurately and efficiently. This protocol screens out high-energy conformers using lower levels of theory and outlines a general method for predicting small peptide structures.

Assessment of interatomic potentials for atomistic analysis of static and dynamic properties of screw dislocations in W

Screw dislocations in bcc metals display non-planar cores at zero temperature which result in high lattice friction and thermally-activated strain rate behavior. In bcc W, electronic structure molecular statics calculations reveal a compact, non-degenerate core with an associated Peierls stress between 1.7 and 2.8 GPa. However, a full picture of the dynamic behavior of dislocations can only be gained by using more efficient atomistic simulations based on semiempirical interatomic potentials. In this paper we assess the suitability of five different potentials in terms of static properties relevant to screw dislocations in pure W. Moreover, we perform molecular dynamics simulations of stress-assisted glide using all five potentials to study the dynamic behavior of screw dislocations under shear stress. Dislocations are seen to display thermally-activated motion in most of the applied stress range, with a gradual transition to a viscous damping regime at high stresses. We find that one potential predicts a core transformation from compact to dissociated at finite temperature that affects the energetics of kink-pair production and impacts the mechanism of motion. We conclude that a modified embedded-atom potential achieves the best compromise in terms of static and dynamic screw dislocation properties, although at an expense of about ten-fold compared to central potentials.

Raman microscopy of formamide-intercalated kaolinites treated by controlled-rate thermal analysis technology

Raman spectroscopy of formamide-intercalated kaolinites treated using controlled-rate thermal analysis technology (CRTA), allowing the separation of adsorbed formamide from intercalated formamide in formamide-intercalated kaolinites, is reported. The Raman spectra of the CRTA-treated formamide-intercalated kaolinites are significantly different from those of the intercalated kaolinites, which display a combination of both intercalated and adsorbed formamide. An intense band is observed at 3629 cm-1, attributed to the inner surface hydroxyls hydrogen bonded to the formamide. Broad bands are observed at 3600 and 3639 cm-1, assigned to the inner surface hydroxyls, which are hydrogen bonded to the adsorbed water molecules. The hydroxyl-stretching band of the inner hydroxyl is observed at 3621 cm-1 in the Raman spectra of the CRTA-treated formamide-intercalated kaolinites. The results of thermal analysis show that the amount of intercalated formamide between the kaolinite layers is independent of the presence of water. Significant differences are observed in the CO stretching region between the adsorbed and intercalated formamide.

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.