In silico infrared and Raman spectroscopy under pressure: The case of CaSnO3 perovskite

The CaSnO3 perovskite is investigated under geochemical pressure, up to 25 GPa, by means of periodic ab initio calculations performed at B3LYP level with local Gaussian-type orbital basis sets. Structural, elastic, and spectroscopic (phonon wave-numbers, infrared and Raman intensities) properties are fully characterized and discussed. The evolution of the Raman spectrum of CaSnO3 under pressure is reported to remarkably agree with a recent experimental determination [J. Kung, Y. J. Lin, and C. M. Lin, J. Chem. Phys. 135, 224507 (2011)] as regards both wave-number shifts and intensity changes. All phonon modes are symmetry-labeled and bands assigned. The single-crystal total spectrum is symmetry-decomposed into the six directional spectra related to the components of the polarizability tensor. The infrared spectrum at increasing pressure is reported for the first time and its main features discussed. All calculations are performed using the CRYSTAL14 program, taking advantage of the new implementation of analytical infrared and Raman intensities for crystalline materials. (C) 2015 AIP Publishing LLC.

The generator coordinate method in the unrestricted Hartree-Fock formalism

The generator coordinate method was implemented in the unrestricted Hartree-Fock formalism. Weight functions were built from Gaussian generator functions for 1s, 2s, and 2p orbitals of carbon and oxygen atoms. These weight functions show a similar behavior to those found in the generator coordinate restricted Hartree-Fock method, i.e., they are smooth, continuous, and tend to zero in the limits of integration. Moreover, the weight functions obtained are different for spin-up and spin-down electrons what is a result from spin polarization. (C) 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012

The radical SeCl: A theoretical contribution to the characterization of its low-lying electronic states

All doublet and quartet electronic states correlating with the first dissociation channel of SeCl and some Rydberg states are investigated theoretically at the CASSCF/MRCI level of theory using extended basis sets, including the contribution of spin-orbit effects. The similarity of the potential energy curves with those of SeF suggests that spectroscopic constants for the ground (X (2)Pi) and the first excited quartet (a(4)Sigma) of SeCl could also be determined via an emission resulting from the reaction of selenium with atomic chlorine. The coupling constant of the ground state at R-e is estimated as -1610 cm (1). The potential energy curves calculated and the derived spectroscopic constants do not support the interpretation and assignment of the scarce transitions recorded experimentally as due to (2)Pi-(2)Pi emissions. That the few observed lines might arise from transitions from the state b(4)Sigma(-)(1/2) to a very high vibrational level of the state a(4)Sigma(-)(1/2) is an open possibility, however, the number of vibrational states and the calculated Delta G(1/2) differ significantly from the reported ones. (C) 2012 Elsevier B. V. All rights reserved.

Semiempirical Study of the Bergman Reaction: Towards a Computationally Efficient and Accurate Method for Modeling Enediyne Anticancer Antibiotics

Over the past 7 years, the enediyne anticancer antibiotics have been widely studied due to their DNA cleaving ability. The focus of these antibiotics, represented by kedarcidin chromophore, neocarzinostatin chromophore, calicheamicin, esperamicin A, and dynemicin A, is on the enediyne moiety contained within each of these antibiotics. In its inactive form, the moiety is benign to its environment. Upon suitable activation, the system undergoes a Bergman cycloaromatization proceeding through a 1,4-dehydrobenzene diradical intermediate. It is this diradical intermediate that is thought to cleave double-stranded dna through hydrogen atom abstraction. Semiempirical, semiempiricalci, Hartree–Fock ab initio, and mp2 electron correlation methods have been used to investigate the inactive hex-3-ene-1,5-diyne reactant, the 1,4-dehydrobenzene diradical, and a transition state structure of the Bergman reaction. Geometries calculated with different basis sets and by semiempirical methods have been used for single-point calculations using electron correlation methods. These results are compared with the best experimental and theoretical results reported in the literature. Implications of these results for computational studies of the enediyne anticancer antibiotics are discussed.

Rotational constants and structure of para-difluorobenzene determined by femtosecond Raman coherence spectroscopy: A new transient type

Femtosecond Raman rotational coherence spectroscopy (RCS) detected by degenerate four-wave mixing is a background-free method that allows to determine accurate gas-phase rotational constants of non-polar molecules. Raman RCS has so far mostly been applied to the regular coherence patterns of symmetric-top molecules, while its application to nonpolar asymmetric tops has been hampered by the large number of RCS transient types, the resulting variability of the RCS patterns, and the 10³–10⁴ times larger computational effort to simulate and fit rotational Raman RCS transients. We present the rotational Raman RCS spectra of the nonpolar asymmetric top 1,4-difluorobenzene (para-difluorobenzene, p-DFB) measured in a pulsed Ar supersonic jet and in a gas cell over delay times up to ~2.5 ns. p-DFB exhibits rotational Raman transitions with ΔJ = 0, 1, 2 and ΔK = 0, 2, leading to the observation of J −, K −, A −, and C–type transients, as well as a novel transient (S–type) that has not been characterized so far. The jet and gas cell RCS measurements were fully analyzed and yield the ground-state (v = 0) rotational constants Aₒ = 5637.68(20) MHz, Bₒ = 1428.23(37) MHz, and Cₒ = 1138.90(48) MHz (1σ uncertainties). Combining the Aₒ, Bₒ, and Cₒ constants with coupled-cluster with single-, double- and perturbatively corrected triple-excitation calculations using large basis sets allows to determine the semi-experimental equilibrium bond lengths rₑ(C₁–C₂) = 1.3849(4) Å, rₑ(C₂–C³) = 1.3917(4) Å, rₑ(C–F) = 1.3422(3) Å, and rₑ(C₂–H₂) = 1.0791(5) Å.

Gaussian operator bases for correlated fermions

We formulate a general multi-mode Gaussian operator basis for fermions, to enable a positive phase-space representation of correlated Fermi states. The Gaussian basis extends existing bosonic phase-space methods to Fermi systems and thus allows first-principles dynamical or equilibrium calculations in quantum many-body Fermi systems. We prove the completeness of the basis and derive differential forms for products with one- and two-body operators. Because the basis satisfies fermionic superselection rules, the resulting phase space involves only c-numbers, without requiring anticommuting Grassmann variables. Furthermore, because of the overcompleteness of the basis, the phase-space distribution can always be chosen positive. This has important consequences for the sign problem in fermion physics.

Robust automatic mapping algorithms in a network monitoring scenario

Automatically generating maps of a measured variable of interest can be problematic. In this work we focus on the monitoring network context where observations are collected and reported by a network of sensors, and are then transformed into interpolated maps for use in decision making. Using traditional geostatistical methods, estimating the covariance structure of data collected in an emergency situation can be difficult. Variogram determination, whether by method-of-moment estimators or by maximum likelihood, is very sensitive to extreme values. Even when a monitoring network is in a routine mode of operation, sensors can sporadically malfunction and report extreme values. If this extreme data destabilises the model, causing the covariance structure of the observed data to be incorrectly estimated, the generated maps will be of little value, and the uncertainty estimates in particular will be misleading. Marchant and Lark [2007] propose a REML estimator for the covariance, which is shown to work on small data sets with a manual selection of the damping parameter in the robust likelihood. We show how this can be extended to allow treatment of large data sets together with an automated approach to all parameter estimation. The projected process kriging framework of Ingram et al. [2007] is extended to allow the use of robust likelihood functions, including the two component Gaussian and the Huber function. We show how our algorithm is further refined to reduce the computational complexity while at the same time minimising any loss of information. To show the benefits of this method, we use data collected from radiation monitoring networks across Europe. We compare our results to those obtained from traditional kriging methodologies and include comparisons with Box-Cox transformations of the data. We discuss the issue of whether to treat or ignore extreme values, making the distinction between the robust methods which ignore outliers and transformation methods which treat them as part of the (transformed) process. Using a case study, based on an extreme radiological events over a large area, we show how radiation data collected from monitoring networks can be analysed automatically and then used to generate reliable maps to inform decision making. We show the limitations of the methods and discuss potential extensions to remedy these.

Modern Computational Chemistry Methods for Prediction of Ground- and Excited-State Properties in Open-Shell Systems

The main goal of the research presented in this work is to provide some important insights about computational modeling of open-shell species. Such projects are: the investigation of the size-extensivity error in Equation-of-Motion Coupled Cluster methods, the analysis of the Long-Range corrected scheme in predicting UV-Vis spectra of Cu(II) complexes with the 4-imidazole acetate and its ethylated derivative, and the exploration of the importance of choosing a proper basis set for the description of systems such as the lithium monoxide anion. The most significant findings of this research are: (i) The contribution of the left operator to the size-extensivity error of the CR-EOMCC(2,3) approach, (ii) The cause of d-d shifts when varying the range-separation parameter and the amount of the exact exchange arising from the imbalanced treatment of localized vs. delocalized orbitals via the "tuned" CAM-B3LYP* functional, (iii) The proper acidity trend of the first-row hydrides and their lithiated analogs that may be reversed if the basis sets are not correctly selected.

Core-valence Correlation Effects On Ir Calculations: The Bf3 And Bcl3 Cases.

The first theoretical results of core-valence correlation effects are presented for the infrared wavenumbers and intensities of the BF3 and BCl3 molecules, using (double- and triple-zeta) Dunning core-valence basis sets at the CCSD(T) level. The results are compared with those calculated in the frozen core approximation with standard Dunning basis sets at the same correlation level and with the experimental values. The general conclusion is that the effect of core-valence correlation is, for infrared wavenumbers and intensities, smaller than the effect of adding augmented diffuse functions to the basis set, e.g., cc-pVTZ to aug-cc-pVTZ. Moreover, the trends observed in the data are mainly related to the augmented functions rather than the core-valence functions added to the basis set. The results obtained here confirm previous studies pointing out the large descrepancy between the theoretical and experimental intensities of the stretching mode for BCl3.

Analyzing the n→π* Electronic Transition of Formaldehyde in Water. A Sequential Monte Carlo/Time-Dependent Density Functional Theory

The n→π* absorption transition of formaldehyde in water is analyzed using combined and sequential classical Monte Carlo (MC) simulations and quantum mechanics (QM) calculations. MC simulations generate the liquid solute-solvent structures for subsequent QM calculations. Using time-dependent density functional theory in a localized set of gaussian basis functions (TD-DFT/6-311++G(d,p)) calculations are made on statistically relevant configurations to obtain the average solvatochromic shift. All results presented here use the electrostatic embedding of the solvent. The statistically converged average result obtained of 2300 cm-1 is compared to previous theoretical results available. Analysis is made of the effective dipole moment of the hydrogen-bonded shell and how it could be held responsible for the polarization of the solvent molecules in the outer solvation shells.

Systematic investigation of a family of gradient-dependent functionals for solids

Eleven density functionals are compared with regard to their performance for the lattice constants of solids. We consider standard functionals, such as the local-density approximation and the Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA), as well as variations of PBE GGA, such as PBEsol and similar functionals, PBE-type functionals employing a tighter Lieb-Oxford bound, and combinations thereof. On a test set of 60 solids, we perform a system-by-system analysis for selected functionals and a full statistical analysis for all of them. The impact of restoring the gradient expansion and of tightening the Lieb-Oxford bound is discussed, and confronted with previous results obtained from other codes, functionals or test sets. No functional is uniformly good for all investigated systems, but surprisingly, and pleasingly, the simplest possible modifications to PBE turn out to have the most beneficial effect on its performance. The atomization energy of molecules was also considered and on a testing set of six molecules, we found that the PBE functional is clearly the best, the others leading to strong overbinding.

Assessment of body composition by bioelectrical impedance analysis without the need for measurement of height

Conventional whole-body single frequency bioelectrical impedance analysis (BIA) of body composition typically uses height as a surrogate measure of conductor length. A new method of BIA analysis for the prediction of body cell mass (BCM) and extracellular water (ECW, as % body weight) not using height has been introduced-the Soft Tissue Analyser (STA(TM), Akern Sri, Florence, Italy)-making it ideal for use in subjects where measurement of height is difficult or impossible. The performance of the new analytical method in predicting BCM and ECW in 139 normal control subjects was assessed by comparison with reference data obtained from a four-component (4-C) model of body composition and with predictions obtained from conventional BIA analysis. Both predicted BCM and ECW were strongly (r = 0.82, SEE = 6.3 kg and 0.89, SEE = 1.3 kg respectively) correlated with the corresponding 4-C model measurements although differing significantly from the lines of identity (P < 0.0001). Fat-free mass, calculated from STA estimates of BCM and ECW, was better predicted (r = 0.91, SEE = 5.6 kg). The significant differences in STA-group mean values for BCM and ECW and wide limits of agreement compared with the reference data indicate that the method cannot be used with confidence for prediction of these body compartments despite the obvious advantage of not requiring an accurate measurement of height. (C) 2001 Harcourt Publishers Ltd.

Effective atomic orbitals for fuzzy atoms

The method of extracting effective atomic orbitals and effective minimal basis sets from molecular wave function characterizing the state of an atom in a molecule is developed in the framework of the "fuzzy" atoms. In all cases studied, there were as many effective orbitals that have considerable occupation numbers as orbitals in the classical minimal basis. That is considered to be of high conceptual importance

Estimates of electronic coupling for excess electron transfer in DNA

Electronic coupling Vda is one of the key parameters that determine the rate of charge transfer through DNA. While there have been several computational studies of Vda for hole transfer, estimates of electronic couplings for excess electron transfer (ET) in DNA remain unavailable. In the paper, an efficient strategy is established for calculating the ET matrix elements between base pairs in a π stack. Two approaches are considered. First, we employ the diabatic-state (DS) method in which donor and acceptor are represented with radical anions of the canonical base pairs adenine-thymine (AT) and guanine-cytosine (GC). In this approach, similar values of Vda are obtained with the standard 6-31 G* and extended 6-31++ G* basis sets. Second, the electronic couplings are derived from lowest unoccupied molecular orbitals (LUMOs) of neutral systems by using the generalized Mulliken-Hush or fragment charge methods. Because the radical-anion states of AT and GC are well reproduced by LUMOs of the neutral base pairs calculated without diffuse functions, the estimated values of Vda are in good agreement with the couplings obtained for radical-anion states using the DS method. However, when the calculation of a neutral stack is carried out with diffuse functions, LUMOs of the system exhibit the dipole-bound character and cannot be used for estimating electronic couplings. Our calculations suggest that the ET matrix elements Vda for models containing intrastrand thymine and cytosine bases are essentially larger than the couplings in complexes with interstrand pyrimidine bases. The matrix elements for excess electron transfer are found to be considerably smaller than the corresponding values for hole transfer and to be very responsive to structural changes in a DNA stack

Accurate calculations of intermolecular interaction energies using explicitly correlated coupled cluster wave functions and a dispersion-weighted MP2 method.

Explicitly correlated coupled-cluster calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented. Results obtained with the recently proposed CCSD(T)-F12a method and augmented double-zeta basis sets are found to be in very close agreement with basis set extrapolated conventional CCSD(T) results. Furthermore, we propose a dispersion-weighted MP2 (DW-MP2) approximation that combines the good accuracy of MP2 for complexes with predominately electrostatic bonding and SCS-MP2 for dispersion-dominated ones. The MP2-F12 and SCS-MP2-F12 correlation energies are weighted by a switching function that depends on the relative HF and correlation contributions to the interaction energy. For the S22 set, this yields a mean absolute deviation of 0.2 kcal/mol from the CCSD(T)-F12a results. The method, which allows obtaining accurate results at low cost, is also tested for a number of dimers that are not in the training set.