Relativistic mean-field interaction with density-dependent meson-nucleon vertices based on microscopical calculations

Although ab initio calculations of relativistic Brueckner theory lead to large scalar isovector fields in nuclear matter, at present, successful versions of covariant density functional theory neglect the interactions in this channel. A new high-precision density functional DD-MEδ is presented which includes four mesons, σ, ω, δ, and ρ, with density-dependent meson-nucleon couplings. It is based to a large extent on microscopic ab initiocalculations in nuclear matter. Only four of its parameters are determined by adjusting to binding energies and charge radii of finite nuclei. The other parameters, in particular the density dependence of the meson-nucleon vertices, are adjusted to nonrelativistic and relativistic Brueckner calculations of symmetric and asymmetric nuclear matter. The isovector effective mass mp*−mn* derived from relativistic Brueckner theory is used to determine the coupling strength of the δ meson and its density dependence.

Giant monopole energies from a constrained relativistic mean-field approach

Background:Average energies of nuclear collective modes may be efficiently and accurately computed using a nonrelativistic constrained approach without reliance on a random phase approximation (RPA). Purpose: To extend the constrained approach to the relativistic domain and to establish its impact on the calibration of energy density functionals. Methods: Relativistic RPA calculations of the giant monopole resonance (GMR) are compared against the predictions of the corresponding constrained approach using two accurately calibrated energy density functionals. Results: We find excellent agreement at the 2% level or better between the predictions of the relativistic RPA and the corresponding constrained approach for magic (or semimagic) nuclei ranging from 16 O to 208 Pb. Conclusions: An efficient and accurate method is proposed for incorporating nuclear collective excitations into the calibration of future energy density functionals.

Counterpoise-corrected geometries and harmonic frequencies of N-body clusters: application to (HF)n (n=3,4)

A study was conducted on the methods of basis set superposition error (BSSE)-free geometry optimization and frequency calculations in clusters larger than a dimer. In particular, three different counterpoise schemes were critically examined. It was shown that the counterpoise-corrected supermolecule energy can be easily obtained in all the cases by using the many-body partitioning of energy

Metastable liquid-liquid phase transition in a single-component system with only one crystal phase and no density anomaly

We investigate the phase behavior of a single-component system in three dimensions with spherically-symmetric, pairwise-additive, soft-core interactions with an attractive well at a long distance, a repulsive soft-core shoulder at an intermediate distance, and a hard-core repulsion at a short distance, similar to potentials used to describe liquid systems such as colloids, protein solutions, or liquid metals. We showed [Nature (London) 409, 692 (2001)] that, even with no evidence of the density anomaly, the phase diagram has two first-order fluid-fluid phase transitions, one ending in a gas¿low-density-liquid (LDL) critical point, and the other in a gas¿high-density-liquid (HDL) critical point, with a LDL-HDL phase transition at low temperatures. Here we use integral equation calculations to explore the three-parameter space of the soft-core potential and perform molecular dynamics simulations in the interesting region of parameters. For the equilibrium phase diagram, we analyze the structure of the crystal phase and find that, within the considered range of densities, the structure is independent of the density. Then, we analyze in detail the fluid metastable phases and, by explicit thermodynamic calculation in the supercooled phase, we show the absence of the density anomaly. We suggest that this absence is related to the presence of only one stable crystal structure.

Inflation dynamics and the New Keynesian Phillips curve in EU-4

The paper seeks to shed light on inflation dynamics of four new EU member states: the Czech Republic, Hungary, Poland and Slovakia. To this end, the New Keynesian Phillips curve augmented for open economies is estimated and additional statistical tests applied. We find the following. (1) The claim of New Keynesians that the real marginal cost is the main inflation-forcing variable is fragile. (2) Inflation seems to be driven by external factors. (3) Although inflation holds a forward-looking component, the backward-looking component is substantial. An intuitive explanation for higher inflation persistence may be rather adaptive than rational price setting of local firms.

Detecció, identificació de Mycobacterium tuberculosis i caracterització molecular de la resistència als principals antituberculosos en mostra clínica directa mitjançant piroseqüenciació: aplicabilitat en l’estudi de la clonalitat de les soques.

L’augment de soques de Mycobacerium tuberculosis multiresistents i l’aparició de soques extremadament resistents, ha posat de manifest la necessitat de disposar de nous mètodes de detecció precoç de M.tuberculosis i les seves resistències als principals fàrmacs antituberculosos, així com l’estudi eficaç dels brots, per a poder dur a terme un control eficient de la malaltia. L’objectiu de l’estudi va ser determinar la capacitat de la piroseqüenciació per a detectar i identificar micobacteris a partir d’aïllats clínics i mostra directa, determinar el seu patró de resistència a Rifampicina, Isoniacida, Etambutol, Fluoroquinolones, Canamicina i Capreomicina i explorar la seva potencialitat per a ser emprada en el tipatge molecular de les soques. Mitjançant la optimització de la tecnologia de la piroseqüenciació pels gens analitzats, i el disseny de primers específics, hem verificat la utilitat de la piroseqüenciació per al maneig de la tuberculosi. S’ha pogut estudiar la presència de mutacions relacionades amb resistències a fàrmacs de primera línia en el 96.5% de les posicions analitzades en soca clínica i en el 70.4% en mostra directa. Van concordar amb el resultats obtinguts per altres mètodes fenotípics i/o fenotípics el 97.1% i el 98.2% del resultats obtinguts en soca clínica i mostra directa respectivament. La piroseqüenciació ens ha permet analitzar la presència de mutacions relacionades amb resistències a antituberculosos de segona línia, servint com a mètode de confirmació en l’anàlisi d’altres mètodes genotípics. La tècnica ens ha permès també identificar els principals micobacteris no tuberculosos implicats en la infecció humana, sòls o en coinfecció. Resultats preliminars mostren que l’anàlisi amb piroseqüenciació pot ser d’utilitat per l’estudi clonal de les soques de M.tuberculosis. Les nostres observacions mostren la piroseqüenciació com una eina valuosa en l’àmbit clínic, ja que permet una reducció de la demora del diagnòstic, estudi filogenètic i detecció de resistències M.tuberculosis, i per tant una millora de l’aplicació del tractament adequat, ajudant al control de la malaltia.

Paper de l’ Haemophilus influenzae en les exacerbacions de la Malaltia Pulmonar Obstructiva Crònica (MPOC): tipificació molecular i factors de virulència

Objetivos: 1.-Identificar los factores clínicos y microbiológicos que ayuden a predecir la aparición de exacerbaciones en la EPOC. 2.-Diagnóstico y cuantificación de las especies bacterianas aisladas en esputo (fase de exacerbación y estable) .3.- Tipificación genotípica secuencial de las cepas de H. influenzae y P. aeruginosa. 4.- Impacto del tratamiento antibiótico en la aparición de resistencias en estos patógenos. 5.- Diseño: Estudio prospectivo (3 años). Ámbito del estudio: Hospital Universitario de tercer nivel. Pacientes con EPOC grave atendidos en la Consulta Monográfica de EPOC del Servicio de Neumología. Métodos microbiológicos: Cuantificación de la carga bacteriana en muestras respiratorias en fase estable y en exacerbación. Estudio de la sensibilidad “in vitro”. Tipificación molecular (PFGE y MLST) de H. influenzae y P. aeruginosa. Estudio de los genes de virulencia de H. influenzae mediante PCR. Resultados: Desde Febrero de 2010 a Julio de 2011 se han incluido 77 pacientes. Los microorganismos más frecuentemente aislados en fase de exacerbación fueron: P. aeruginosa (29.3%), H. influenzae (15.92%), M. catarrhalis (12.74%), S. pneumoniae (10.19%) y S. aureus (5.10%). En los 88 episodios por P. aeruginosa se detectaron 38 genotipos diferentes. En los 41 episodios por H. influenzae se detectaron 39 genotipos diferentes. El 10% de los episodios fueron polimicrobianos. Los episodios de EAEPOC y de fase estable tuvieron una distribución de microorganismos similar. Sin embargo, cuando se cuantificaron las cargas bacterianas fueron mayores en EAEPOC (intervalo 4x107 -2x108) que en fase estable (intervalo 2x105 -4x107). Conclusiones: El genotipo de las cepas de P. aeruginosa y H. influenzae aisladas en EAEPOC difieren de un paciente a otro, sin embargo la mayoría de los episodios de cada paciente están causados por un genotipo único.

Ab initio absorption spectrum of NiO combining molecular dynamics with the embedded cluster approach in a discrete reaction field

We developed a procedure that combines three complementary computational methodologies to improve the theoretical description of the electronic structure of nickel oxide. The starting point is a Car-Parrinello molecular dynamics simulation to incorporate vibrorotational degrees of freedom into the material model. By means ofcomplete active space self-consistent field second-order perturbation theory (CASPT2) calculations on embedded clusters extracted from the resulting trajectory, we describe localized spectroscopic phenomena on NiO with an efficient treatment of electron correlation. The inclusion of thermal motion into the theoretical description allowsus to study electronic transitions that, otherwise, would be dipole forbidden in the ideal structure and results in a natural reproduction of the band broadening. Moreover, we improved the embedded cluster model by incorporating self-consistently at the complete active space self-consistent field (CASSCF) level a discrete (or direct) reaction field (DRF) in the cluster surroundings. The DRF approach offers an efficient treatment ofelectric response effects of the crystalline embedding to the electronic transitions localized in the cluster. We offer accurate theoretical estimates of the absorption spectrum and the density of states around the Fermi level of NiO, and a comprehensive explanation of the source of the broadening and the relaxation of the charge transferstates due to the adaptation of the environment

Selected configuration interaction with truncation energy error and application to the Ne atom

Selected configuration interaction (SCI) for atomic and molecular electronic structure calculations is reformulated in a general framework encompassing all CI methods. The linked cluster expansion is used as an intermediate device to approximate CI coefficients BK of disconnected configurations (those that can be expressed as products of combinations of singly and doubly excited ones) in terms of CI coefficients of lower-excited configurations where each K is a linear combination of configuration-state-functions (CSFs) over all degenerate elements of K. Disconnected configurations up to sextuply excited ones are selected by Brown's energy formula, ΔEK=(E-HKK)BK2/(1-BK2), with BK determined from coefficients of singly and doubly excited configurations. The truncation energy error from disconnected configurations, Δdis, is approximated by the sum of ΔEKS of all discarded Ks. The remaining (connected) configurations are selected by thresholds based on natural orbital concepts. Given a model CI space M, a usual upper bound ES is computed by CI in a selected space S, and EM=E S+ΔEdis+δE, where δE is a residual error which can be calculated by well-defined sensitivity analyses. An SCI calculation on Ne ground state featuring 1077 orbitals is presented. Convergence to within near spectroscopic accuracy (0.5 cm-1) is achieved in a model space M of 1.4× 109 CSFs (1.1 × 1012 determinants) containing up to quadruply excited CSFs. Accurate energy contributions of quintuples and sextuples in a model space of 6.5 × 1012 CSFs are obtained. The impact of SCI on various orbital methods is discussed. Since ΔEdis can readily be calculated for very large basis sets without the need of a CI calculation, it can be used to estimate the orbital basis incompleteness error. A method for precise and efficient evaluation of ES is taken up in a companion paper

Basis set superposition error-counterpoise corrected potential energy surfaces : application to hydrogen peroxide ... X (X=F-, Cl-, Br-, Li+, Na+) complexes

Møller-Plesset (MP2) and Becke-3-Lee-Yang-Parr (B3LYP) calculations have been used to compare the geometrical parameters, hydrogen-bonding properties, vibrational frequencies and relative energies for several X- and X+ hydrogen peroxide complexes. The geometries and interaction energies were corrected for the basis set superposition error (BSSE) in all the complexes (1-5), using the full counterpoise method, yielding small BSSE values for the 6-311 + G(3df,2p) basis set used. The interaction energies calculated ranged from medium to strong hydrogen-bonding systems (1-3) and strong electrostatic interactions (4 and 5). The molecular interactions have been characterized using the atoms in molecules theory (AIM), and by the analysis of the vibrational frequencies. The minima on the BSSE-counterpoise corrected potential-energy surface (PES) have been determined as described by S. Simón, M. Duran, and J. J. Dannenberg, and the results were compared with the uncorrected PES

Electronic couplings and on-site energies for hole transfer in DNA: systematic quantum mechanical/molecular dynamic study

The electron hole transfer (HT) properties of DNA are substantially affected by thermal fluctuations of the π stack structure. Depending on the mutual position of neighboring nucleobases, electronic coupling V may change by several orders of magnitude. In the present paper, we report the results of systematic QM/molecular dynamic (MD) calculations of the electronic couplings and on-site energies for the hole transfer. Based on 15 ns MD trajectories for several DNA oligomers, we calculate the average coupling squares 〈 V2 〉 and the energies of basepair triplets X G+ Y and X A+ Y, where X, Y=G, A, T, and C. For each of the 32 systems, 15 000 conformations separated by 1 ps are considered. The three-state generalized Mulliken-Hush method is used to derive electronic couplings for HT between neighboring basepairs. The adiabatic energies and dipole moment matrix elements are computed within the INDO/S method. We compare the rms values of V with the couplings estimated for the idealized B -DNA structure and show that in several important cases the couplings calculated for the idealized B -DNA structure are considerably underestimated. The rms values for intrastrand couplings G-G, A-A, G-A, and A-G are found to be similar, ∼0.07 eV, while the interstrand couplings are quite different. The energies of hole states G+ and A+ in the stack depend on the nature of the neighboring pairs. The X G+ Y are by 0.5 eV more stable than X A+ Y. The thermal fluctuations of the DNA structure facilitate the HT process from guanine to adenine. The tabulated couplings and on-site energies can be used as reference parameters in theoretical and computational studies of HT processes in DNA

Two complementary molecular energy decomposition schemes: the Mayer and Ziegler-Rauk methods in comparison

In the present paper we discuss and compare two different energy decomposition schemes: Mayer's Hartree-Fock energy decomposition into diatomic and monoatomic contributions [Chem. Phys. Lett. 382, 265 (2003)], and the Ziegler-Rauk dissociation energy decomposition [Inorg. Chem. 18, 1558 (1979)]. The Ziegler-Rauk scheme is based on a separation of a molecule into fragments, while Mayer's scheme can be used in the cases where a fragmentation of the system in clearly separable parts is not possible. In the Mayer scheme, the density of a free atom is deformed to give the one-atom Mulliken density that subsequently interacts to give rise to the diatomic interaction energy. We give a detailed analysis of the diatomic energy contributions in the Mayer scheme and a close look onto the one-atom Mulliken densities. The Mulliken density ρA has a single large maximum around the nuclear position of the atom A, but exhibits slightly negative values in the vicinity of neighboring atoms. The main connecting point between both analysis schemes is the electrostatic energy. Both decomposition schemes utilize the same electrostatic energy expression, but differ in how fragment densities are defined. In the Mayer scheme, the electrostatic component originates from the interaction of the Mulliken densities, while in the Ziegler-Rauk scheme, the undisturbed fragment densities interact. The values of the electrostatic energy resulting from the two schemes differ significantly but typically have the same order of magnitude. Both methods are useful and complementary since Mayer's decomposition focuses on the energy of the finally formed molecule, whereas the Ziegler-Rauk scheme describes the bond formation starting from undeformed fragment densities

Understanding the two-component sensing System of virulent bacteria

Report for the scientific sojourn carried out at the l’ Institute for Computational Molecular Science of the Temple University, United States, from 2010 to 2012. Two-component systems (TCS) are used by pathogenic bacteria to sense the environment within a host and activate mechanisms related to virulence and antimicrobial resistance. A prototypical example is the PhoQ/PhoP system, which is the major regulator of virulence in Salmonella. Hence, PhoQ is an attractive target for the design of new antibiotics against foodborne diseases. Inhibition of the PhoQ-mediated bacterial virulence does not result in growth inhibition, presenting less selective pressure for the generation of antibiotic resistance. Moreover, PhoQ is a histidine kinase (HK) and it is absent in animals. Nevertheless, the design of satisfactory HK inhibitors has been proven to be a challenge. To compete with the intracellular ATP concentrations, the affinity of a HK inhibidor must be in the micromolar-nanomolar range, whereas the current lead compounds have at best millimolar affinities. Moreover, the drug selectivity depends on the conformation of a highly variable loop, referred to as the “ATP-lid, which is difficult to study by X-Ray crystallography due to its flexibility. I have investigated the binding of different HK inhibitors to PhoQ. In particular, all-atom molecular dynamics simulations have been combined with enhanced sampling techniques in order to provide structural and dynamic information of the conformation of the ATP-lid. Transient interactions between these drugs and the ATP-lid have been identified and the free energy of the different binding modes has been estimated. The results obtained pinpoint the importance of protein flexibility in the HK-inhibitor binding, and constitute a first step in developing more potent and selective drugs. The computational resources of the hosting institution as well as the experience of the members of the group in drug binding and free energy methods have been crucial to carry out this work.