Science.gov

Sample records for large basis ab

  1. Large basis ab initio shell model investigation of {sup 9}Be and {sup 11}Be

    SciTech Connect

    Forssen, C.; Navratil, P.; Ormand, W.E.; Caurier, E.

    2005-04-01

    We present the first ab initio structure investigation of the loosely bound {sup 11}Be nucleus, together with a study of the lighter isotope {sup 9}Be. The nuclear structure of these isotopes is particularly interesting because of the appearance of a parity-inverted ground state in {sup 11}Be. Our study is performed in the framework of the ab initio no-core shell model. Results obtained using four different, high-precision two-nucleon interactions, in model spaces up to 9({Dirac_h}/2{pi}){omega}, are shown. For both nuclei, and all potentials, we reach convergence in the level ordering of positive- and negative-parity spectra separately. Concerning their relative position, the positive-parity states are always too high in excitation energy, but a fast drop with respect to the negative-parity spectrum is observed when the model space is increased. This behavior is most dramatic for {sup 11}Be. In the largest model space we were able to reach, the 1/2{sup +} level has dropped down to become either the first or the second excited state, depending on which interaction we use. We also observe a contrasting behavior in the convergence patterns for different two-nucleon potentials and argue that a three-nucleon interaction is needed to explain the parity inversion. Furthermore, large-basis calculations of {sup 13}C and {sup 11}B are performed. This allows us to study the systematics of the position of the first unnatural-parity state in the N=7 isotone and the A=11 isobar. The {sup 11}B run in the 9({Dirac_h}/2{pi}){omega} model space involves a matrix with dimension exceeding 1.1x10{sup 9}, and is our largest calculation so far. We present results on binding energies, excitation spectra, level configurations, radii, electromagnetic observables, and {sup 10}Be+n overlap functions.

  2. Large-scale ab initio configuration interaction calculations for light nuclei

    NASA Astrophysics Data System (ADS)

    Maris, Pieter; Metin Aktulga, H.; Caprio, Mark A.; Çatalyürek, Ümit V.; Ng, Esmond G.; Oryspayev, Dossay; Potter, Hugh; Saule, Erik; Sosonkina, Masha; Vary, James P.; Yang, Chao; Zhou, Zheng

    2012-12-01

    In ab-initio Configuration Interaction calculations, the nuclear wavefunction is expanded in Slater determinants of single-nucleon wavefunctions and the many-body Schrodinger equation becomes a large sparse matrix problem. The challenge is to reach numerical convergence to within quantified numerical uncertainties for physical observables using finite truncations of the infinite-dimensional basis space. We discuss strategies for constructing and solving the resulting large sparse matrix eigenvalue problems on current multicore computer architectures. Several of these strategies have been implemented in the code MFDn, a hybrid MPI/OpenMP Fortran code for ab-initio nuclear structure calculations that can scale to 100,000 cores and more. Finally, we will conclude with some recent results for 12C including emerging collective phenomena such as rotational band structures using SRG evolved chiral N3LO interactions.

  3. Ab initio no core configuration interaction calculations in the natural orbital basis

    NASA Astrophysics Data System (ADS)

    Constantinou, Chrysovalantis; Caprio, Mark A.; Vary, James P.; Maris, Pieter

    2015-10-01

    The natural orbital basis has been successfully used in the past in atomic and molecular structure calculations. The natural orbitals used in those calculations are calculated by diagonalizing the electron one-body density matrix. Here we develop natural orbitals for nuclear no-core configuration interaction (NCCI) calculations. A NCCI calculation using an initial single particle basis, such as the harmonic oscillator basis, must first be performed in order to obtain a one-body density matrix. The eigenvectors of the one-body density matrix are the natural orbitals, and the corresponding eigenvalues are the occupations of these natural orbitals in the nuclear wave function. According to these occupancies, the most important natural orbitals, in the sense of the most occupied, can then be selected and used in a NCCI calculation. We discuss ab initio nuclear NCCI calculations for light nuclei and assess their ability to provide faster convergence. Supported by the US DOE (under Grants DE-FG02-95ER-40934, DESC0008485 SciDAC/NUCLEI, and DE-FG02-87ER40371), and the US NSF (under Grant 0904782). Computational resources provided by NERSC (supported by US DOE Contract DE-AC02-05CH11231), and NDCRC.

  4. Highly localized quasiatomic minimal basis orbitals for Mo from ab initio calculations

    NASA Astrophysics Data System (ADS)

    Chan, T.-L.; Yao, Y. X.; Wang, C. Z.; Lu, W. C.; Li, J.; Qian, X. F.; Yip, S.; Ho, K. M.

    2007-11-01

    A minimal basis set of localized quasiatomic orbitals for Mo is constructed using the fully converged eigenstates from first-principles calculations with a large basis set. The orbitals, although similar in shape to those of a free atom, are slightly deformed such that it can reproduce all the occupied-state electronic properties of the system. They are very useful for analyzing chemical bonding by calculating the Mulliken overlap population and bond order index between atoms. In addition, the transferability of tight-binding parametrizations can be evaluated, for example, the effect of the two-center approximation.

  5. On the feasibility of ab initio electronic structure calculations for Cu using a single s orbital basis

    SciTech Connect

    Hegde, Ganesh Bowen, R. Chris

    2015-10-15

    The accuracy of a single s-orbital representation of Cu towards enabling multi-thousand atom ab initio calculations of electronic structure is evaluated in this work. If an electrostatic compensation charge of 0.3 electron per atom is used in this basis representation, the electronic transmission in bulk and nanocrystalline Cu can be made to compare accurately to that obtained with a Double Zeta Polarized basis set. The use of this representation is analogous to the use of single band effective mass representation for semiconductor electronic structure. With a basis of just one s-orbital per Cu atom, the representation is extremely computationally efficient and can be used to provide much needed ab initio insight into electronic transport in nanocrystalline Cu interconnects at realistic dimensions of several thousand atoms.

  6. High Yield Chemical Vapor Deposition Growth of High Quality Large-Area AB Stacked Bilayer Graphene

    PubMed Central

    Liu, Lixin; Zhou, Hailong; Cheng, Rui; Yu, Woo Jong; Liu, Yuan; Chen, Yu; Shaw, Jonathan; Zhong, Xing; Huang, Yu; Duan, Xiangfeng

    2012-01-01

    Bernal stacked (AB stacked) bilayer graphene is of significant interest for functional electronic and photonic devices due to the feasibility to continuously tune its band gap with a vertical electrical field. Mechanical exfoliation can be used to produce AB stacked bilayer graphene flakes but typically with the sizes limited to a few micrometers. Chemical vapor deposition (CVD) has been recently explored for the synthesis of bilayer graphene but usually with limited coverage and a mixture of AB and randomly stacked structures. Herein we report a rational approach to produce large-area high quality AB stacked bilayer graphene. We show that the self-limiting effect of graphene growth on Cu foil can be broken by using a high H2/CH4 ratio in a low pressure CVD process to enable the continued growth of bilayer graphene. A high temperature and low pressure nucleation step is found to be critical for the formation of bilayer graphene nuclei with high AB stacking ratio. A rational design of a two-step CVD process is developed for the growth of bilayer graphene with high AB stacking ratio (up to 90 %) and high coverage (up to 99 %). The electrical transport studies demonstrated that devices made of the as-grown bilayer graphene exhibit typical characteristics of AB stacked bilayer graphene with the highest carrier mobility exceeding 4,000 cm2/V·s at room temperature, comparable to that of the exfoliated bilayer graphene. PMID:22906199

  7. Ab Initio Calculation of Nuclear Magnetic Resonance Chemical Shift Anisotropy Tensors 1. Influence of Basis Set on the Calculation of 31P Chemical Shifts

    SciTech Connect

    Alam, T.M.

    1998-09-01

    The influence of changes in the contracted Gaussian basis set used for ab initio calculations of nuclear magnetic resonance (NMR) phosphorous chemical shift anisotropy (CSA) tensors was investigated. The isotropic chemical shitl and chemical shift anisotropy were found to converge with increasing complexity of the basis set at the Hartree-Fock @IF) level. The addition of d polarization function on the phosphorous nucIei was found to have a major impact of the calculated chemical shi~ but diminished with increasing number of polarization fimctions. At least 2 d polarization fimctions are required for accurate calculations of the isotropic phosphorous chemical shift. The introduction of density fictional theory (DFT) techniques through tie use of hybrid B3LYP methods for the calculation of the phosphorous chemical shift tensor resulted in a poorer estimation of the NMR values, even though DFT techniques result in improved energy and force constant calculations. The convergence of the W parametem with increasing basis set complexity was also observed for the DFT calculations, but produced results with consistent large deviations from experiment. The use of a HF 6-31 l++G(242p) basis set represents a good compromise between accuracy of the simulation and the complexity of the calculation for future ab initio calculations of 31P NMR parameters in larger complexes.

  8. Auxiliary basis expansions for large-scale electronic structure calculations

    SciTech Connect

    Jung, Yousung; Sodt, Alexander; Gill, Peter W.M.; Head-Gordon, Martin

    2005-04-04

    One way to reduce the computational cost of electronic structure calculations is to employ auxiliary basis expansions to approximate 4 center integrals in terms of 2 and 3-center integrals, usually using the variationally optimum Coulomb metric to determine the expansion coefficients. However the long-range decay behavior of the auxiliary basis expansion coefficients has not been characterized. We find that this decay can be surprisingly slow. Numerical experiments on linear alkanes and a toy model both show that the decay can be as slow as 1/r in the distance between the auxiliary function and the fitted charge distribution. The Coulomb metric fitting equations also involve divergent matrix elements for extended systems treated with periodic boundary conditions. An attenuated Coulomb metric that is short-range can eliminate these oddities without substantially degrading calculated relative energies. The sparsity of the fit coefficients is assessed on simple hydrocarbon molecules, and shows quite early onset of linear growth in the number of significant coefficients with system size using the attenuated Coulomb metric. This means it is possible to design linear scaling auxiliary basis methods without additional approximations to treat large systems.

  9. An extended basis set {ital ab} {ital initio} study of alkali metal cation--water clusters

    SciTech Connect

    Feller, D.; Glendening, E.D.; Woon, D.E.; Feyereisen, M.W.

    1995-09-01

    Ionic clusters comprised of a single alkali metal cation and up to eight water molecules were studied at the Hartree--Fock and correlated levels of theory using the correlation consistent sequence of basis sets. Estimates of the degree of convergence in the computed properties with respect to the complete basis set limit were facilitated by the underlying systematic manner in which the correlation consistent sets approach completeness. In favorable cases, improved property values could be obtained by fitting finite basis set results with a simple analytical expression in order to extrapolate to the complete basis set limit. The sensitivity of structures and binding energies were analyzed with regard to the inclusion of valence and core-valence correlation recovery at the MP2, MP4, and CCSD(T) levels of theory. The replacement of metal core electrons and the introduction of relativistic contributions via effective core potentials was compared to corresponding all-electron results. {copyright} {ital 1995} {ital American} {ital Institute} {ital of} {ital Physics}.

  10. Ab initio molecular dynamics study of water at constant pressure using converged basis sets and empirical dispersion corrections

    NASA Astrophysics Data System (ADS)

    Ma, Zhonghua; Zhang, Yanli; Tuckerman, Mark E.

    2012-07-01

    It is generally believed that studies of liquid water using the generalized gradient approximation to density functional theory require dispersion corrections in order to obtain reasonably accurate structural and dynamical properties. Here, we report on an ab initio molecular dynamics study of water in the isothermal-isobaric ensemble using a converged discrete variable representation basis set and an empirical dispersion correction due to Grimme [J. Comp. Chem. 27, 1787 (2006)], 10.1002/jcc.20495. At 300 K and an applied pressure of 1 bar, the density obtained without dispersion corrections is approximately 0.92 g/cm3 while that obtained with dispersion corrections is 1.07 g/cm3, indicating that the empirical dispersion correction overestimates the density by almost as much as it is underestimated without the correction for this converged basis. Radial distribution functions exhibit a loss of structure in the second solvation shell. Comparison of our results with other studies using the same empirical correction suggests the cause of the discrepancy: the Grimme dispersion correction is parameterized for use with a particular basis set; this parameterization is sensitive to this choice and, therefore, is not transferable to other basis sets.

  11. Ab initio modeling of carbohydrates: on the proper selection of computational methods and basis sets

    Technology Transfer Automated Retrieval System (TEKTRAN)

    With the development of faster computer hardware and quantum mechanical software it has become more feasible to study large carbohydrate molecules via quantum mechanical methods. In the past, studies of carbohydrates were restricted to empirical/semiempirical methods and Hartree Fock. In the last ...

  12. Ab initio nuclear many-body perturbation calculations in the Hartree-Fock basis

    NASA Astrophysics Data System (ADS)

    Hu, B. S.; Xu, F. R.; Sun, Z. H.; Vary, J. P.; Li, T.

    2016-07-01

    Starting from realistic nuclear forces, the chiral N3LO and JISP16, we have applied many-body perturbation theory (MBPT) to the structure of closed-shell nuclei, 4He and 16O. The two-body N3LO interaction is softened by a similarity renormalization group transformation while JISP16 is adopted without renormalization. The MBPT calculations are performed within the Hartree-Fock (HF) bases. The angular momentum coupled scheme is used, which can reduce the computational task. Corrections up to the third order in energy and up to the second order in radius are evaluated. Higher-order corrections in the HF basis are small relative to the leading-order perturbative result. Using the antisymmetrized Goldstone diagram expansions of the wave function, we directly correct the one-body density for the calculation of the radius, rather than calculate corrections to the occupation probabilities of single-particle orbits as found in other treatments. We compare our results with other methods where available and find good agreement. This supports the conclusion that our methods produce reasonably converged results with these interactions. We also compare our results with experimental data.

  13. Large scale production of Bacillus thuringiensis PS149B1 insecticidal proteins Cry34Ab1 and Cry35Ab1 from Pseudomonas fluorescens.

    PubMed

    Huang, Ke-Xue; Badger, Monty; Haney, Keith; Evans, Steve L

    2007-06-01

    The 14kDa (Cry34Ab1) and 44kDa (Cry35Ab1) binary insecticidal proteins are produced naturally by Bacillus thuringiensis PS149B1 as parasporal inclusion bodies. Here, we show production of these two insecticidal proteins in recombinant Pseudomonas fluorescens and their subsequent purification to near homogeneity to provide large quantities of protein for safety-assessment studies associated with the registration of transgenic corn plants. The gene sequence specific for each protein was expressed in P. fluorescens and fermented at the 75-L scale. For Cry34Ab1, the protein accumulated as insoluble inclusion bodies, and was purified by extraction directly from the cell pastes at pH 3.4 with a sodium acetate buffer, selective precipitation at pH 7.0, and differential centrifugation. For Cry35Ab1, the protein was extracted from the purified inclusion bodies with sodium acetate buffer (pH 3.5) containing 0.5M urea, followed by diafiltration. No chromatography steps were required to produce over 30g of lyophilized protein powder with purity greater than 98%, while retaining full insecticidal activity against Western corn rootworm larvae. The proteins were further characterized to assure identity and suitability for use in safety-assessment studies. PMID:17337206

  14. Recent progress with large-scale ab initio calculations: the CONQUEST code

    NASA Astrophysics Data System (ADS)

    Bowler, D. R.; Choudhury, R.; Gillan, M. J.; Miyazaki, T.

    While the success of density functional theory (DFT) has led to its use in a wide variety of fields such as physics, chemistry, materials science and biochemistry, it has long been recognised that conventional methods are very inefficient for large complex systems, because the memory requirements scale as N 2 and the cpu requirements as N 3 (where N is the number of atoms). The principles necessary to develop methods with linear scaling of the cpu and memory requirements with system size (O(N ) methods) have been established for more than ten years, but only recently have practical codes showing this scaling for DFT started to appear. We report recent progress in the development of the Conquest code, which performs O(N ) DFT calculations on parallel computers, and has a demonstrated ability to handle systems of over 10000 atoms. The code can be run at different levels of precision, ranging from empirical tight-binding, through ab initio tight-binding, to full ab initio , and techniques for calculating ionic forces in a consistent way at all levels of precision will be presented. Illustrations are given of practical Conquest calculations in the strained Ge/Si(001) system.

  15. High-yield chemical vapor deposition growth of high-quality large-area AB-stacked bilayer graphene.

    PubMed

    Liu, Lixin; Zhou, Hailong; Cheng, Rui; Yu, Woo Jong; Liu, Yuan; Chen, Yu; Shaw, Jonathan; Zhong, Xing; Huang, Yu; Duan, Xiangfeng

    2012-09-25

    Bernal-stacked (AB-stacked) bilayer graphene is of significant interest for functional electronic and photonic devices due to the feasibility to continuously tune its band gap with a vertical electric field. Mechanical exfoliation can be used to produce AB-stacked bilayer graphene flakes but typically with the sizes limited to a few micrometers. Chemical vapor deposition (CVD) has been recently explored for the synthesis of bilayer graphene but usually with limited coverage and a mixture of AB- and randomly stacked structures. Herein we report a rational approach to produce large-area high-quality AB-stacked bilayer graphene. We show that the self-limiting effect of graphene growth on Cu foil can be broken by using a high H(2)/CH(4) ratio in a low-pressure CVD process to enable the continued growth of bilayer graphene. A high-temperature and low-pressure nucleation step is found to be critical for the formation of bilayer graphene nuclei with high AB stacking ratio. A rational design of a two-step CVD process is developed for the growth of bilayer graphene with high AB stacking ratio (up to 90%) and high coverage (up to 99%). The electrical transport studies demonstrate that devices made of the as-grown bilayer graphene exhibit typical characteristics of AB-stacked bilayer graphene with the highest carrier mobility exceeding 4000 cm(2)/V · s at room temperature, comparable to that of the exfoliated bilayer graphene. PMID:22906199

  16. Introduction of a fully relativistic capable basis set in the ab initio orthogonalized linear combination of atomic orbitals method

    NASA Astrophysics Data System (ADS)

    Thomas, Patrick Ryan

    Large simulation cell sizes, relativistic effects, and the need to correctly model excited state properties are major impediments to the accurate prediction of the optical properties of candidate materials for solid-state laser crystal and luminescent applications. To overcome these challenges, new methods must be created to improve the electron orbital wavefunction and interactions. In this work, a method has been developed to create new analytical four-component, fully-relativistic and single-component scalar relativistic descriptions of the atomic orbital wave functions from Grasp2K numerically represented atomic orbitals. In addition, adapted theory for the calculation of the relativistic kinetic energy contribution to Hamiltonian which bypasses directly solving the Dirac equation has been explicated. The orbital description improvements are tested against YAG, YBCO, SnO2 and BiF3. The improvements to the basis set reflect an improvement in both computational speed and accuracy.

  17. Systematic chemical variations in large 3AB iron meteorites: Clues to core crystallization

    NASA Technical Reports Server (NTRS)

    Haack, H.; Scott, E. R. D.; Rubio, G. S.; Gutierrez, D. F.; Lewis, C. F.; Wasson, J. T.; Brooks, R. R.; Guo, X.; Ryan, D. E.; Holzbecher, J.

    1993-01-01

    Analysis of numerous individual iron meteorites have shown that fractional crystallization of iron cores result in variations in chemical concentration of the solid core which span several orders of magnitude. The magnitude and direction of the resulting spatial gradients in the core can provide clues to the physical nature of the core crystallization process. We have analyzed suites of samples from three large 3AB irons (Cape York, 58t; Chupaderos, 24t; Morito, 10t) in order to estimate local chemical gradients. Initial results for the concentrations of Ge, Pd, Pt (Massey group), Ir, Au, As, Co, Os, and Rh (Dalhouse group), and P (Arizona group) show significant ranges among the Cape York and Chupaderos samples and marginally significant ranges among the Morito samples. Measurements of Au, Ir, Co, Ni, Cu, Ga, As, W, Re (from UCLA) and Ni and Co (Arizona group) are in progress. We find a spatial Ir gradient in Chupaderos with a magnitude similar to the one reported for Agpalilik (Cape York iron) by Esbensen et al.

  18. Investigation of new phases in the Ba-Si phase diagram under high pressure using ab initio structural search.

    PubMed

    Shi, Jingming; Cui, Wenwen; Flores-Livas, José A; San-Miguel, Alfonso; Botti, Silvana; Marques, Miguel A L

    2016-03-01

    Barium silicides are versatile materials that have attracted attention for a variety of applications in electronics and optoelectronics. Using an unbiased structural search based on a particle-swarm optimization algorithm combined with density functional theory calculations, we investigate systematically the ground-state phase stability and the structural diversity of Ba-Si binaries under high pressure. The phase diagram turns out to be quite intricate, with several compositions stabilizing/destabilizing as a function of pressure. In particular, we identify novel phases of BaSi, BaSi2, BaSi3, and BaSi5 that might be synthesizable experimentally over a wide range of pressures. Our results not only clarify and complete the previously known structural phase diagram, but also provide new insights for understanding the Ba-Si binary system. PMID:26923068

  19. The Genetic Basis for Altered Blood Vessel Function in Disease: Large Artery Stiffening

    PubMed Central

    Agrotis, Alex

    2005-01-01

    The progressive stiffening of the large arteries in humans that occurs during aging constitutes a potential risk factor for increased cardiovascular morbidity and mortality, and is accompanied by an elevation in systolic blood pressure and pulse pressure. While the underlying basis for these changes remains to be fully elucidated, factors that are able to influence the structure and composition of the extracellular matrix and the way it interacts with arterial smooth muscle cells could profoundly affect the properties of the large arteries. Thus, while age and sex represent important factors contributing to large artery stiffening, the variation in growth-stimulating factors and those that modulate extracellular production and homeostasis are also being increasingly recognized to play a key role in the process. Therefore, elucidating the contribution that genetic variation makes to large artery stiffening could ultimately provide the basis for clinical strategies designed to regulate the process for therapeutic benefit. PMID:17315605

  20. On basis set superposition error corrected stabilization energies for large n-body clusters.

    PubMed

    Walczak, Katarzyna; Friedrich, Joachim; Dolg, Michael

    2011-10-01

    In this contribution, we propose an approximate basis set superposition error (BSSE) correction scheme for the site-site function counterpoise and for the Valiron-Mayer function counterpoise correction of second order to account for the basis set superposition error in clusters with a large number of subunits. The accuracy of the proposed scheme has been investigated for a water cluster series at the CCSD(T), CCSD, MP2, and self-consistent field levels of theory using Dunning's correlation consistent basis sets. The BSSE corrected stabilization energies for a series of water clusters are presented. A study regarding the possible savings with respect to computational resources has been carried out as well as a monitoring of the basis set dependence of the approximate BSSE corrections. PMID:21992293

  1. Effect of noncovalent interactions on conformers of the n-butylbenzene monomer studied by mass analyzed threshold ionization spectroscopy and basis-set convergent ab initio computations.

    PubMed

    Tong, Xin; Cerný, Jirí; Müller-Dethlefs, Klaus; Dessent, Caroline E H

    2008-07-01

    Two conformational isomers of the aromatic hydrocarbon n-butylbenzene have been studied using two-color MATI (mass analyzed threshold ionization) spectroscopy to explore the effect of conformation on ionization dynamics. Cationic states of g auche-conformer III and anti- conformers IV were selectively produced by two-color excitation via the respective S 1 origins. Adiabatic ionization potentials of the gauche- and anti-conformations were determined to be 70146 and 69872 +/- 5 cm (-1) respectively. Spectral features and vibrational modes are interpreted with the aid of MP2/cc-pVDZ ab initio calculations, and ionization-induced changes in the molecular conformations are discussed. Complete basis set (CBS) ab initio studies at MP2 level reveal reliable energetics for all four n-butylbenzene conformers observed in earlier two-color REMPI (resonance enhanced multiphoton ionization) spectra. For the S 0 state, the energies of conformer III, IV and V are above conformer I by 130, 289, 73 cm (-1), respectively. Furthermore, the combination of the CBS calculations with the measured REMPI, MATI spectra allowed the determination of the energetics of all four conformers in the S 1 and D 0 states. PMID:18533642

  2. Large scale ab initio calculations of extended defects in materials: screw dislocations in bcc metals

    NASA Astrophysics Data System (ADS)

    Dézerald, Lucile; Ventelon, Lisa; Willaime, François; Clouet, Emmanuel; Rodney, David

    2014-06-01

    Ab initio methods, based on the Density Functional Theory (DFT), have been extensively used to study point defects and defect clusters in materials. Present HPC resources and DFT codes now allow similar investigations to be performed on dislocations. The study of these extended defects requires not only larger simulation cells but also a higher accuracy because the energy differences, which are involved, are rather small, typically 50-to-100 meV for supercells containing 50-to-500 atoms. The topology of the Peierls potential of screw dislocations with 1/2 <111>Burgers vector, i.e. the 2D energy landscape seen by these dislocations, is being completely revisited by DFT calculations. From results obtained in all body-centered cubic (bcc) transition metals, except Cr (V, Nb, Ta, Mo, W and Fe), using the PWSCF code, which is part of the Quantum-Espresso package, we concluded that the 2D Peierls potentials have two common features: the single-hump shape of the barrier between two minima of the potential, and the presence of a maximum - and not a minimum as predicted by most empirical potentials - around the split core. In iron, the topology of the Peierls potential is reversed compared to the classical sinusoidal picture: the location of the saddle point and the maximum are indeed inverted with unexpected flat regions. The first results obtained within the framework of the PRACE project, DIMAIM (DIslocations in Metals using Ab Initio Methods), started at the beginning of 2013, will also be presented. In particular, in order to address the twinning-antitwinning asymmetry often observed in bcc metals and regarded as the major contribution to the breakdown of Schmid's law, we have determined the crystal orientation dependence of the Peierls stress, i.e. the critical stress required for dislocation motion. These computationally most expensive simulations were performed on the PRACE Tier-0 system at Barcelona Supercomputing Center (Marenostrum III). The scalability results

  3. Use of reduced basis technique in the inverse dynamics of large space cranes

    NASA Technical Reports Server (NTRS)

    Das, S. K.; Utku, S.; Wada, B. K.

    1990-01-01

    The inverse dynamics of adaptive structures used as space cranes can prove computationally expensive in the case of large structures, due to the large number of degrees of freedom involved. Consequently, reduced basis techniques (reduction techniques) are frequently used to reduce the problem size to a time manageable level (for possible use in real time control). A reduced basis technique is proposed which is different from, but related to, the path-derivatives reduction technique. A linearly independent set of deflection n-tuples is used, chosen at the beginning of the time range in which it is wished to reduce the equations, in whose subspace it is assumed that the deflection vectors of the unreduced problem will lie (approximately).

  4. Accurate ab initio-based adiabatic global potential energy surface for the 22A″ state of NH2 by extrapolation to the complete basis set limit

    NASA Astrophysics Data System (ADS)

    Li, Y. Q.; Ma, F. C.; Sun, M. T.

    2013-10-01

    A full three-dimensional global potential energy surface is reported first time for the title system, which is important for the photodissociation processes. It is obtained using double many-body expansion theory and an extensive set of accurate ab initio energies extrapolated to the complete basis set limit. Such a work can be recommended for dynamics studies of the N(2D) + H2 reaction, a reliable theoretical treatment of the photodissociation dynamics and as building blocks for constructing the double many-body expansion potential energy surface of larger nitrogen/hydrogen containing systems. In turn, a preliminary theoretical study of the reaction N(^2D)+H_2(X^1Σ _g^+)(ν =0,j=0)rArr NH(a^1Δ )+H(^2S) has been carried out with the method of quasi-classical trajectory on the new potential energy surface. Integral cross sections and thermal rate constants have been calculated, providing perhaps the most reliable estimate of the integral cross sections and the rate constants known thus far for such a reaction.

  5. Accurate ab initio-based adiabatic global potential energy surface for the 2(2)A" state of NH2 by extrapolation to the complete basis set limit.

    PubMed

    Li, Y Q; Ma, F C; Sun, M T

    2013-10-21

    A full three-dimensional global potential energy surface is reported first time for the title system, which is important for the photodissociation processes. It is obtained using double many-body expansion theory and an extensive set of accurate ab initio energies extrapolated to the complete basis set limit. Such a work can be recommended for dynamics studies of the N((2)D) + H2 reaction, a reliable theoretical treatment of the photodissociation dynamics and as building blocks for constructing the double many-body expansion potential energy surface of larger nitrogen/hydrogen containing systems. In turn, a preliminary theoretical study of the reaction N((2)D)+H2(X(1)Σg (+))(ν=0,j=0)→NH(a(1)Δ)+H((2)S) has been carried out with the method of quasi-classical trajectory on the new potential energy surface. Integral cross sections and thermal rate constants have been calculated, providing perhaps the most reliable estimate of the integral cross sections and the rate constants known thus far for such a reaction. PMID:24160511

  6. Efficient solution of Poisson's equation using discrete variable representation basis sets for Car-Parrinello ab initio molecular dynamics simulations with cluster boundary conditions.

    PubMed

    Lee, Hee-Seung; Tuckerman, Mark E

    2008-12-14

    An efficient computational approach to perform Car-Parrinello ab initio molecular dynamics (CPAIMD) simulations under cluster (free) boundary conditions is presented. The general approach builds upon a recent real-space CPAIMD formalism using discrete variable representation (DVR) basis sets [Y. Liu et al., Phys. Rev. B 12, 125110 (2003); H.-S. Lee and M. E. Tuckerman, J. Phys. Chem. A 110, 5549 (2006)]. In order to satisfy cluster boundary conditions, a DVR based on sinc functions is utilized to expand the Kohn-Sham orbitals and electron density. Poisson's equation is solved in order to calculate the Hartree potential via an integral representation of the 1/r singularity. Excellent convergence properties are achieved with respect to the number of grid points (or DVR functions) and the size of the simulation cell. A straightforward implementation of the present approach leads to near linear scaling [O(N(4/3))] of the computational cost with respect to the system size (N) for the solution of Poisson's equation. The accuracy and stability of CPAIMD simulations based on sinc DVR are tested for a model problem as well as for N(2) and a water dimer. PMID:19071908

  7. Ghost transmission: How large basis sets can make electron transport calculations worse

    SciTech Connect

    Herrmann, Carmen; Solomon, Gemma C.; Subotnik, Joseph E.; Mujica, Vladimiro; Ratner, Mark A.

    2010-01-01

    The Landauer approach has proven to be an invaluable tool for calculating the electron transport properties of single molecules, especially when combined with a nonequilibrium Green’s function approach and Kohn–Sham density functional theory. However, when using large nonorthogonal atom-centered basis sets, such as those common in quantum chemistry, one can find erroneous results if the Landauer approach is applied blindly. In fact, basis sets of triple-zeta quality or higher sometimes result in an artificially high transmission and possibly even qualitatively wrong conclusions regarding chemical trends. In these cases, transport persists when molecular atoms are replaced by basis functions alone (“ghost atoms”). The occurrence of such ghost transmission is correlated with low-energy virtual molecular orbitals of the central subsystem and may be interpreted as a biased and thus inaccurate description of vacuum transmission. An approximate practical correction scheme is to calculate the ghost transmission and subtract it from the full transmission. As a further consequence of this study, it is recommended that sensitive molecules be used for parameter studies, in particular those whose transmission functions show antiresonance features such as benzene-based systems connected to the electrodes in meta positions and other low-conducting systems such as alkanes and silanes.

  8. Efficient preconditioning of the electronic structure problem in large scale ab initio molecular dynamics simulations.

    PubMed

    Schiffmann, Florian; VandeVondele, Joost

    2015-06-28

    We present an improved preconditioning scheme for electronic structure calculations based on the orbital transformation method. First, a preconditioner is developed which includes information from the full Kohn-Sham matrix but avoids computationally demanding diagonalisation steps in its construction. This reduces the computational cost of its construction, eliminating a bottleneck in large scale simulations, while maintaining rapid convergence. In addition, a modified form of Hotelling's iterative inversion is introduced to replace the exact inversion of the preconditioner matrix. This method is highly effective during molecular dynamics (MD), as the solution obtained in earlier MD steps is a suitable initial guess. Filtering small elements during sparse matrix multiplication leads to linear scaling inversion, while retaining robustness, already for relatively small systems. For system sizes ranging from a few hundred to a few thousand atoms, which are typical for many practical applications, the improvements to the algorithm lead to a 2-5 fold speedup per MD step. PMID:26133420

  9. Efficient preconditioning of the electronic structure problem in large scale ab initio molecular dynamics simulations

    NASA Astrophysics Data System (ADS)

    Schiffmann, Florian; VandeVondele, Joost

    2015-06-01

    We present an improved preconditioning scheme for electronic structure calculations based on the orbital transformation method. First, a preconditioner is developed which includes information from the full Kohn-Sham matrix but avoids computationally demanding diagonalisation steps in its construction. This reduces the computational cost of its construction, eliminating a bottleneck in large scale simulations, while maintaining rapid convergence. In addition, a modified form of Hotelling's iterative inversion is introduced to replace the exact inversion of the preconditioner matrix. This method is highly effective during molecular dynamics (MD), as the solution obtained in earlier MD steps is a suitable initial guess. Filtering small elements during sparse matrix multiplication leads to linear scaling inversion, while retaining robustness, already for relatively small systems. For system sizes ranging from a few hundred to a few thousand atoms, which are typical for many practical applications, the improvements to the algorithm lead to a 2-5 fold speedup per MD step.

  10. Efficient preconditioning of the electronic structure problem in large scale ab initio molecular dynamics simulations

    SciTech Connect

    Schiffmann, Florian; VandeVondele, Joost

    2015-06-28

    We present an improved preconditioning scheme for electronic structure calculations based on the orbital transformation method. First, a preconditioner is developed which includes information from the full Kohn-Sham matrix but avoids computationally demanding diagonalisation steps in its construction. This reduces the computational cost of its construction, eliminating a bottleneck in large scale simulations, while maintaining rapid convergence. In addition, a modified form of Hotelling’s iterative inversion is introduced to replace the exact inversion of the preconditioner matrix. This method is highly effective during molecular dynamics (MD), as the solution obtained in earlier MD steps is a suitable initial guess. Filtering small elements during sparse matrix multiplication leads to linear scaling inversion, while retaining robustness, already for relatively small systems. For system sizes ranging from a few hundred to a few thousand atoms, which are typical for many practical applications, the improvements to the algorithm lead to a 2-5 fold speedup per MD step.

  11. Large-scale ab initio simulations of binary transition metal clusters for storage media materials.

    PubMed

    Entel, P; Gruner, M E

    2009-02-11

    In the quest for ultra-high-density magnetic recording, new materials in the nanometre range have attracted much interest over the last decade involving intense studies of L1(0) phases of contemporary or future storage media materials like FePt or CoPt nanoparticles. Based on large-scale density functional theory calculations, we provide a systematic overview of the structural and magnetic properties of various morphologies of FePt and CoPt nanoclusters with diameters up to 3 nm. In this size range, the ordered multiply twinned morphologies are energetically favoured over the nanoclusters with the desired layer type L1(0) and high magnetocrystalline anisotropy. Other nanoparticles of interest, like FePd, also show a preference for multiply twinned structures or exhibit, as in the case of MnPt nanoclusters, a strong tendency for antiferromagnetic ordering instead of ferromagnetic order. The compositional trends of the various nanoparticles can be traced back to differences in the partial electronic density of states of the 3d element. PMID:21715930

  12. Application of the MELCOR code to design basis PWR large dry containment analysis.

    SciTech Connect

    Phillips, Jesse; Notafrancesco, Allen; Tills, Jack Lee

    2009-05-01

    The MELCOR computer code has been developed by Sandia National Laboratories under USNRC sponsorship to provide capability for independently auditing analyses submitted by reactor manufactures and utilities. MELCOR is a fully integrated code (encompassing the reactor coolant system and the containment building) that models the progression of postulated accidents in light water reactor power plants. To assess the adequacy of containment thermal-hydraulic modeling incorporated in the MELCOR code for application to PWR large dry containments, several selected demonstration designs were analyzed. This report documents MELCOR code demonstration calculations performed for postulated design basis accident (DBA) analysis (LOCA and MSLB) inside containment, which are compared to other code results. The key processes when analyzing the containment loads inside PWR large dry containments are (1) expansion and transport of high mass/energy releases, (2) heat and mass transfer to structural passive heat sinks, and (3) containment pressure reduction due to engineered safety features. A code-to-code benchmarking for DBA events showed that MELCOR predictions of maximum containment loads were equivalent to similar predictions using a qualified containment code known as CONTAIN. This equivalency was found to apply for both single- and multi-cell containment models.

  13. Structural basis for the superior activity of the large isoform of snow flea antifreeze protein.

    PubMed

    Mok, Yee-Foong; Lin, Feng-Hsu; Graham, Laurie A; Celik, Yeliz; Braslavsky, Ido; Davies, Peter L

    2010-03-23

    The snow flea (Hypogastrum harveyi) is protected from freezing at sub-zero temperatures by a glycine-rich antifreeze protein (AFP) that binds to seed ice crystals and prevents them from growing larger. This AFP is hyperactive and comprises two isoforms [Graham, L. A., and Davies, P. L. (2005) Science 310, 461]. The larger isoform (15.7 kDa) exhibits several-fold higher activity than the smaller isoform (6.5 kDa), although it is considerably less abundant. To establish the molecular basis for this difference in activity, we determined the sequence of the large isoform. The primary sequences of these two isoforms are surprisingly divergent. However, both contain tripeptide repeats and turn motifs that enabled us to build a three-dimensional model of the large isoform based upon the six-polyproline helix structure of the small isoform. Our model contains 13 polyproline type II helices connected by proline-containing loops stacked into two flat sheets oriented antiparallel to one another. The structure is strictly amphipathic, with a hydrophilic surface on one side and a hydrophobic, putative ice-binding surface on the other. The putative ice-binding site is approximately twice as large in area as that of the small isoform, providing an explanation for the difference in activity that is consistent with other examples noted. By tagging the recombinant AFP with green fluorescent protein, we observed its binding to multiple planes of ice, especially the basal plane. This finding supports the correlation between AFP hyperactivity and basal plane binding first observed with spruce budworm AFP. PMID:20158269

  14. Structural response of large penetrations and closures for containment vessels subjected to loadings beyond design basis

    SciTech Connect

    Kulak, R.F.; Hsieh, B.J.; Ash, J.E.; Kennedy, J.M.; McLennan, G.A.; Pan, Y.C.

    1985-02-01

    The Reactor Analysis and Safety Division (RAS) and the Components Technology Division (CT) of Argonne National Laboratory (ANL) are performing analytical/numerical simulations of the response of selected large penetrations and closures, which use some type of seal or gasketed joint, for containment vessels subject to pressure and thermal loads that are beyond the design basis (BDB). The objectives of this task were to identify the methodology required to simulate the structural response of selected penetrations/closures to BDB loadings and to apply this methodology to representative penetrations/closures. Section II discusses a detailed study conducted to determine the structural response of an equipment hatch for a PWR with a steel containment vessel. The macro-deformations of the gasketed joint were computed and then used in a leakage analysis. In Section III, the methodology used to assess the structural integrity of a BWR Mark II containment head is presented. Section IV describes the approach used to obtain upper and lower bounds for the maximum allowable internal pressure and deflection for a representative bellows connection.

  15. The molecular basis for ANE syndrome revealed by the large ribosomal subunit processome interactome

    PubMed Central

    McCann, Kathleen L; Teramoto, Takamasa; Zhang, Jun; Tanaka Hall, Traci M; Baserga, Susan J

    2016-01-01

    ANE syndrome is a ribosomopathy caused by a mutation in an RNA recognition motif of RBM28, a nucleolar protein conserved to yeast (Nop4). While patients with ANE syndrome have fewer mature ribosomes, it is unclear how this mutation disrupts ribosome assembly. Here we use yeast as a model system and show that the mutation confers growth and pre-rRNA processing defects. Recently, we found that Nop4 is a hub protein in the nucleolar large subunit (LSU) processome interactome. Here we demonstrate that the ANE syndrome mutation disrupts Nop4’s hub function by abrogating several of Nop4’s protein-protein interactions. Circular dichroism and NMR demonstrate that the ANE syndrome mutation in RRM3 of human RBM28 disrupts domain folding. We conclude that the ANE syndrome mutation generates defective protein folding which abrogates protein-protein interactions and causes faulty pre-LSU rRNA processing, thus revealing one aspect of the molecular basis of this human disease. DOI: http://dx.doi.org/10.7554/eLife.16381.001 PMID:27077951

  16. Phase 1 trial of IL-15 trans presentation blockade using humanized Mik-Beta-1 mAb in patients with T-cell large granular lymphocytic leukemia

    PubMed Central

    Conlon, Kevin C.; Stewart, Donn M.; Worthy, TatYana A.; Janik, John E.; Fleisher, Thomas A.; Albert, Paul S.; Figg, William D.; Spencer, Shawn D.; Raffeld, Mark; Decker, Jean R.; Goldman, Carolyn K.; Bryant, Bonita R.; Petrus, Michael N.; Creekmore, Stephen P.; Morris, John C.

    2013-01-01

    In the present study, Hu-Mikβ1, a humanized mAb directed at the shared IL-2/IL-15Rβ subunit (CD122) was evaluated in patients with T-cell large granular lymphocytic (T-LGL) leukemia. Hu-Mikβ1 blocked the trans presentation of IL-15 to T cells expressing IL-2/IL-15Rβ and the common γ-chain (CD132), but did not block IL-15 action in cells that expressed the heterotrimeric IL-15 receptor in cis. There was no significant toxicity associated with Hu-Mikβ1 administration in patients with T-LGL leukemia, but no major clinical responses were observed. One patient who had previously received murine Mikβ1 developed a measurable Ab response to the infused Ab. Nevertheless, the safety profile of this first in-human study of the humanized mAb to IL-2/IL-15Rβ (CD122) supports its evaluation in disorders such as refractory celiac disease, in which IL-15 and its receptor have been proposed to play a critical role in the pathogenesis and maintenance of disease activity. The protocol is registered with www.clinicaltrials.gov as number NCT 00076180. PMID:23212516

  17. Complete Genome Sequence of the Multiresistant Acinetobacter baumannii Strain AbH12O-A2, Isolated during a Large Outbreak in Spain

    PubMed Central

    Merino, M.; Alvarez-Fraga, L.; Gómez, M. J.; Aransay, A. M.; Lavín, J. L.; Chaves, F.

    2014-01-01

    We report the complete genome sequence of Acinetobacter baumannii strain AbH12O-A2, isolated during a large outbreak in Spain. The genome has 3,875,775 bp and 3,526 coding sequences, with 39.4% G+C content. The availability of this genome will facilitate the study of the pathogenicity of the Acinetobacter species. PMID:25395646

  18. Complete Genome Sequence of the Multiresistant Acinetobacter baumannii Strain AbH12O-A2, Isolated during a Large Outbreak in Spain.

    PubMed

    Merino, M; Alvarez-Fraga, L; Gómez, M J; Aransay, A M; Lavín, J L; Chaves, F; Bou, G; Poza, M

    2014-01-01

    We report the complete genome sequence of Acinetobacter baumannii strain AbH12O-A2, isolated during a large outbreak in Spain. The genome has 3,875,775 bp and 3,526 coding sequences, with 39.4% G+C content. The availability of this genome will facilitate the study of the pathogenicity of the Acinetobacter species. PMID:25395646

  19. Molecular basis for the inhibition of β-hydroxyacyl-ACP dehydratase HadAB complex from Mycobacterium tuberculosis by flavonoid inhibitors.

    PubMed

    Dong, Yu; Qiu, Xiaodi; Shaw, Neil; Xu, Yueyang; Sun, Yuna; Li, Xuemei; Li, Jun; Rao, Zihe

    2015-07-01

    Dehydration is one of the key steps in the biosynthesis of mycolic acids and is vital to the growth of Mycobacterium tuberculosis (Mtb). Consequently, stalling dehydration cures tuberculosis (TB). Clinically used anti-TB drugs like thiacetazone (TAC) and isoxyl (ISO) as well as flavonoids inhibit the enzyme activity of the β-hydroxyacyl-ACP dehydratase HadAB complex. How this inhibition is exerted, has remained an enigma for years. Here, we describe the first crystal structures of the MtbHadAB complex bound with flavonoid inhibitor butein, 2',4,4'-trihydroxychalcone or fisetin. Despite sharing no sequence identity from Blast, HadA and HadB adopt a very similar hotdog fold. HadA forms a tight dimer with HadB in which the proteins are sitting side-by-side, but are oriented anti-parallel. While HadB contributes the catalytically critical His-Asp dyad, HadA binds the fatty acid substrate in a long channel. The atypical double hotdog fold with a single active site formed by MtbHadAB gives rise to a long, narrow cavity that vertically traverses the fatty acid binding channel. At the base of this cavity lies Cys61, which upon mutation to Ser confers drug-resistance in TB patients. We show that inhibitors bind in this cavity and protrude into the substrate binding channel. Thus, inhibitors of MtbHadAB exert their effect by occluding substrate from the active site. The unveiling of this mechanism of inhibition paves the way for accelerating development of next generation of anti-TB drugs. PMID:26081470

  20. Solving the Dirac equation, using the large component only, in a Dirac-type Slater orbital basis set

    NASA Astrophysics Data System (ADS)

    van Lenthe, E.; Baerends, E. J.; Snijders, J. G.

    1995-04-01

    We solve the Dirac equation by solving the two-component energy-dependent equation for the large component that results from the elimination of the small component. This requires for every occupied orbital the diagonalization of a Hamiltonian. Advantages are, however, that these Hamiltonians are all bounded from below, unlike the Dirac Hamiltonian, and that only a basis set for the large component is needed. We use Dirac-type Slater orbitals, adapted from solutions to the hydrogen-like atom. This offers the perspective of performing relativistic calculations to the same accuracy as non-relativistic ones, with a comparable number of basis functions.

  1. On the performance of large Gaussian basis sets for the computation of total atomization energies

    NASA Technical Reports Server (NTRS)

    Martin, J. M. L.

    1992-01-01

    The total atomization energies of a number of molecules have been computed using an augmented coupled-cluster method and (5s4p3d2f1g) and 4s3p2d1f) atomic natural orbital (ANO) basis sets, as well as the correlation consistent valence triple zeta plus polarization (cc-pVTZ) correlation consistent valence quadrupole zeta plus polarization (cc-pVQZ) basis sets. The performance of ANO and correlation consistent basis sets is comparable throughout, although the latter can result in significant CPU time savings. Whereas the inclusion of g functions has significant effects on the computed Sigma D(e) values, chemical accuracy is still not reached for molecules involving multiple bonds. A Gaussian-1 (G) type correction lowers the error, but not much beyond the accuracy of the G1 model itself. Using separate corrections for sigma bonds, pi bonds, and valence pairs brings down the mean absolute error to less than 1 kcal/mol for the spdf basis sets, and about 0.5 kcal/mol for the spdfg basis sets. Some conclusions on the success of the Gaussian-1 and Gaussian-2 models are drawn.

  2. Toward accurate thermochemical models for transition metals : G3large basis sets for atoms Sc-Zn.

    SciTech Connect

    Mayhall, N. J.; Raghavachari, K.; Redfern, P. C.; Curtiss, L. A.; Rassolov, V.; Indiana Univ.; Univ. of South Carolina

    2008-04-01

    An augmented valence triple-zeta basis set, referred to as G3Large, is reported for the first-row transition metal elements Sc through Zn. The basis set is constructed in a manner similar to the G3Large basis set developed previously for other elements (H-Ar, K, Ca, Ga-Kr) and used as a key component in Gaussian-3 theory. It is based on a contraction of a set of 15s13p5d Gaussian primitives to 8s7p3d, and also includes sets of f and g polarization functions, diffuse spd functions, and core df polarization functions. The basis set is evaluated with triples-augmented coupled cluster [CCSD(T)] and Brueckner orbital [BD(T)] methods for a small test set involving energies of atoms, atomic ions, and diatomic hydrides. It performs well for the low-lying s{yields}d excitation energies of atoms, atomic ionization energies, and the dissociation energies of the diatomic hydrides. The Brueckner orbital-based BD(T) method performs substantially better than Hartree-Fock-based CCSD(T) for molecules such as NiH, where the starting unrestricted Hartree-Fock wavefunction suffers from a high degree of spin contamination. Comparison with available data for geometries of transition metal hydrides also shows good agreement. A smaller basis set without core polarization functions, G3MP2Large, is also defined.

  3. TECHNICAL BASIS REPORT FOR LARGE FIRE ACCIDENTS INVOLVING ABOVEGROUND TANKS & VESSELS

    SciTech Connect

    MARCHESE, A.R.

    2005-03-03

    This document analyzes large fire accidents involving aboveground tanks and vessels during Demonstration Bulk Vitrification System (DBVS) operations. The fire accident scenarios are consistent with RPP-22461,''Preliminary Fire Hazard Analysis (PFHA) for DBVS''. The radiological and toxicological consequences are determined for a wide spectrum of fire sizes to bracket the range of possible consequences resulting from large fires involving aboveground tanks/vessels that are part of DBVS.

  4. Possible calcium centers for hydrogen storage applications: An accurate many-body study by AFQMC calculations with large basis sets

    NASA Astrophysics Data System (ADS)

    Purwanto, Wirawan; Krakauer, Henry; Zhang, Shiwei; Virgus, Yudistira

    2011-03-01

    Weak H2 physisorption energies present a significant challenge to first-principle theoretical modeling and prediction of materials for H storage. There has been controversy regarding the accuracy of DFT on systems involving Ca cations. We use the auxiliary-field quantum Monte Carlo (AFQMC) method to accurately predict the binding energy of Ca + , - 4{H}2 . AFQMC scales as Nbasis3and has demonstrated accuracy similar to or better than the gold-standard coupled cluster CCSD(T) method. We apply a modified Cholesky decomposition to achieve efficient Hubbard-Stratonovich transformation in AFQMC at large basis sizes. We employ the largest correlation consistent basis sets available, up to Ca/cc-pCV5Z, to extrapolate to the complete basis limit. The calculated potential energy curve exhibits binding with a double-well structure. Supported by DOE and NSF. Calculations were performed at OLCF Jaguar and CPD.

  5. Defect-Enhanced Charge Transfer by Ion-Solid Interactions in SiC using Large-Scale Ab Initio Molecular Dynamics Simulations

    SciTech Connect

    Gao, Fei; Xiao, H. Y.; Zu, Xiaotao T.; Posselt, Matthias; Weber, William J.

    2009-07-10

    Large-scale ab initio molecular dynamics simulations of ion-solid interactions in SiC reveal that significant charge-transfer occurs between atoms and defects can enhance charge transfer to surrounding atoms. The results demonstrate that charge transfer to and from recoiling atoms can alter the energy barriers and dynamics for stable defect formation. The present simulations illustrate in detail the dynamic processes for charged defect formation. The averaged values of displacement threshold energies along four main crystallographic directions are smaller than those determined by empirical potentials due to charge transfer effects on recoil atoms.

  6. Defect-Enhanced Charge Transfer by Ion-Solid Interactions in SiC using Large-Scale Ab Initio Molecular Dynamics Simulations

    SciTech Connect

    Gao Fei; Weber, William J.; Xiao Haiyan; Zu Xiaotao; Posselt, Matthias

    2009-07-10

    Large-scale ab initio molecular dynamics simulations of ion-solid interactions in SiC reveal that significant charge transfer occurs between atoms, and defects can enhance charge transfer to surrounding atoms. The results demonstrate that charge transfer to and from recoiling atoms can alter the energy barriers and dynamics for stable defect formation. The present simulations illustrate in detail the dynamic processes for charged defect formation. The averaged values of displacement threshold energies along four main crystallographic directions are smaller than those determined by empirical potentials due to charge-transfer effects on recoil atoms.

  7. Large Tanker Motion Model Identification Using Generalized Ellipsoidal Basis Function-Based Fuzzy Neural Networks.

    PubMed

    Wang, Ning; Er, Meng Joo; Han, Min

    2015-12-01

    In this paper, the motion dynamics of a large tanker is modeled by the generalized ellipsoidal function-based fuzzy neural network (GEBF-FNN). The reference model of tanker motion dynamics in the form of nonlinear difference equations is established to generate training data samples for the GEBF-FNN algorithm which begins with no hidden neuron. In the sequel, fuzzy rules associated with the GEBF-FNN-based model can be online self-constructed by generation criteria and parameter estimation, and can dynamically capture essential motion dynamics of the large tanker with high prediction accuracy. Simulation studies and comprehensive comparisons are conducted on typical zig-zag maneuvers with moderate and extreme steering, and demonstrate that the GEBF-FNN-based model of tanker motion dynamics achieves superior performance in terms of both approximation and prediction. PMID:25561605

  8. Structural basis for translational surveillance by the large ribosomal subunit-associated protein quality control complex

    PubMed Central

    Lyumkis, Dmitry; Oliveira dos Passos, Dario; Tahara, Erich B.; Webb, Kristofor; Bennett, Eric J.; Vinterbo, Staal; Potter, Clinton S.; Carragher, Bridget; Joazeiro, Claudio A. P.

    2014-01-01

    All organisms have evolved mechanisms to manage the stalling of ribosomes upon translation of aberrant mRNA. In eukaryotes, the large ribosomal subunit-associated quality control complex (RQC), composed of the listerin/Ltn1 E3 ubiquitin ligase and cofactors, mediates the ubiquitylation and extraction of ribosome-stalled nascent polypeptide chains for proteasomal degradation. How RQC recognizes stalled ribosomes and performs its functions has not been understood. Using single-particle cryoelectron microscopy, we have determined the structure of the RQC complex bound to stalled 60S ribosomal subunits. The structure establishes how Ltn1 associates with the large ribosomal subunit and properly positions its E3-catalytic RING domain to mediate nascent chain ubiquitylation. The structure also reveals that a distinguishing feature of stalled 60S particles is an exposed, nascent chain-conjugated tRNA, and that the Tae2 subunit of RQC, which facilitates Ltn1 binding, is responsible for selective recognition of stalled 60S subunits. RQC components are engaged in interactions across a large span of the 60S subunit surface, connecting the tRNA in the peptidyl transferase center to the distally located nascent chain tunnel exit. This work provides insights into a mechanism linking translation and protein degradation that targets defective proteins immediately after synthesis, while ignoring nascent chains in normally translating ribosomes. PMID:25349383

  9. Large variation among photoreceptors as the basis of visual flexibility in the common backswimmer

    PubMed Central

    Immonen, Esa-Ville; Ignatova, Irina; Gislen, Anna; Warrant, Eric; Vähäsöyrinki, Mikko; Weckström, Matti; Frolov, Roman

    2014-01-01

    The common backswimmer, Notonecta glauca, uses vision by day and night for functions such as underwater prey animal capture and flight in search of new habitats. Although previous studies have identified some of the physiological mechanisms facilitating such flexibility in the animal's vision, neither the biophysics of Notonecta photoreceptors nor possible cellular adaptations are known. Here, we studied Notonecta photoreceptors using patch-clamp and intracellular recording methods. Photoreceptor size (approximated by capacitance) was positively correlated with absolute sensitivity and acceptance angles. Information rate measurements indicated that large and more sensitive photoreceptors performed better than small ones. Our results suggest that backswimmers are adapted for vision in both dim and well-illuminated environments by having open-rhabdom eyes with large intrinsic variation in absolute sensitivity among photoreceptors, exceeding those found in purely diurnal or nocturnal species. Both electrophysiology and microscopic analysis of retinal structure suggest two retinal subsystems: the largest peripheral photoreceptors provide vision in dim light and the smaller peripheral and central photoreceptors function primarily in sunlight, with light-dependent pigment screening further contributing to adaptation in this system by dynamically recruiting photoreceptors with varying sensitivity into the operational pool. PMID:25274359

  10. Large variation among photoreceptors as the basis of visual flexibility in the common backswimmer.

    PubMed

    Immonen, Esa-Ville; Ignatova, Irina; Gislen, Anna; Warrant, Eric; Vähäsöyrinki, Mikko; Weckström, Matti; Frolov, Roman

    2014-11-22

    The common backswimmer, Notonecta glauca, uses vision by day and night for functions such as underwater prey animal capture and flight in search of new habitats. Although previous studies have identified some of the physiological mechanisms facilitating such flexibility in the animal's vision, neither the biophysics of Notonecta photoreceptors nor possible cellular adaptations are known. Here, we studied Notonecta photoreceptors using patch-clamp and intracellular recording methods. Photoreceptor size (approximated by capacitance) was positively correlated with absolute sensitivity and acceptance angles. Information rate measurements indicated that large and more sensitive photoreceptors performed better than small ones. Our results suggest that backswimmers are adapted for vision in both dim and well-illuminated environments by having open-rhabdom eyes with large intrinsic variation in absolute sensitivity among photoreceptors, exceeding those found in purely diurnal or nocturnal species. Both electrophysiology and microscopic analysis of retinal structure suggest two retinal subsystems: the largest peripheral photoreceptors provide vision in dim light and the smaller peripheral and central photoreceptors function primarily in sunlight, with light-dependent pigment screening further contributing to adaptation in this system by dynamically recruiting photoreceptors with varying sensitivity into the operational pool. PMID:25274359

  11. Large Earthquakes at the Ibero-Maghrebian Region: Basis for an EEWS

    NASA Astrophysics Data System (ADS)

    Buforn, Elisa; Udías, Agustín; Pro, Carmen

    2015-09-01

    Large earthquakes (Mw > 6, Imax > VIII) occur at the Ibero-Maghrebian region, extending from a point (12ºW) southwest of Cape St. Vincent to Tunisia, with different characteristics depending on their location, which cause considerable damage and casualties. Seismic activity at this region is associated with the boundary between the lithospheric plates of Eurasia and Africa, which extends from the Azores Islands to Tunisia. The boundary at Cape St. Vincent, which has a clear oceanic nature in the westernmost part, experiences a transition from an oceanic to a continental boundary, with the interaction of the southern border of the Iberian Peninsula, the northern border of Africa, and the Alboran basin between them, corresponding to a wide area of deformation. Further to the east, the plate boundary recovers its oceanic nature following the northern coast of Algeria and Tunisia. The region has been divided into four zones with different seismic characteristics. From west to east, large earthquake occurrence, focal depth, total seismic moment tensor, and average seismic slip velocities for each zone along the region show the differences in seismic release of deformation. This must be taken into account in developing an EEWS for the region.

  12. Optically driven nanostructures as the basis for large-scale quantum computing

    NASA Astrophysics Data System (ADS)

    Tsukanov, Alexander V.

    2008-03-01

    We propose a large-scale quantum computer architecture based upon the regular arrays of dopant atoms implanted into the semiconductor host matrix. The singly-ionized pairs of donors represent charge qubits on which arbitrary quantum operations can be achieved by application of two strongly detuned laser pulses. The implementation of two-qubit operations as well as the qubit read-out utilize the intermediate circuit containing a probe electron that is able to shuttle along the array of ionized ancilla donors providing the indirect conditional coupling between the qubits. The quantum bus strategy enables us to handle the qubits connected in parallel and enhances the efficiency of the quantum information processing. We demonstrate that non-trivial multi-qubit operations in the quantum register (e.g., an entanglement generation) can be accomplished by the sequence of the optical pulses combined with an appropriate voltage gate pattern.

  13. Accurate ab initio double many-body expansion potential energy surface for ground-state H2S by extrapolation to the complete basis set limit.

    PubMed

    Song, Y Z; Varandas, A J C

    2009-04-01

    A single-sheeted potential energy surface is reported for the electronic ground-state of H(2)S by fitting accurate multireference configuration interaction energies calculated using aug-cc-pVTZ and aug-cc-pVQZ basis sets with extrapolation of the electron correlation energy to the complete basis set limit, plus extrapolation to the complete basis set limit of the complete-active-space self-consistent field energy. A switching function formalism has been used to warrant the correct behavior at the H(2)(X (1)Sigma(g) (+))+S((1)D) and SH(X (2)Pi)+H((2)S) dissociation limits. The topographical features of the novel global potential energy surface are examined in detail, with the former being used for exploratory quasiclassical trajectory calculations of the thermal rate constant for the S((1)D)+H(2), S((1)D)+D(2), and S((1)D)+HD reactions at room temperature. A comparison with other available potential energy surfaces as well as kinetics data is also provided. PMID:19355742

  14. Automatic location of L/H transition times for physical studies with a large statistical basis

    NASA Astrophysics Data System (ADS)

    González, S.; Vega, J.; Murari, A.; Pereira, A.; Dormido-Canto, S.; Ramírez, J. M.; contributors, JET-EFDA

    2012-06-01

    Completely automatic techniques to estimate and validate L/H transition times can be essential in L/H transition analyses. The generation of databases with hundreds of transition times and without human intervention is an important step to accomplish (a) L/H transition physics analysis, (b) validation of L/H theoretical models and (c) creation of L/H scaling laws. An entirely unattended methodology is presented in this paper to build large databases of transition times in JET using time series. The proposed technique has been applied to a dataset of 551 JET discharges between campaigns C21 and C26. A prediction with discharges that show a clear signature in time series is made through the locating properties of the wavelet transform. It is an accurate prediction and the uncertainty interval is ±3.2 ms. The discharges with a non-clear pattern in the time series use an L/H mode classifier based on discharges with a clear signature. In this case, the estimation error shows a distribution with mean and standard deviation of 27.9 ms and 37.62 ms, respectively. Two different regression methods have been applied to the measurements acquired at the transition times identified by the automatic system. The obtained scaling laws for the threshold power are not significantly different from those obtained using the data at the transition times determined manually by the experts. The automatic methods allow performing physical studies with a large number of discharges, showing, for example, that there are statistically different types of transitions characterized by different scaling laws.

  15. Structural basis of ligand interactions of the large extracellular domain of tetraspanin CD81.

    PubMed

    Rajesh, Sundaresan; Sridhar, Pooja; Tews, Birke Andrea; Fénéant, Lucie; Cocquerel, Laurence; Ward, Douglas G; Berditchevski, Fedor; Overduin, Michael

    2012-09-01

    Hepatitis C virus (HCV) causes chronic liver disease, cirrhosis, and primary liver cancer. Despite 130 million people being at risk worldwide, no vaccine exists, and effective therapy is limited by drug resistance, toxicity, and high costs. The tetraspanin CD81 is an essential entry-level receptor required for HCV infection of hepatocytes and represents a critical target for intervention. In this study, we report the first structural characterization of the large extracellular loop of CD81, expressed in mammalian cells and studied in physiological solutions. The HCV E2 glycoprotein recognizes CD81 through a dynamic loop on the helical bundle, which was shown by nuclear magnetic resonance (NMR) spectroscopy to adopt a conformation distinct from that seen in crystals. A novel membrane binding interface was revealed adjacent to the exposed HCV interaction site in the extracellular loop of CD81. The binding pockets for two proposed inhibitors of the CD81-HCV interaction, namely, benzyl salicylate and fexofenadine, were shown to overlap the HCV and membrane interaction sites. Although the dynamic loop region targeted by these compounds presents challenges for structure-based design, the NMR assignments enable realistic screening and validation of ligands. Together, these data provide an improved avenue for developing potent agents that specifically block CD81-HCV interaction and also pave a way for elucidating the recognition mechanisms of diverse tetraspanins. PMID:22740401

  16. Overcoming the existent computational challenges in the ab initio calculations of the two-photon circular dichroism spectra of large molecules using a fragment-recombination approach

    NASA Astrophysics Data System (ADS)

    Diaz, Carlos; Echevarria, Lorenzo; Hernández, Florencio E.

    2013-05-01

    Herein we report on the development of a fragment-recombination approach (FRA) that allows overcoming the computational limitations found in the ab initio calculation of the two-photon circular dichroism (TPCD) spectra of large optically active molecules. Through the comparative analysis of the corresponding theoretical TPCD spectra of the fragments and that of the entire molecule, we prove that TPCD is an additive property. We also demonstrate that the same property apply to two-photon absorption (TPA). TPCD-FRA is expected to find great applications in the structural-analysis of large catalysts and polypeptides due to its reduced computational complexity, cost and time, and to reveal fingerprints in the obscure spectral region between the near and far UV.

  17. Accurate adiabatic potential energy surface for 12A' state of FH2 based on ab initio data extrapolated to the complete basis set limit

    NASA Astrophysics Data System (ADS)

    Li, Yong-Qing; Song, Yu-Zhi; Joaquim de Campos Varandas, António

    2015-01-01

    An accurate single-sheeted double many-body expansion potential energy surface is reported for the title system. It is obtained by using the aug-cc-pVTZ and aug-cc-pVQZ basis sets with extrapolation of the electron correlation energy to the complete basis set limit, plus extrapolation to the complete basis set limit of the complete-active-space self-consistent field energy. The collinear and bending barrier heights of the new global potential energy surface is 2.301 and 1.768 kcal mol-1, in very good agreement with the values of 2.222 and 1.770 kcal mol-1 from the current best potential energy surface. In particular, the new potential energy surface describes well the important van der Waals interactions which is very useful for investigating the dynamics of the title system. Thus, the new potential energy surface can both be recommended for dynamics studies of the F + H2 reaction and as building block for constructing the potential energy surfaces of larger fluorine/hydrogen containing systems. Based on the new potential energy surface, a preliminary theoretical study of the reaction F(2P) + H2 (X1 Σg+) → FH(X1Σ+) + H(2S) has been carried out with the methods of quasi-classical trajectory and quantum mechanical. The results have shown that the new PES is suitable for any kind of dynamics studies. Supplementary material in the form of one pdf file available from the Journal web page at http://dx.doi.org/10.1140/epjd/e2014-50445-3

  18. Euglena light-harvesting chlorophyll A/B binding protein (LHCP) synthesized as an unusually large precursor

    SciTech Connect

    Rikin, A.; Meyer, A.; Schwartzbach, S.

    1987-04-01

    Light increased the rate of LHCP synthesis as measured by pulse-labeling with /sup 35/SO/sub 4/ and immunoprecipitation with antibody specific for Euglena LHCP. In addition to the mature LHCP, 26,000 daltons, the LHCP specific antibody immunoprecipitated large amounts of several proteins having molecular weights of approximately 100,000. On immunoblots of immunoprecipitated unlabeled protein, the antibody only detected the mature LHCP suggesting that the high molecular weight proteins are not LHCP aggregates produced during immunoprecipitation. After a 10 min pulse with /sup 35/SO/sub 4/, the 100,000 dalton proteins constituted over 80% of the immunoprecipitated material. In a subsequent chase, the radioactivity in the 100,000 dalton proteins decreased and the radioactivity in the mature LHCP increased suggesting a precursor-product relationship. After a 35 minute chase, the mature LHCP was the major radioactive protein immunoprecipitated. Peptide mapping and in vitro translation are being used to clarify the structural and functional relationships, if any, between the 100,000 and 26,000 dalton immunoprecipitation products.

  19. Accurate ab initio-based adiabatic global potential energy surface for the 2{sup 2}A″ state of NH{sub 2} by extrapolation to the complete basis set limit

    SciTech Connect

    Li, Y. Q.; Ma, F. C.; Sun, M. T.

    2013-10-21

    A full three-dimensional global potential energy surface is reported first time for the title system, which is important for the photodissociation processes. It is obtained using double many-body expansion theory and an extensive set of accurate ab initio energies extrapolated to the complete basis set limit. Such a work can be recommended for dynamics studies of the N({sup 2}D) + H{sub 2} reaction, a reliable theoretical treatment of the photodissociation dynamics and as building blocks for constructing the double many-body expansion potential energy surface of larger nitrogen/hydrogen containing systems. In turn, a preliminary theoretical study of the reaction N({sup 2}D)+H{sub 2}(X{sup 1}Σ{sub g}{sup +})(ν=0,j=0)→NH(a{sup 1}Δ)+H({sup 2}S) has been carried out with the method of quasi-classical trajectory on the new potential energy surface. Integral cross sections and thermal rate constants have been calculated, providing perhaps the most reliable estimate of the integral cross sections and the rate constants known thus far for such a reaction.

  20. Remembering AB

    NASA Astrophysics Data System (ADS)

    Belyayev, S. T.

    2013-06-01

    In 1947 I became a second-year student at Moscow State University's Physics and Engineering Department, where a part of the week's classes were taught at base organizations. Our group's base was the future Kurchatov Institute, at that time known as the mysterious "Laboratory N^circ 2," and later as LIPAN. . Besides group lectures and practical work at the experimental laboratories, we also had access to the general seminars which Igor Vasilyevich Kurchatov tried to hold, with Leonid Vasilyevich Groshev filling in when he was absent. At the seminar, theorists spoke as welcome co-presenters and commentators. In 1949 I felt ready to approach A. B. Migdal to ask if I could transfer to his theoretical sector. In response, he suggested a number of simple qualitative problems, which I then successfully solved. (Incidentally, AB used the very same "introductory problems" for screening many generations of students.) So I wound up among AB's students. From 1952 on (for 10 years) I also served as an employee of the Migdal Sector. My memoirs here are mainly inspired by these years of constant communication with AB. After my departure for Novosibirsk in 1962, although our meetings still took place, they became occasional....

  1. Evaluation of the Genetic Basis of Familial Aggregation of Pacemaker Implantation by a Large Next Generation Sequencing Panel

    PubMed Central

    Steiner, Hillel A.; Uradu, Andrea; Lynnes, Ty C.; Groh, William J.; Miller, John M.; Lin, Hai; Gao, Hongyu; Wang, Zhiping; Liu, Yunlong; Chen, Peng-Sheng; Vatta, Matteo

    2015-01-01

    Background The etiology of conduction disturbances necessitating permanent pacemaker (PPM) implantation is often unknown, although familial aggregation of PPM (faPPM) suggests a possible genetic basis. We developed a pan-cardiovascular next generation sequencing (NGS) panel to genetically characterize a selected cohort of faPPM. Materials and Methods We designed and validated a custom NGS panel targeting the coding and splicing regions of 246 genes with involvement in cardiac pathogenicity. We enrolled 112 PPM patients and selected nine (8%) with faPPM to be analyzed by NGS. Results Our NGS panel covers 95% of the intended target with an average of 229x read depth at a minimum of 15-fold depth, reaching a SNP true positive rate of 98%. The faPPM patients presented with isolated cardiac conduction disease (ICCD) or sick sinus syndrome (SSS) without overt structural heart disease or identifiable secondary etiology. Three patients (33.3%) had heterozygous deleterious variants previously reported in autosomal dominant cardiac diseases including CCD: LDB3 (p.D117N) and TRPM4 (p.G844D) variants in patient 4; TRPM4 (p.G844D) and ABCC9 (p.V734I) variants in patient 6; and SCN5A (p.T220I) and APOB (p.R3527Q) variants in patient 7. Conclusion FaPPM occurred in 8% of our PPM clinic population. The employment of massive parallel sequencing for a large selected panel of cardiovascular genes identified a high percentage (33.3%) of the faPPM patients with deleterious variants previously reported in autosomal dominant cardiac diseases, suggesting that genetic variants may play a role in faPPM. PMID:26636822

  2. Authorization basis requirements comparison report

    SciTech Connect

    Brantley, W.M.

    1997-08-18

    The TWRS Authorization Basis (AB) consists of a set of documents identified by TWRS management with the concurrence of DOE-RL. Upon implementation of the TWRS Basis for Interim Operation (BIO) and Technical Safety Requirements (TSRs), the AB list will be revised to include the BIO and TSRs. Some documents that currently form part of the AB will be removed from the list. This SD identifies each - requirement from those documents, and recommends a disposition for each to ensure that necessary requirements are retained when the AB is revised to incorporate the BIO and TSRs. This SD also identifies documents that will remain part of the AB after the BIO and TSRs are implemented. This document does not change the AB, but provides guidance for the preparation of change documentation.

  3. Formation of large-grain-sized BaSi2 epitaxial layers grown on Si(111) by molecular beam epitaxy

    NASA Astrophysics Data System (ADS)

    Baba, M.; Toh, K.; Toko, K.; Hara, K. O.; Usami, N.; Saito, N.; Yoshizawa, N.; Suemasu, T.

    2013-09-01

    BaSi2 epitaxial films were grown on Si(111) substrates by a two-step growth method including reactive deposition epitaxy (RDE) and molecular beam epitaxy (MBE). To enlarge the grain size of BaSi2, the Ba deposition rate and duration were varied from 0.25 to 1.0 nm/min and from 5 to 120 min during RDE, respectively. The effect of post-annealing was also investigated at 760 °C for 10 min. Plan-view transmission electron micrographs indicated that the grain size in the MBE-grown BaSi2 was significantly increased up to approximately 4.0 μm, which is much larger than 0.2 μm, reported previously.

  4. Spin-orbit decomposition of ab initio nuclear wave functions

    NASA Astrophysics Data System (ADS)

    Johnson, Calvin W.

    2015-03-01

    Although the modern shell-model picture of atomic nuclei is built from single-particle orbits with good total angular momentum j , leading to j -j coupling, decades ago phenomenological models suggested that a simpler picture for 0 p -shell nuclides can be realized via coupling of the total spin S and total orbital angular momentum L . I revisit this idea with large-basis, no-core shell-model calculations using modern ab initio two-body interactions and dissect the resulting wave functions into their component L - and S -components. Remarkably, there is broad agreement with calculations using the phenomenological Cohen-Kurath forces, despite a gap of nearly 50 years and six orders of magnitude in basis dimensions. I suggest that L -S decomposition may be a useful tool for analyzing ab initio wave functions of light nuclei, for example, in the case of rotational bands.

  5. Ab initio nuclear structure theory

    NASA Astrophysics Data System (ADS)

    Negoita, Gianina Alina

    Ab initio no core methods have become major tools for understanding the properties of light nuclei based on realistic nucleon-nucleon (NN) and three-nucleon (NNN) interactions. A brief description is provided for the inter-nucleon interactions that fit two-body scattering and bound state data, as well as NNN interactions. Major new progress, including the goal of applying these interactions to solve for properties of nuclei, is limited by convergence issues. That is, with the goal of obtaining high precision solutions of the nuclear many-body Hamiltonian with no core methods (all nucleons treated on the same footing), one needs to proceed to very large basis spaces to achieve a convergence pattern suitable for extrapolation to the exact result. This thesis investigates (1) the similarity renormalization group (SRG) approach to soften the interaction, while preserving its phase shift properties, and (2) adoption of a realistic basis space using Woods-Saxon (WS) single-particle wavefunctions. Both have their advantages and limitations, discussed here. For (1), SRG was demonstrated by applying it to a realistic NN interaction, JISP16, in a harmonic oscillator (HO) representation. The degree of interaction softening achieved through a regulator parameter is examined. For (2), new results are obtained with the realistic JISP16 NN interaction in ab initio calculations of light nuclei 4He, 6He and 12C, using a WS basis optimized to minimize the ground-state energy within the truncated no core shell model. These are numerically-intensive many-body calculations. Finally, to gain insight into the potential for no core investigations of heavier nuclei, an initial investigation was obtained for the odd mass A = 47 - 49 region nuclei straddling 48Ca. The motivation for selecting these nuclei stems from the aim of preparing for nuclear double beta-decay studies of 48Ca. In these heavier systems, phenomenological additions to the realistic NN interaction determined by previous

  6. Approaching the complete basis set limit of CCSD(T) for large systems by the third-order incremental dual-basis set zero-buffer F12 method

    NASA Astrophysics Data System (ADS)

    Zhang, Jun; Dolg, Michael

    2014-01-01

    The third-order incremental dual-basis set zero-buffer approach was combined with CCSD(T)-F12x (x = a, b) theory to develop a new approach, i.e., the inc3-db-B0-CCSD(T)-F12 method, which can be applied as a black-box procedure to efficiently obtain the near complete basis set (CBS) limit of the CCSD(T) energies also for large systems. We tested this method for several cases of different chemical nature: four complexes taken from the standard benchmark sets S66 and X40, the energy difference between isomers of water hexamer and the rotation barrier of biphenyl. The results show that our method has an error relative to the best estimation of CBS energy of only 0.2 kcal/mol or less. By parallelization, our method can accomplish the CCSD(T)-F12 calculations of about 60 correlated electrons and 800 basis functions in only several days, which by standard implementation are impossible for ordinary hardware. We conclude that the inc3-db-B0-CCSD(T)-F12a/AVTZ method, which is of CCSD(T)/AV5Z quality, is close to the limit of accuracy that one can achieve for large systems currently.

  7. Approaching the complete basis set limit of CCSD(T) for large systems by the third-order incremental dual-basis set zero-buffer F12 method

    SciTech Connect

    Zhang, Jun Dolg, Michael

    2014-01-28

    The third-order incremental dual-basis set zero-buffer approach was combined with CCSD(T)-F12x (x = a, b) theory to develop a new approach, i.e., the inc3-db-B0-CCSD(T)-F12 method, which can be applied as a black-box procedure to efficiently obtain the near complete basis set (CBS) limit of the CCSD(T) energies also for large systems. We tested this method for several cases of different chemical nature: four complexes taken from the standard benchmark sets S66 and X40, the energy difference between isomers of water hexamer and the rotation barrier of biphenyl. The results show that our method has an error relative to the best estimation of CBS energy of only 0.2 kcal/mol or less. By parallelization, our method can accomplish the CCSD(T)-F12 calculations of about 60 correlated electrons and 800 basis functions in only several days, which by standard implementation are impossible for ordinary hardware. We conclude that the inc3-db-B0-CCSD(T)-F12a/AVTZ method, which is of CCSD(T)/AV5Z quality, is close to the limit of accuracy that one can achieve for large systems currently.

  8. Changes in the Large-Scale Atmospheric Circulation over Romania Between 1961 and 2010 on Seasonal Basis

    NASA Astrophysics Data System (ADS)

    Barbu, Nicu; Burada, Cristina; Stefan, Sabina; Georgescu, Florinela

    2016-04-01

    The aim of this paper is to investigate the trends and shifts of the circulation types over Romania for 50-year period (1961-2010) on seasonal basis. In order to achieve this, two objective catalogues, namely GWT and WLK, from COST733 Action were employed. Daily circulation types were grouped according to the cyclonicity and anticyclonicity and were used to calculate the seasonal occurrence frequency of cyclonic and anticyclonic types. The trend of seasonal time series was investigated by using Mann-Kendall test and the shifts points were determined by using Pettitt test. The results reveal that the occurrence frequency of anticyclonic types increases in summer and winter seasons and the occurrence frequency of cyclonic ones decreases for the summer season (for alpha = 0.05).

  9. Single particle calculations for a Woods-Saxon potential with triaxial deformations, and large Cartesian oscillator basis

    NASA Astrophysics Data System (ADS)

    Mohammed-Azizi, B.; Medjadi, D. E.

    2004-01-01

    We present a computer program which solves the Schrodinger equation of the stationary states for an average nuclear potential of Woods-Saxon type. In this work, we take specifically into account triaxial (i.e. ellipsoidal) nuclear surfaces. The deformation is specified by the usual Bohr parameters. The calculations are carried out in two stages. In the first, one calculates the representative matrix of the Hamiltonian in the Cartesian oscillator basis. In the second stage one diagonalizes this matrix with the help of subroutines of the Eispack library. If it is wished, one can calculate all eigenvalues, or only the part of the eigenvalues that are contained in a fixed interval defined in advance. In this latter case the eigenvectors are given conjointly. The program is very rapid, and the run-time is mainly used for the diagonalization. Thus, it is possible to use a significant number of the basis states in order to insure a best convergence of the results. Program summaryProgram obtainable from:CPC Program Library, Queen's University of Belfast, N. Ireland Title of program:Triaxial Catalogue number:ADSK Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADSK Licensing provisions:None Computer:PC. AMD Athlon 1000 MHz Hard disk:40 Go Ram:256 Mo Swap file:4 Go Operating system:WINDOWS XP Software used:Microsoft Visual Fortran 5.0A (with full optimizations in the settings project options) Programming language:Fortran 77/90 (double precision) Number of bits in a word:32 Number of lines:7662 No. of bytes in distributed program, including test data, etc.:174 601 Distribution format:tar gzip file Nature of the problem: The single particle energies and the single particle wave functions are calculated from one-body Hamiltonian including a central field of Woods-Saxon type, a spin-orbit interaction, and the Coulomb potential for the protons. We consider only ellipsoidal (triaxial) shapes. The deformation of the nuclear shape is fixed by the usual Bohr parameters ( β,

  10. Hybrid functionals for large periodic systems in an all-electron, numeric atom-centered basis framework

    NASA Astrophysics Data System (ADS)

    Levchenko, Sergey V.; Ren, Xinguo; Wieferink, Jürgen; Johanni, Rainer; Rinke, Patrick; Blum, Volker; Scheffler, Matthias

    2015-07-01

    We describe a framework to evaluate the Hartree-Fock exchange operator for periodic electronic-structure calculations based on general, localized atom-centered basis functions. The functionality is demonstrated by hybrid-functional calculations of properties for several semiconductors. In our implementation of the Fock operator, the Coulomb potential is treated either in reciprocal space or in real space, where the sparsity of the density matrix can be exploited for computational efficiency. Computational aspects, such as the rigorous avoidance of on-the-fly disk storage, and a load-balanced parallel implementation, are also discussed. We demonstrate linear scaling of our implementation with system size by calculating the electronic structure of a bulk semiconductor (GaAs) with up to 1,024 atoms per unit cell without compromising the accuracy.

  11. Single particle calculations for a Woods Saxon potential with triaxial deformations, and large Cartesian oscillator basis (new version code)

    NASA Astrophysics Data System (ADS)

    Mohammed-Azizi, B.; Medjadi, D. E.

    2007-05-01

    We present a new version of the computer program which solves the Schrödinger equation of the stationary states for an average nuclear potential of Woods-Saxon type. In this work, we take specifically into account triaxial (i.e. ellipsoidal) nuclear surfaces. The deformation is specified by the usual Bohr parameters. The calculations are carried out in two stages. In the first, one calculates the representative matrix of the Hamiltonian in the Cartesian oscillator basis. In the second stage one diagonalizes this matrix with the help of subroutines of the EISPACK library. This new version calculates all the eigenvalues up to a given cutoff energy, and gives the components of the corresponding eigenfunctions. For a more convenient handling, these results are stored simultaneously in the computer memory, and on a files. Program summaryTitle of program:Triaxial2007 Catalogue identifier:ADSK_v2_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADSK_v2_0 Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Summary of revision:One input file instead two. Reduced number of input parameters. Storage of eigenvalues and eigenvectors in memory in a very simple way which makes the code very convenient to the user. Reasons for the new version: More convenient handling of the eigenvectors Catalogue number old version: ADSK Catalogue number new version:ADSK_v2_0 Journal: Computer Physics Commun. 156 (2004) 241-282 Licensing provisions: none Computer: PC Pentium 4, 2600 MHz Hard disk: 40 Gb RAM: 256 Mb Swap file: 4 Gb Operating system: WINDOWS XP Software used: Compaq Visual FORTRAN (with full optimizations in the settings project options) Programming language used:Fortran 77/90 (double precision) Number of bits in a word: 32 No. of lines in distributed program, including test data, etc.:4058 No. of bytes in distributed program, including test data, etc.:75 590 Distribution format:tar.gz Nature of the problem: The single particle energies

  12. Large-Scale Atmospheric Circulation Patterns Associated with Temperature Extremes as a Basis for Model Evaluation: Methodological Overview and Results

    NASA Astrophysics Data System (ADS)

    Loikith, P. C.; Broccoli, A. J.; Waliser, D. E.; Lintner, B. R.; Neelin, J. D.

    2015-12-01

    Anomalous large-scale circulation patterns often play a key role in the occurrence of temperature extremes. For example, large-scale circulation can drive horizontal temperature advection or influence local processes that lead to extreme temperatures, such as by inhibiting moderating sea breezes, promoting downslope adiabatic warming, and affecting the development of cloud cover. Additionally, large-scale circulation can influence the shape of temperature distribution tails, with important implications for the magnitude of future changes in extremes. As a result of the prominent role these patterns play in the occurrence and character of extremes, the way in which temperature extremes change in the future will be highly influenced by if and how these patterns change. It is therefore critical to identify and understand the key patterns associated with extremes at local to regional scales in the current climate and to use this foundation as a target for climate model validation. This presentation provides an overview of recent and ongoing work aimed at developing and applying novel approaches to identifying and describing the large-scale circulation patterns associated with temperature extremes in observations and using this foundation to evaluate state-of-the-art global and regional climate models. Emphasis is given to anomalies in sea level pressure and 500 hPa geopotential height over North America using several methods to identify circulation patterns, including self-organizing maps and composite analysis. Overall, evaluation results suggest that models are able to reproduce observed patterns associated with temperature extremes with reasonable fidelity in many cases. Model skill is often highest when and where synoptic-scale processes are the dominant mechanisms for extremes, and lower where sub-grid scale processes (such as those related to topography) are important. Where model skill in reproducing these patterns is high, it can be inferred that extremes are

  13. Discriminating the molecular basis of hepatotoxicity using the large-scale characteristic molecular signatures of toxicants by expression profiling analysis.

    PubMed

    Eun, Jung Woo; Ryu, So Yeon; Noh, Ji Heon; Lee, Min-Jae; Jang, Ja-Jun; Ryu, Jae Chun; Jung, Kwang Hwa; Kim, Jeong Kyu; Bae, Hyun Jin; Xie, Hongjian; Kim, Su Young; Lee, Sug Hyung; Park, Won Sang; Yoo, Nam Jin; Lee, Jung Young; Nam, Suk Woo

    2008-07-30

    Predicting the potential human health risk posed by chemical stressors has long been a major challenge for toxicologists, and the use of microarrays to measure responses to toxicologically relevant genes, and to identify selective, sensitive biomarkers of toxicity is a major application of predictive and discovery toxicology. To investigate this possibility, we investigated whether carcinogens (at doses known to induce liver tumors in chronic exposure bioassays) deregulate characteristic sets of genes in mice. Male C3H/He mice were dosed with two hepatocarcinogens (vinyl chloride (VC, 50-25 mg/kg), aldrin (AD, 0.8-0.4 mg/kg)), or two non-hepatocarcinogens (copper sulfate (CS, 150-60 mg/kg), 2,4,5-trichlorophenoxyacetic acid (2,4,5-T, 150-60 mg/kg)). Large-scale molecular changes elicited by these four hepatotoxicants in liver tissues were analyzed using DNA microarray. Three days after administration, no significant phenotypic changes were induced by these four different hepatotoxicants in terms of histological examination or blood biochemical assay. However, unsupervised hierarchical analysis of gene expressional changes induced by hepatotoxicants resulted in two major gene subclusters on dendrogram, i.e., a carcinogen (VN, AD) and non-carcinogen group (CS, 2,4,5-T), and also revealed that distinct molecular signatures exist. These signatures were founded on well-defined functional gene categories and may differentiate genotoxic and non-genotoxic carcinogens. Furthermore, Venn diagram analysis allowed us to identify carcinogen and non-carcinogen-associated molecular signatures. Using statistical methods, we analyzed outlier genes for four different classes (genotoxic-, non-genotoxic-carcinogen, genotoxic-, non-genotoxic non-carcinogen) in terms of their potential to predict different modes-of-action. In conclusion, the identification of large-scale molecular changes in different hepatocarcinogen exposure models revealed that different types of hepatotoxicants are

  14. Small Atomic Orbital Basis Set First-Principles Quantum Chemical Methods for Large Molecular and Periodic Systems: A Critical Analysis of Error Sources.

    PubMed

    Sure, Rebecca; Brandenburg, Jan Gerit; Grimme, Stefan

    2016-04-01

    In quantum chemical computations the combination of Hartree-Fock or a density functional theory (DFT) approximation with relatively small atomic orbital basis sets of double-zeta quality is still widely used, for example, in the popular B3LYP/6-31G* approach. In this Review, we critically analyze the two main sources of error in such computations, that is, the basis set superposition error on the one hand and the missing London dispersion interactions on the other. We review various strategies to correct those errors and present exemplary calculations on mainly noncovalently bound systems of widely varying size. Energies and geometries of small dimers, large supramolecular complexes, and molecular crystals are covered. We conclude that it is not justified to rely on fortunate error compensation, as the main inconsistencies can be cured by modern correction schemes which clearly outperform the plain mean-field methods. PMID:27308221

  15. Small Atomic Orbital Basis Set First‐Principles Quantum Chemical Methods for Large Molecular and Periodic Systems: A Critical Analysis of Error Sources

    PubMed Central

    Sure, Rebecca; Brandenburg, Jan Gerit

    2015-01-01

    Abstract In quantum chemical computations the combination of Hartree–Fock or a density functional theory (DFT) approximation with relatively small atomic orbital basis sets of double‐zeta quality is still widely used, for example, in the popular B3LYP/6‐31G* approach. In this Review, we critically analyze the two main sources of error in such computations, that is, the basis set superposition error on the one hand and the missing London dispersion interactions on the other. We review various strategies to correct those errors and present exemplary calculations on mainly noncovalently bound systems of widely varying size. Energies and geometries of small dimers, large supramolecular complexes, and molecular crystals are covered. We conclude that it is not justified to rely on fortunate error compensation, as the main inconsistencies can be cured by modern correction schemes which clearly outperform the plain mean‐field methods. PMID:27308221

  16. Physiological basis for large differences in resistance to nitrite among freshwater and freshwater-acclimated euryhaline fishes.

    PubMed

    Tomasso, Joseph R; Grosell, Martin

    2005-01-01

    Uptake of environmental NO2- by most freshwater fishes occurs at the gills where NO2- is actively transported into the blood by the Cl- uptake pathway. Some freshwater fishes do not concentrate NO2- in their plasma, regardless of environmental NO2- exposure and exhibit a high degree of resistance to NO2-. Recent studies indicate that freshwater-adapted killifish (Fundulidae: Fundulus heteroclitus) and European eel (Anguillidae: Anguilla anguilla) have no or minimal Cl- uptake activity at the gills relative to most freshwater fishes; rather, Cl- requirements are met in other ways (probably dietary). We hypothesized that different rates of Cl- uptake by the gill may explain the observed differences in NO2- uptake and consequent toxicity among freshwater fishes. Cl- influx rates of channel catfish (Ictaluridae: Ictalurus punctatus), a species that concentrates NO2- in the plasma and is sensitive to NO2-, and bluegill (Centrarchidae: Lepomis macrochirus), a species that does not concentrate NO2- in the plasma and is resistant to NO2-, were determined over a range of environmental Cl- concentrations. Channel catfish actively transported chloride into the plasma (Km = 155.6+/-101.2 micromol/L Cl-; Jmax = 414.9+/-51.4 nmol/g/h; +/-SEM). In contrast, bluegill exhibited no observable Cl- uptake. We placed our results and previously reported results in a phylogenetic context and concluded that differences in Cl- uptake mechanisms among groups of freshwater fishes may explain, in large part, the wide range of sensitivity to environmental NO2-. NO2- uptake determinations may also prove to be an easy screening method when studying the phylogenetic distribution and nature of Cl- uptake mechanisms in the gills of fishes. PMID:15667081

  17. Ab initio Potential Energy Surface for H-H2

    NASA Technical Reports Server (NTRS)

    Partridge, Harry; Bauschlicher, Charles W., Jr.; Stallcop, James R.; Levin, Eugene

    1993-01-01

    Ab initio calculations employing large basis sets are performed to determine an accurate potential energy surface for H-H2 interactions for a broad range of separation distances. At large distances, the spherically averaged potential determined from the calculated energies agrees well with the corresponding results determined from dispersion coefficients; the van der Waals well depth is predicted to be 75 +/- (mu)E(sub h). Large basis sets have also been applied to reexamine the accuracy of theoretical repulsive potential energy surfaces. Multipolar expansions of the computed H-H2 potential energy surface are reported for four internuclear separation distances (1.2, 1.401, 1.449, and 1.7a(sub 0) of the hydrogen molecule. The differential elastic scattering cross section calculated from the present results is compared with the measurements from a crossed beam experiment.

  18. The implementation of ab initio quantum chemistry calculations on transporters.

    PubMed

    Cooper, M D; Hillier, I H

    1991-06-01

    The RHF and geometry optimization sections of the ab initio quantum chemistry code, GAMESS, have been optimized for a network of parallel microprocessors, Inmos T800-20 transputers, using both indirect and direct SCF techniques. The results indicate great scope for implementation of such codes on small parallel computer systems, very high efficiencies having been achieved, particularly in the cases of direct SCF and geometry optimization with large basis sets. The work, although performed upon one particular parallel system, the Meiko Computing Surface, is applicable to a wide range of parallel systems with both shared and distributed memory. PMID:1919615

  19. A geometrical correction for the inter- and intra-molecular basis set superposition error in Hartree-Fock and density functional theory calculations for large systems.

    PubMed

    Kruse, Holger; Grimme, Stefan

    2012-04-21

    chemistry yields MAD=0.68 kcal/mol, which represents a huge improvement over plain B3LYP/6-31G* (MAD=2.3 kcal/mol). Application of gCP-corrected B97-D3 and HF-D3 on a set of large protein-ligand complexes prove the robustness of the method. Analytical gCP gradients make optimizations of large systems feasible with small basis sets, as demonstrated for the inter-ring distances of 9-helicene and most of the complexes in Hobza's S22 test set. The method is implemented in a freely available FORTRAN program obtainable from the author's website. PMID:22519309

  20. Crystal structure of the large fragment of Thermus aquaticus DNA polymerase I at 2.5-A resolution: structural basis for thermostability.

    PubMed

    Korolev, S; Nayal, M; Barnes, W M; Di Cera, E; Waksman, G

    1995-09-26

    The crystal structure of the large fragment of the Thermus aquaticus DNA polymerase (Klentaq1), determined at 2.5-A resolution, demonstrates a compact two-domain architecture. The C-terminal domain is identical in fold to the equivalent region of the Klenow fragment of Escherichia coli DNA polymerase I (Klenow pol I). Although the N-terminal domain of Klentaq1 differs greatly in sequence from its counterpart in Klenow pol I, it has clearly evolved from a common ancestor. The structure of Klentaq1 reveals the strategy utilized by this protein to maintain activity at high temperatures and provides the structural basis for future improvements of the enzyme. PMID:7568114

  1. Optical and other material properties of SiO2 from ab initio studies

    NASA Astrophysics Data System (ADS)

    Warmbier, Robert; Mohammed, Faris; Quandt, Alexander

    2014-07-01

    The optical properties of photonic devices largely depend on the dielectric properties of the underlying materials. We apply modern ab initio methods to study crystalline SiO2 phases, which serve as toy models for amorphous glass. We discuss the dielectric response from the infrared to the VIS/UV, which is crucial for glass based photonic applications. Low density silica, like cristobalite, may provide a good basis for high transmission optical devices.

  2. Memories of AB

    NASA Astrophysics Data System (ADS)

    Vaks, V. G.

    2013-06-01

    I had the good fortune to be a student of A. B. Migdal - AB, as we called him in person or in his absence - and to work in the sector he headed at the Kurchatov Institute, along with his other students and my friends, including Vitya Galitsky, Spartak Belyayev and Tolya Larkin. I was especially close with AB in the second half of the 1950s, the years most important for my formation, and AB's contribution to this formation was very great. To this day, I've often quoted AB on various occasions, as it's hard to put things better or more precisely than he did; I tell friends stories heard from AB, because these stories enhance life as AB himself enhanced it; my daughter is named Tanya after AB's wife Tatyana Lvovna, and so on. In what follows, I'll recount a few episodes in my life in which AB played an important or decisive role, and then will share some other memories of AB...

  3. Accurate Ab Initio and Template-Based Prediction of Short Intrinsically-Disordered Regions by Bidirectional Recurrent Neural Networks Trained on Large-Scale Datasets

    PubMed Central

    Volpato, Viola; Alshomrani, Badr; Pollastri, Gianluca

    2015-01-01

    Intrinsically-disordered regions lack a well-defined 3D structure, but play key roles in determining the function of many proteins. Although predictors of disorder have been shown to achieve relatively high rates of correct classification of these segments, improvements over the the years have been slow, and accurate methods are needed that are capable of accommodating the ever-increasing amount of structurally-determined protein sequences to try to boost predictive performances. In this paper, we propose a predictor for short disordered regions based on bidirectional recurrent neural networks and tested by rigorous five-fold cross-validation on a large, non-redundant dataset collected from MobiDB, a new comprehensive source of protein disorder annotations. The system exploits sequence and structural information in the forms of frequency profiles, predicted secondary structure and solvent accessibility and direct disorder annotations from homologous protein structures (templates) deposited in the Protein Data Bank. The contributions of sequence, structure and homology information result in large improvements in predictive accuracy. Additionally, the large scale of the training set leads to low false positive rates, making our systems a robust and efficient way to address high-throughput disorder prediction. PMID:26307973

  4. Reactivity of the binuclear non-heme iron active site of Δ⁹ desaturase studied by large-scale multireference ab initio calculations.

    PubMed

    Chalupský, Jakub; Rokob, Tibor András; Kurashige, Yuki; Yanai, Takeshi; Solomon, Edward I; Rulíšek, Lubomír; Srnec, Martin

    2014-11-12

    The results of density matrix renormalization group complete active space self-consistent field (DMRG-CASSCF) and second-order perturbation theory (DMRG-CASPT2) calculations are presented on various structural alternatives for the O-O and first C-H activating step of the catalytic cycle of the binuclear nonheme iron enzyme Δ(9) desaturase. This enzyme is capable of inserting a double bond into an alkyl chain by double hydrogen (H) atom abstraction using molecular O2. The reaction step studied here is presumably associated with the highest activation barrier along the full pathway; therefore, its quantitative assessment is of key importance to the understanding of the catalysis. The DMRG approach allows unprecedentedly large active spaces for the explicit correlation of electrons in the large part of the chemically important valence space, which is apparently conditio sine qua non for obtaining well-converged reaction energetics. The derived reaction mechanism involves protonation of the previously characterized 1,2-μ peroxy Fe(III)Fe(III) (P) intermediate to a 1,1-μ hydroperoxy species, which abstracts an H atom from the C10 site of the substrate. An Fe(IV)-oxo unit is generated concomitantly, supposedly capable of the second H atom abstraction from C9. In addition, several popular DFT functionals were compared to the computed DMRG-CASPT2 data. Notably, many of these show a preference for heterolytic C-H cleavage, erroneously predicting substrate hydroxylation. This study shows that, despite its limitations, DMRG-CASPT2 is a significant methodological advancement toward the accurate computational treatment of complex bioinorganic systems, such as those with the highly open-shell diiron active sites. PMID:25313991

  5. Towards Accurate Ab Initio Predictions of the Spectrum of Methane

    NASA Technical Reports Server (NTRS)

    Schwenke, David W.; Kwak, Dochan (Technical Monitor)

    2001-01-01

    We have carried out extensive ab initio calculations of the electronic structure of methane, and these results are used to compute vibrational energy levels. We include basis set extrapolations, core-valence correlation, relativistic effects, and Born- Oppenheimer breakdown terms in our calculations. Our ab initio predictions of the lowest lying levels are superb.

  6. Ab initio analytical Raman intensities for periodic systems through a coupled perturbed Hartree-Fock/Kohn-Sham method in an atomic orbital basis. II. Validation and comparison with experiments

    NASA Astrophysics Data System (ADS)

    Maschio, Lorenzo; Kirtman, Bernard; Rérat, Michel; Orlando, Roberto; Dovesi, Roberto

    2013-10-01

    In this work, we validate a new, fully analytical method for calculating Raman intensities of periodic systems, developed and presented in Paper I [L. Maschio, B. Kirtman, M. Rérat, R. Orlando, and R. Dovesi, J. Chem. Phys. 139, 164101 (2013)]. Our validation of this method and its implementation in the CRYSTAL code is done through several internal checks as well as comparison with experiment. The internal checks include consistency of results when increasing the number of periodic directions (from 0D to 1D, 2D, 3D), comparison with numerical differentiation, and a test of the sum rule for derivatives of the polarizability tensor. The choice of basis set as well as the Hamiltonian is also studied. Simulated Raman spectra of α-quartz and of the UiO-66 Metal-Organic Framework are compared with the experimental data.

  7. An expanded calibration study of the explicitly correlated CCSD(T)-F12b method using large basis set standard CCSD(T) atomization energies

    NASA Astrophysics Data System (ADS)

    Feller, David; Peterson, Kirk A.

    2013-08-01

    The effectiveness of the recently developed, explicitly correlated coupled cluster method CCSD(T)-F12b is examined in terms of its ability to reproduce atomization energies derived from complete basis set extrapolations of standard CCSD(T). Most of the standard method findings were obtained with aug-cc-pV7Z or aug-cc-pV8Z basis sets. For a few homonuclear diatomic molecules it was possible to push the basis set to the aug-cc-pV9Z level. F12b calculations were performed with the cc-pVnZ-F12 (n = D, T, Q) basis set sequence and were also extrapolated to the basis set limit using a Schwenke-style, parameterized formula. A systematic bias was observed in the F12b method with the (VTZ-F12/VQZ-F12) basis set combination. This bias resulted in the underestimation of reference values associated with small molecules (valence correlation energies <0.5 Eh) and an even larger overestimation of atomization energies for bigger systems. Consequently, caution should be exercised in the use of F12b for high accuracy studies. Root mean square and mean absolute deviation error metrics for this basis set combination were comparable to complete basis set values obtained with standard CCSD(T) and the aug-cc-pVDZ through aug-cc-pVQZ basis set sequence. However, the mean signed deviation was an order of magnitude larger. Problems partially due to basis set superposition error were identified with second row compounds which resulted in a weak performance for the smaller VDZ-F12/VTZ-F12 combination of basis sets.

  8. Computational Evidence of Inversion of (1)La and (1)Lb-Derived Excited States in Naphthalene Excimer Formation from ab Initio Multireference Theory with Large Active Space: DMRG-CASPT2 Study.

    PubMed

    Shirai, Soichi; Kurashige, Yuki; Yanai, Takeshi

    2016-05-10

    The naphthalene molecule has two important lowest-lying singlet excited states, denoted (1)La and (1)Lb. Association of the excited and ground state monomers yields a metastable excited dimer (excimer), which emits characteristic fluorescence. Here, we report a first computational result based on ab initio theory to corroborate that the naphthalene excimer fluorescence is (1)La parentage, resulting from inversion of (1)La and (1)Lb-derived dimer states. This inversion was hypothesized by earlier experimental studies; however, it has not been confirmed rigorously. In this study, the advanced multireference (MR) theory based on the density matrix renormalization group that enables using unprecedented large-size active space for describing significant electron correlation effects is used to provide accurate potential energy curves (PECs) of the excited states. The results evidenced the inversion of the PECs and accurately predicted transition energies for excimer fluorescence and monomer absorption. Traditional MR calculations with smaller active spaces and single-reference theory calculations exhibit serious inconsistencies with experimental observations. PMID:27082241

  9. Ab initio potential energy surface and rovibrational states of HBO

    NASA Astrophysics Data System (ADS)

    Ha, Tae-Kyu; Makarewicz, Jan

    1999-01-01

    The potential energy surface describing the large-amplitude motion of H around the BO core in the HBO molecule has been determined from ab initio calculations. This surface has been sampled by a set of 170 grid points from a two-dimensional space defined by the stretching and the bending coordinates of the H nucleus. At each grid point, the BO bond length has been optimized using the second-order Møller-Plesset perturbation theory with the basis set aug-cc-pVTZ. The surface has a local minimum for the linear as well as the bent configuration of HBO. A low energy barrier to the linear configuration BOH causes a large-amplitude motion and a strong rovibrational interaction in the molecule. Its rovibrational dynamics is different from the dynamics in bent or quasilinear triatomic molecules.

  10. Spatial distributed hydrological-hydraulic modelling as a basis for the evaluation of flood mitigation measures in a large alpine catchment

    NASA Astrophysics Data System (ADS)

    Gems, B.; Achleitner, St.; Huttenlau, M.; Thieken, A.; Aufleger, M.

    2009-04-01

    Background Precipitation and discharge characteristics have changed in alpine regions within the last decades. Moreover, the changes in the characteristic and superposition of heavy precipitation and snow melt processes as well as modified land use patterns have induced a higher frequency of flooding events. To cope with modified flow characteristics and increasing flow peaks, a sustainable flood control management has to verify, adapt and enhance existing protection measurements. Structures with multifunctional purpose for flood protection as well as hydropower utilization are possible and represent appropriate measurements. Especially when considering large-scale catchment areas, a spatial distributed rainfall runoff modelling including mitigation measures is of importance. Situation The presented contribution deals with the integrated evaluation and upgrade of flood protection structures in an alpine valley located in the Tyrol (Austria). The investigated catchment has an area of ~660 km² and a total main stream length of 22 km. It is situated to the south of the Inntal. It's elevation ranges from 1100 m.a.s.l. up to 3774 m.a.s.l.. Roughly 20% of the total catchment area is permanently glaciated. The used modelling approach incorporates rainfall-runoff models as well as 1D- and 2D-numerical river flow models. Sediment transport processes are going to be considered respectively. More than 30 tributaries to the main receiving water course are modelled using the SCS-procedure. Totally 20 rain gauges with temporal resolutions of at least one hour being available in and in close proximity to the catchment are used as input to the models. Additional five gauges with daily measurements are used. Torrent control measures and sediment transport on the alluvial cone thereby are included using an empirical approach. The 1D-numerical model set up represents the main river in the valley. The dynamic inputs to the model are obtained from the hydrological models. Both - the

  11. The innovative use of a large-scale industry biomedical consortium to research the genetic basis of drug induced serious adverse events.

    PubMed

    Holden, Arthur L

    2007-01-01

    The International Serious Adverse Event Consortium (SAEC) is a pharmaceutical industry and FDA led international (501 c3 non-profit) consortium, focused on identifying and validating DNA-variants useful in predicting the risk of drug induced, rare serious adverse events (SAEs). As such, it functions with the explicit purpose of enhancing the 'public good'. Its members are (i) organizations engaged principally in the business of discovering, developing and marketing pharmaceutical products, or (ii) a charitable, governmental, or other non-profit organization with an interest in researching the molecular basis of drug response.Drug-induced, rare SAEs present significant health issues for patients; and pose challenges for the safe use of approved drugs and the development of new drugs. Examples of drug-induced, rare SAEs include hepatotoxicity, QT prolongation, rhabdomyolosis, serious skin rashes (e.g. SJS), edema, acute renal failure, acute hypersensitivity, anemias/neutropenias, excessive weigh gain, retinopathy, vasculitis, among others. The rarity of such drug induced SAEs and the absence of effective government surveillance/research networks, makes it extremely difficult for any one company or research entity to accrue enough SAE cases and controls to conduct effective whole genome studies. Central to the notion of the SAEC is industry, government and health care providers can join forces to make use of a variety of sample and data resources in researching the genetic basis of these events.The purpose of the SAEC is threefold:•To carry out research directed toward the discovery of DNA-variants clinically useful in understanding and predicting the risk of drug induced serious adverse events and similar scientific research.•To ensure the widespread availability of the results of such research to the scientific research community and the public at large for no charge through publication and web-based methods; and•To educate the scientific research and medical

  12. Ab initio methods for nuclear properties - a computational physics approach

    NASA Astrophysics Data System (ADS)

    Maris, Pieter

    2011-04-01

    A microscopic theory for the structure and reactions of light nuclei poses formidable challenges for high-performance computing. Several ab-initio methods have now emerged that provide nearly exact solutions for some nuclear properties. The ab-initio no-core full configuration (NCFC) approach is based on basis space expansion methods and uses Slater determinants of single-nucleon basis functions to express the nuclear wave function. In this approach, the quantum many-particle problem becomes a large sparse matrix eigenvalue problem. The eigenvalues of this matrix give us the binding energies, and the corresponding eigenvectors the nuclear wave functions. These wave functions can be employed to evaluate experimental quantities. In order to reach numerical convergence for fundamental problems of interest, the matrix dimension often exceeds 1 billion, and the number of nonzero matrix elements may saturate available storage on present-day leadership class facilities. I discuss different strategies for distributing and solving this large sparse matrix on current multicore computer architectures, including methods to deal with with memory bottleneck. Several of these strategies have been implemented in the code MFDn, which is a parallel fortran code for nuclear structure calculations. I will show scaling behavior and compare the performance of the pure MPI version with the hybrid MPI/OpenMP code on Cray XT4 and XT5 platforms. For large core counts (typically 5,000 and above), the hybrid version is more efficient than pure MPI. With this code, we have been able to predict properties of the unstable nucleus 14F, which have since been confirmed by experiments. I will also give an overview of other recent results for nuclei in the A = 6 to 16 range with 2- and 3-body interactions. Supported in part by US DOE Grant DE-FC02-09ER41582.

  13. Omega-AB

    Energy Science and Technology Software Center (ESTSC)

    2007-05-01

    A hierarchical, modular modeling environment for hybrid simulations of sequential-modular, systems dynamics, discrete-event, and agent-based paradigms Omega-AB models contain a hierarchically-defined module tree that specifies the execution logic for the simulation, and a multi-network graph that defines the environment within which the simulation occurs. Modules are the fundamental buildinig blocks of an Omega-AB model and can define anything from a basic mathematical operation to a complex behavioral response model. Modules rely on the "plug-in" conceptmore » which allows developers to build independent module libraries that are gathered, linked, and instantiated by the Omega-AB engine at run time. Inter-module communication occurs through two complimentary systems: pull-based "ports" for general computation patterns and push-based "plugs" for event processing. The simulation environment is an abstract graph of nodes and links. Agents (module sub-trees headed up by an Agent module) reside at nodes and relate to their neighbors through typed links. To facilitate the construction and visualization of complex, interacting networks with dramatically different structure, Omega-AB provides a system for organizing the nodes into hierarchica trees that describe "slices" of the overall network.« less

  14. Ab Initio Potential Energy Surface for H-H2

    NASA Technical Reports Server (NTRS)

    Patridge, Harry; Bauschlicher, Charles W., Jr.; Stallcop, James R.; Levin, Eugene

    1993-01-01

    Ab initio calculations employing large basis sets are performed to determine an accurate potential energy surface for H-H2 interactions for a broad range of separation distances. At large distances, the spherically averaged potential determined from the calculated energies agrees well with the corresponding results determined from dispersion coefficients; the van der Waals well depth is predicted to be 75 +/- 3 micro E(h). Large basis sets have also been applied to reexamine the accuracy of theoretical repulsive potential energy surfaces (25-70 kcal/mol above the H-H2 asymptote) at small interatomic separations; the Boothroyd, Keogh, Martin, and Peterson (BKMP) potential energy surface is found to agree with results of the present calculations within the expected uncertainty (+/- 1 kcal/mol) of the fit. Multipolar expansions of the computed H-H2 potential energy surface are reported for four internuclear separation distances (1.2, 1.401, 1.449, and 1.7a(0)) of the hydrogen molecule. The differential elastic scattering cross section calculated from the present results is compared with the measurements from a crossed beam experiment.

  15. In pursuit of the ab initio limit for conformational energy prototypes

    NASA Astrophysics Data System (ADS)

    Császár, Attila G.; Allen, Wesley D.; Schaefer, Henry F.

    1998-06-01

    The convergence of ab initio predictions to the one- and n-particle limits has been systematically explored for several conformational energy prototypes: the inversion barriers of ammonia, water, and isocyanic acid, the torsional barrier of ethane, the E/Z rotamer separation of formic acid, and the barrier to linearity of silicon dicarbide. Explicit ab initio results were obtained with atomic-orbital basis sets as large as [7s6p5d4f3g2h1i/6s5p4d3f2g1h] and electron correlation treatments as extensive as fifth-order Møller-Plesset perturbation theory (MP5), the full coupled-cluster method through triple excitations (CCSDT), and Brueckner doubles theory including perturbational corrections for both triple and quadruple excitations [BD(TQ)]. Subsequently, basis set and electron correlation extrapolation schemes were invoked to gauge any further variations in arriving at the ab initio limit. Physical effects which are tacitly neglected in most theoretical work have also been quantified by computations of non-Born-Oppenheimer (BODC), relativistic, and core correlation shifts of relative energies. Instructive conclusions are drawn for the pursuit of spectroscopic accuracy in theoretical conformational analyses, and precise predictions for the key energetic quantities of the molecular prototypes are advanced.

  16. Ab initio potentials of F+Li{sub 2} accessible at ultracold temperatures

    SciTech Connect

    Wright, K. W. A.; Lane, Ian C.

    2010-09-15

    Ab initio calculations for the strongly exoergic Li{sub 2}+F harpoon reaction are presented using density-functional theory, complete active space self-consistent field, and multireference configuration interaction methods to argue that this reaction would be an ideal candidate for investigation with ultracold molecules. The lowest six states are calculated with the aug-correlation-consistent polarized valence triple-zeta basis set and at least two can be accessed by a ground rovibronic Li{sub 2} molecule with zero collision energy at all reaction geometries. The large reactive cross section (characteristic of harpoon reactions) and chemiluminescent products are additional attractive features of these reactions.

  17. The ab initio potential energy surface and spectroscopic constants of HOCl

    NASA Astrophysics Data System (ADS)

    Koput, Jacek; Peterson, Kirk A.

    1998-02-01

    The potential energy surface of hypochlorous acid, HOCl, has been determined from large-scale ab initio calculations using the coupled-cluster method CCSD(T), with basis sets of quadruple- and quintuple-zeta quality. The effect of core-electron correlation on the calculated structural parameters has been investigated. The vibrational-rotational energy levels of the three isotopic species of HOCl have then been calculated using the variational method and have been further characterized by the spectroscopic constants determined using the perturbational approach. The spectroscopic constants determined, are found to be in excellent agreement with experimental data.

  18. Ab initio potentials of F+Li2 accessible at ultracold temperatures

    NASA Astrophysics Data System (ADS)

    Wright, K. W. A.; Lane, Ian C.

    2010-09-01

    Ab initio calculations for the strongly exoergic Li2+F harpoon reaction are presented using density-functional theory, complete active space self-consistent field, and multireference configuration interaction methods to argue that this reaction would be an ideal candidate for investigation with ultracold molecules. The lowest six states are calculated with the aug-correlation-consistent polarized valence triple-zeta basis set and at least two can be accessed by a ground rovibronic Li2 molecule with zero collision energy at all reaction geometries. The large reactive cross section (characteristic of harpoon reactions) and chemiluminescent products are additional attractive features of these reactions.

  19. An ab initio molecular dynamics study of the roaming mechanism of the H2+HOC+ reaction

    NASA Astrophysics Data System (ADS)

    Yu, Hua-Gen

    2011-08-01

    We report here a direct ab initio molecular dynamics study of the p-/o-H2+HOC+ reaction on the basis of the accurate SAC-MP2 potential energy surface. The quasi-classical trajectory method was employed. This work largely focuses on the study of reaction mechanisms. A roaming mechanism was identified for this molecular ion-molecule reaction. The driving forces behind the roaming mechanism were thoroughly investigated by using a trajectory dynamics approach. In addition, the thermal rate coefficients of the H2+HOC+ reaction were calculated in the temperature range [25, 300] K and are in good agreement with experiments.

  20. Single particle calculations for a Woods-Saxon potential with triaxial deformations, and large Cartesian oscillator basis (TRIAXIAL 2014, Third version of the code Triaxial)

    NASA Astrophysics Data System (ADS)

    Mohammed-Azizi, B.; Medjadi, D. E.

    2014-11-01

    Theory and FORTRAN program of the first version of this code (TRIAXIAL) have already been described in detail in Computer Physics Comm. 156 (2004) 241-282. A second version of this code (TRIAXIAL 2007) has been given in CPC 176 (2007) 634-635. The present FORTRAN program is the third version (TRIAXIAL 2014) of the same code. Now, It is written in free format. As the former versions, this FORTRAN program solves the same Schrodinger equation of the independent particle model of the atomic nucleus with the same method. However, the present version is much more convenient. In effect, it is characterized by the fact that the eigenvalues and the eigenfunctions can be given by specific subroutines. The latters did not exist in the old versions (2004 and 2007). In addition, it is to be noted that in the previous versions, the eigenfunctions were only given by their coefficients of their expansion onto the harmonic oscillator basis. This method is needed in some cases. But in other cases, it is preferable to treat the eigenfunctions directly in configuration space. For this reason, we have implemented an additional subroutine for this task. Some other practical subroutines have also been implemented. Moreover, eigenvalues and eigenfunctions are recorded onto several files. All these new features of the code and some important aspects of its structure are explained in the document ‘Triaxial2014 use.pdf’. Catalogue identifier: ADSK_v3_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADSK_v3_0.html Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 13672 No. of bytes in distributed program, including test data, etc.: 217598 Distribution format: tar.gz Programming language: FORTRAN 77/90 (double precision). Computer: PC. Pentium 4, 2600MHz and beyond. Operating system: WINDOWS XP

  1. Solution of the Dirac equation using the Rayleigh-Ritz method: Flexible basis coupling large and small components. Results for one-electron systems

    NASA Astrophysics Data System (ADS)

    Baǧcı, A.; Hoggan, P. E.

    2016-07-01

    An algebraic solution of the Dirac equation is reinvestigated. Slater-type spinor orbitals and their corresponding system of differential equations are defined in two- and four-component formalism. They describe the radial function in components of the wave function of the Dirac equation solution to high accuracy. They constitute the matrix elements arising in a generalized eigenvalue equation. These terms are evaluated through prolate spheroidal coordinates. The corresponding integrals are calculated by the numerical global-adaptive method taking into account the Gauss-Kronrod numerical integration extension. Sample calculations are performed using flexible basis sets generated with both signs of the relativistic angular momentum quantum number κ . Applications to one-electron atoms and diatomics are detailed. Variationally optimum values for orbital parameters are obtained at given nuclear separation. Methods discussed in this work are capable of yielding highly accurate relativistic two-center integrals for all ranges of orbital parameters. This work provides an efficient way to overcome the problems that arise in relativistic calculations.

  2. Solution of the Dirac equation using the Rayleigh-Ritz method: Flexible basis coupling large and small components. Results for one-electron systems.

    PubMed

    Bağcı, A; Hoggan, P E

    2016-07-01

    An algebraic solution of the Dirac equation is reinvestigated. Slater-type spinor orbitals and their corresponding system of differential equations are defined in two- and four-component formalism. They describe the radial function in components of the wave function of the Dirac equation solution to high accuracy. They constitute the matrix elements arising in a generalized eigenvalue equation. These terms are evaluated through prolate spheroidal coordinates. The corresponding integrals are calculated by the numerical global-adaptive method taking into account the Gauss-Kronrod numerical integration extension. Sample calculations are performed using flexible basis sets generated with both signs of the relativistic angular momentum quantum number κ. Applications to one-electron atoms and diatomics are detailed. Variationally optimum values for orbital parameters are obtained at given nuclear separation. Methods discussed in this work are capable of yielding highly accurate relativistic two-center integrals for all ranges of orbital parameters. This work provides an efficient way to overcome the problems that arise in relativistic calculations. PMID:27575231

  3. ABS plastic RPCs

    SciTech Connect

    Ables, E.; Bionta, R.; Olson, H.; Ott, L.; Parker, E.; Wright, D.; Wuest, C

    1996-02-01

    After investigating a number of materials, we discovered that an ABS plastic doped with a conducting polymer performs well as the resistive electrode in a narrow gap RPC (resistive plate chamber). Operating in the streamer mode, we find efficiencies of 90-96% with low noise and low strip multiplicities. We have also studied a variety of operating gases and found that a mixture containing SF{sub 6}, a non-ozone depleting gas, argon and isobutane gives good streamer mode performance, even with isobutane concentrations of 20% or less.

  4. Single particle calculations for a Woods-Saxon potential with triaxial deformations, and large Cartesian oscillator basis (TRIAXIAL 2014, Third version of the code Triaxial)

    NASA Astrophysics Data System (ADS)

    Mohammed-Azizi, B.; Medjadi, D. E.

    2014-11-01

    Theory and FORTRAN program of the first version of this code (TRIAXIAL) have already been described in detail in Computer Physics Comm. 156 (2004) 241-282. A second version of this code (TRIAXIAL 2007) has been given in CPC 176 (2007) 634-635. The present FORTRAN program is the third version (TRIAXIAL 2014) of the same code. Now, It is written in free format. As the former versions, this FORTRAN program solves the same Schrodinger equation of the independent particle model of the atomic nucleus with the same method. However, the present version is much more convenient. In effect, it is characterized by the fact that the eigenvalues and the eigenfunctions can be given by specific subroutines. The latters did not exist in the old versions (2004 and 2007). In addition, it is to be noted that in the previous versions, the eigenfunctions were only given by their coefficients of their expansion onto the harmonic oscillator basis. This method is needed in some cases. But in other cases, it is preferable to treat the eigenfunctions directly in configuration space. For this reason, we have implemented an additional subroutine for this task. Some other practical subroutines have also been implemented. Moreover, eigenvalues and eigenfunctions are recorded onto several files. All these new features of the code and some important aspects of its structure are explained in the document ‘Triaxial2014 use.pdf’. Catalogue identifier: ADSK_v3_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADSK_v3_0.html Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 13672 No. of bytes in distributed program, including test data, etc.: 217598 Distribution format: tar.gz Programming language: FORTRAN 77/90 (double precision). Computer: PC. Pentium 4, 2600MHz and beyond. Operating system: WINDOWS XP

  5. Ab inito investigation of hydrodesulferization

    SciTech Connect

    Tilson, J.L.; Marshall, C.L.; Brenner, J.R.

    1997-12-31

    We utilize MPP and vector computers to model the interaction of large sulfur-containing species bonded with hydrodesulfurization (HDS) catalysts. This work is made possible by the availability large aggregate memory, parallel computers. The ability of modern non-traditional computers to solve large-scale scientific problems has been demonstrated. This success is accomplished, in part, by access to portable low- and user-level software tools which exhibit good control over NUMA. ab initio SCF methods are used to characterize the energies of adsorption of poly-aromatic, sulfur-containing hydrocarbons onto a series of molecular-based models of MoS2. These hydrocarbon include thiophene (TP), benzothiophene (BT) and dibenzothiophene (DBT) which are representative of heavy crude feedstocks. Our work attempts to ascertain if a consistent model of adsorption onto a MoS2 surface exists and to determine how the electronic and steric configuration of each species affect their energies of interaction with the metal surface.

  6. MELCOR 1.8.2 assessment: Aerosol experiments ABCOVE AB5, AB6, AB7, and LACE LA2

    SciTech Connect

    Souto, F.J.; Haskin, F.E.; Kmetyk, L.N.

    1994-10-01

    The MELCOR computer code has been used to model four of the large-scale aerosol behavior experiments conducted in the Containment System Test Facility (CSTF) vessel. Tests AB5, AB6 and AB7 of the ABCOVE program simulate the dry aerosol conditions during a hypothetical severe accident in an LMFBR. Test LA2 of the LACE program simulates aerosol behavior in a condensing steam environment during a postulated severe accident in an LWR with failure to isolate the containment. The comparison of code results to experimental data show that MELCOR is able to correctly predict most of the thermal-hydraulic results in the four tests. MELCOR predicts reasonably well the dry aerosol behavior of the ABCOVE tests, but significant disagreements are found in the aerosol behavior modelling for the LA2 experiment. These results tend to support some of the concerns about the MELCOR modelling of steam condensation onto aerosols expressed in previous works. During these analyses, a limitation in the MELCOR input was detected for the specification of the aerosol parameters for more than one component. A Latin Hypercube Sampling (LHS) sensitivity study of the aerosol dynamic constants is presented for test AB6. The study shows the importance of the aerosol shape factors in the aerosol deposition behavior, and reveals that MELCOR input/output processing is highly labor intensive for uncertainty and sensitivity analyses based on LHS.

  7. Phylogenetic and environmental diversity of DsrAB-type dissimilatory (bi)sulfite reductases

    PubMed Central

    Müller, Albert Leopold; Kjeldsen, Kasper Urup; Rattei, Thomas; Pester, Michael; Loy, Alexander

    2015-01-01

    The energy metabolism of essential microbial guilds in the biogeochemical sulfur cycle is based on a DsrAB-type dissimilatory (bi)sulfite reductase that either catalyzes the reduction of sulfite to sulfide during anaerobic respiration of sulfate, sulfite and organosulfonates, or acts in reverse during sulfur oxidation. Common use of dsrAB as a functional marker showed that dsrAB richness in many environments is dominated by novel sequence variants and collectively represents an extensive, largely uncharted sequence assemblage. Here, we established a comprehensive, manually curated dsrAB/DsrAB database and used it to categorize the known dsrAB diversity, reanalyze the evolutionary history of dsrAB and evaluate the coverage of published dsrAB-targeted primers. Based on a DsrAB consensus phylogeny, we introduce an operational classification system for environmental dsrAB sequences that integrates established taxonomic groups with operational taxonomic units (OTUs) at multiple phylogenetic levels, ranging from DsrAB enzyme families that reflect reductive or oxidative DsrAB types of bacterial or archaeal origin, superclusters, uncultured family-level lineages to species-level OTUs. Environmental dsrAB sequences constituted at least 13 stable family-level lineages without any cultivated representatives, suggesting that major taxa of sulfite/sulfate-reducing microorganisms have not yet been identified. Three of these uncultured lineages occur mainly in marine environments, while specific habitat preferences are not evident for members of the other 10 uncultured lineages. In summary, our publically available dsrAB/DsrAB database, the phylogenetic framework, the multilevel classification system and a set of recommended primers provide a necessary foundation for large-scale dsrAB ecology studies with next-generation sequencing methods. PMID:25343514

  8. Features of the application of the Monte Carlo method to calculations for large RBMK reactors and to model correction on the basis of data from in-core detectors

    SciTech Connect

    Ivanov, I. E. Schukin, N. V.; Bychkov, S. A.; Druzhinin, V. E.; Lysov, D. A.; Shmonin, Yu. V.; Gurevich, M. I.

    2014-12-15

    Statistical errors in sampling neutron fields in physically large systems like an RBMK are analyzed both qualitatively and quantitatively. Recommendations concerning the choice of parameters for calculations are given. A new procedure for Monte Carlo RBMK calculations with model corrections on the basis of data from in-core detectors is proposed. Dedicated software based on the CUDA software and hardware platform is developed for computational research. Results of testing the procedure and software in question via calculations for real RBMK reactors are discussed.

  9. Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics

    SciTech Connect

    Makhov, Dmitry V.; Shalashilin, Dmitrii V.; Glover, William J.; Martinez, Todd J.

    2014-08-07

    We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as “cloning,” in analogy to the “spawning” procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, “trains,” as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions.

  10. Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics

    NASA Astrophysics Data System (ADS)

    Makhov, Dmitry V.; Glover, William J.; Martinez, Todd J.; Shalashilin, Dmitrii V.

    2014-08-01

    We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as "cloning," in analogy to the "spawning" procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, "trains," as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions.

  11. Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics.

    PubMed

    Makhov, Dmitry V; Glover, William J; Martinez, Todd J; Shalashilin, Dmitrii V

    2014-08-01

    We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as "cloning," in analogy to the "spawning" procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, "trains," as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions. PMID:25106573

  12. The computation of ionization potentials for second-row elements by ab initio and density functional theory methods

    SciTech Connect

    Jursic, B.S.

    1996-12-31

    Up to four ionization potentials of elements from the second-row of the periodic table were computed using the ab initio (HF, MP2, MP3, MP4, QCISD, GI, G2, and G2MP2) and DFT (B3LY, B3P86, B3PW91, XALPHA, HFS, HFB, BLYP, BP86, BPW91, BVWN, XAPLY, XAP86, XAPW91, XAVWN, SLYR SP86, SPW91 and SVWN) methods. In all of the calculations, the large 6-311++G(3df,3pd) gaussian type of basis set was used. The computed values were compared with the experimental results and suitability of the ab initio and DFF methods were discussed, in regard to reproducing the experimental data. From the computed ionization potentials of the second-row elements, it can be concluded that the HF ab initio computation is not capable of reproducing the experimental results. The computed ionization potentials are too low. However, by using the ab initio methods that include electron correlation, the computed IPs are becoming much closer to the experimental values. In all cases, with the exception of the first ionization potential for oxygen, the G2 computation result produces ionization potentials that are indistinguishable from the experimental results.

  13. Calibration study of the CCSD(T)-F12a/b methods for C2 and small hydrocarbons

    NASA Astrophysics Data System (ADS)

    Feller, David; Peterson, Kirk A.; Hill, J. Grant

    2010-11-01

    Explicitly correlated CCSD(T)-F12a/b methods combined with basis sets specifically designed for this technique have been tested for their ability to reproduce standard CCSD(T) benchmark data covering 16 small molecules composed of hydrogen and carbon. The standard method calibration set was obtained with very large one-particle basis sets, including some aug-cc-pV7Z and aug-cc-pV8Z results. Whenever possible, the molecular properties (atomization energies, structures, and harmonic frequencies) were extrapolated to the complete basis set limit in order to facilitate a direct comparison of the standard and explicitly correlated approaches without ambiguities arising from the use of different basis sets. With basis sets of triple-ζ quality or better, the F12a variant was found to overshoot the presumed basis set limit, while the F12b method converged rapidly and uniformly. Extrapolation of F12b energies to the basis set limit was found to be very effective at reproducing the best standard method atomization energies. Even extrapolations based on the small cc-pVDZ-F12/cc-pVTZ-F12 combination proved capable of a mean absolute deviation of 0.20 kcal/mol. The accuracy and simultaneous cost savings of the F12b approach are such that it should enable high quality property calculations to be performed on chemical systems that are too large for standard CCSD(T).

  14. Ab initio infrared and Raman spectra

    NASA Astrophysics Data System (ADS)

    Fredkin, Donald R.; Komornicki, Andrew; White, Steven R.; Wilson, Kent R.

    1983-06-01

    We discuss several ways in which molecular absorption and scattering spectra can be computed ab initio, from the fundamental constants of nature. These methods can be divided into two general categories. In the first, or sequential, type of approach, one first solves the electronic part of the Schrödinger equation in the Born-Oppenheimer approximation, mapping out the potential energy, dipole moment vector (for infrared absorption) and polarizability tensor (for Raman scattering) as functions of nuclear coordinates. Having completed the electronic part of the calculation, one then solves the nuclear part of the problem either classically or quantum mechanically. As an example of the sequential ab initio approach, the infrared and Raman rotational and vibrational-rotational spectral band contours for the water molecule are computed in the simplest rigid rotor, normal mode approximation. Quantum techniques are used to calculate the necessary potential energy, dipole moment, and polarizability information at the equilibrium geometry. A new quick, accurate, and easy to program classical technique involving no reference to Euler angles or special functions is developed to compute the infrared and Raman band contours for any rigid rotor, including asymmetric tops. A second, or simultaneous, type of ab initio approach is suggested for large systems, particularly those for which normal mode analysis is inappropriate, such as liquids, clusters, or floppy molecules. Then the curse of dimensionality prevents mapping out in advance the complete potential, dipole moment, and polarizability functions over the whole space of nuclear positions of all atoms, and a solution in which the electronic and nuclear parts of the Born-Oppenheimer approximation are simultaneously solved is needed. A quantum force classical trajectory (QFCT) molecular dynamic method, based on linear response theory, is described, in which the forces, dipole moment, and polarizability are computed quantum

  15. An ab initio potential energy surface and dynamics of the Ar+H2+ → ArH + + H reaction

    NASA Astrophysics Data System (ADS)

    Liu, Xinguo; Liu, Huirong; Zhang, Qinggang

    2011-04-01

    An ab initio potential energy surface (PES) for the ground state (1 2A') of the chemical reaction Ar+H2+ → ArH + + H has been constructed from a set of accurate ab initio data, which we have computed using the coupled-cluster theory including all single and double excitations plus perturbative corrections for the triples UCCSD(T) with a large orbital basis set of aug-cc-pV5Z. The new PES has a root-mean-square (rms) error of 0.5341 kcal/mol. The total integral reaction cross-sections have been calculated at three collision energies by means of the quasi-classical trajectory (QCT) calculation based on the new PES and compared with previous TSH results.

  16. GINGA Observations of AB Doradus

    NASA Astrophysics Data System (ADS)

    Vilhu, O.; Tsuru, T.; Collier Cameron, A.

    We report GINGA observations of the pre main sequence star AB Doradus (HD 36705), performed during 8 - 12 January, 1990. Some rotational modulation might be present. four X-ray flares were detected. Three of these events were similar to the EINSTEIN HRI-flare (Vilhu and Linsky, 1987), with decay times around 25 min. The last flare had long rise and decay times (100 min), resembling the EXOSAT flares observed by Collier Cameron et.al. (1988). The mean flare spectrum can be fitted by a thermal bremstrahlung with temperature 5.0 keV, or by a power-law model with photon index 2.2. The 3 upper limit of the Iron line equivalent width in the flare spectrum is 1 keV, weaker than predicted by thermal models. This Iron line anomaly was first discussed in the case of UX Ari by Tsuru et. al. (1989). However, normal equivalent widths can be derived from several EXOSAT spectra of active cool stars (Pallavicini and Tagliaferri, 1990). We discuss the possibility that the continuum from non-thermal electrons (producing also the microwave emission) could occasionally lower the apparent equivalent width. The mechanism works for reasonably low magnetic field strengths and electon power-law indexes. However, a large population of non-thermal electrons is needed (comparable to the thermal one). Stronger magnetic fields could explain the radio emission with less electrons, but then the non-thermal X-ray continuum remains small.

  17. An accurate {ital ab initio} HOCl potential energy surface, vibrational and rotational calculations, and comparison with experiment

    SciTech Connect

    Skokov, S.; Peterson, K.A.; Bowman, J.M.

    1998-08-01

    Accurate {ital ab initio} multireference configuration interaction (CI) calculations with large correlation-consistent basis sets are performed for HOCl. After extrapolation to the complete basis set limit, the {ital ab initio} data are precisely fit to give a semiglobal three-dimensional potential energy surface to describe HOCl{r_arrow}Cl+OH from high overtone excitation of the OH-stretch. The average absolute deviation between the {ital ab initio} and fitted energies is 4.2thinspcm{sup {minus}1} for energies up to 60 kcal/mol relative to the HOCl minimum. Vibrational energies of HOCl including the six overtones of the OH-stretch are computed using a vibrational-Cl method on the fitted potential and also on a slightly adjusted potential. Near-spectroscopic accuracy is obtained using the adjusted potential; the average absolute deviation between theory and experiment for 19 experimentally reported states is 4.8thinspcm{sup {minus}1}. Very good agreement with experiment is also obtained for numerous rotational energies for the ground vibrational state, the ClO-stretch fundamental, and the fifth overtone of the OH-stretch. {copyright} {ital 1998 American Institute of Physics.}

  18. Ab initio potential energy curves of the valence, Rydberg, and ion-pair states of iodine monochloride, ICl

    SciTech Connect

    Kalemos, Apostolos; Prosmiti, Rita

    2014-09-14

    We present for the first time a coherent ab initio study of 39 states of valence, Rydberg, and ion-pair character of the diatomic interhalogen ICl species through large scale multireference variational methods including spin-orbit effects coupled with quantitative basis sets. Various avoided crossings are responsible for a non-adiabatic behaviour creating a wonderful vista for its theoretical description. Our molecular constants are compared with all available experimental data with the aim to assist experimentalists especially in the high energy regime of up to ∼95 000 cm{sup −1}.

  19. Hamiltonian Light-Front Ffield Theory in a Basis Function Approach

    SciTech Connect

    Vary, J.P.; Honkanen, H.; Li, Jun; Maris, P.; Brodsky, S.J.; Harindranath, A.; de Teramond, G.F.; Sternberg, P.; Ng, E.G.; Yang, C.

    2009-05-15

    Hamiltonian light-front quantum field theory constitutes a framework for the non-perturbative solution of invariant masses and correlated parton amplitudes of self-bound systems. By choosing the light-front gauge and adopting a basis function representation, we obtain a large, sparse, Hamiltonian matrix for mass eigenstates of gauge theories that is solvable by adapting the ab initio no-core methods of nuclear many-body theory. Full covariance is recovered in the continuum limit, the infinite matrix limit. There is considerable freedom in the choice of the orthonormal and complete set of basis functions with convenience and convergence rates providing key considerations. Here, we use a two-dimensional harmonic oscillator basis for transverse modes that corresponds with eigensolutions of the soft-wall AdS/QCD model obtained from light-front holography. We outline our approach, present illustrative features of some non-interacting systems in a cavity and discuss the computational challenges.

  20. Radiative lifetimes of spin forbidden a1Δ → X 3Σ- and spin allowed A3Π → X 3Σ- transitions and complete basis set extrapolated ab initio potential energy curves for the ground and excited states of CH-

    NASA Astrophysics Data System (ADS)

    Srivastava, Saurabh; Sathyamurthy, N.

    2012-12-01

    The spin forbidden transition a1Δ → X3Σ- in CH- has been studied using the Breit-Pauli Hamiltonian for a large number of geometries. This transition acquires intensity through spin-orbit coupling with singlet and triplet Π states. The transition moment matrix including more than one singlet and triplet Π states was calculated at the multi-reference configuration interaction/aug-cc-pV6Z level of theory. The computed radiative lifetime of 5.63 s is in good agreement with the experimental (5.9 s) and other theoretical (6.14 s) results. Transition moment values of the spin allowed A3Π → X3Σ- transition have also been calculated at the same level of theory. Calculations show that the corresponding radiative lifetime is considerably low, 2.4 × 10-7 s. Complete basis set extrapolated potential energy curves for the ground state of CH and the ground state and six low lying excited states (a1Δ, b1Σ+, two 3Π, and two 1Π) of CH- are reported. These curves are then used to calculate the vibrational bound states for CH and CH-. The computed electron affinity of CH supports the electron affinity bounds reported by Okumura et al. [J. Chem. Phys. 85, 1971 (1986), 10.1063/1.451140].

  1. Carbon dioxide hydrate phase equilibrium and cage occupancy calculations using ab initio intermolecular potentials.

    PubMed

    Velaga, Srinath C; Anderson, Brian J

    2014-01-16

    Gas hydrate deposits are receiving increased attention as potential locations for CO2 sequestration, with CO2 replacing the methane that is recovered as an energy source. In this scenario, it is very important to correctly characterize the cage occupancies of CO2 to correctly assess the sequestration potential as well as the methane recoverability. In order to predict accurate cage occupancies, the guest–host interaction potential must be represented properly. Earlier, these potential parameters were obtained by fitting to experimental equilibrium data and these fitted parameters do not match with those obtained by second virial coefficient or gas viscosity data. Ab initio quantum mechanical calculations provide an independent means to directly obtain accurate intermolecular potentials. A potential energy surface (PES) between H2O and CO2 was computed at the MP2/aug-cc-pVTZ level and corrected for basis set superposition error (BSSE), an error caused due to the lower basis set, by using the half counterpoise method. Intermolecular potentials were obtained by fitting Exponential-6 and Lennard-Jones 6-12 models to the ab initio PES, correcting for many-body interactions. We denoted this model as the “VAS” model. Reference parameters for structure I carbon dioxide hydrate were calculated using the VAS model (site–site ab initio intermolecular potentials) as Δμ(w)(0) = 1206 ± 2 J/mol and ΔH(w)(0) = 1260 ± 12 J/mol. With these reference parameters and the VAS model, pure CO2 hydrate equilibrium pressure was predicted with an average absolute deviation of less than 3.2% from the experimental data. Predictions of the small cage occupancy ranged from 32 to 51%, and the large cage is more than 98% occupied. The intermolecular potentials were also tested by calculating the pure CO2 density and diffusion of CO2 in water using molecular dynamics simulations. PMID:24328234

  2. Halo nuclei He6 and He8 with the Coulomb-Sturmian basis

    NASA Astrophysics Data System (ADS)

    Caprio, M. A.; Maris, P.; Vary, J. P.

    2014-09-01

    The rapid Gaussian falloff of the oscillator functions at large radius makes them poorly suited for the description of the asymptotic properties of the nuclear wave function, a problem which becomes particularly acute for halo nuclei. We consider an alternative basis for ab initio no-core configuration interaction (NCCI) calculations, built from Coulomb-Sturmian radial functions, allowing for realistic exponential falloff. NCCI calculations are carried out for the neutron halo nuclei He6,8, as well as for the baseline case 4He, with the JISP16 nucleon-nucleon interaction. Estimates are made for the root-mean-square radii of the proton and matter distributions.

  3. Converging Nuclear Magnetic Shielding Calculations with Respect to Basis and System Size in Protein Systems

    PubMed Central

    Hartman, Joshua D.; Neubauer, Thomas J.; Caulkins, Bethany G.; Mueller, Leonard J.; Beran, Gregory J. O.

    2015-01-01

    Ab initio chemical shielding calculations greatly facilitate the interpretation of nuclear magnetic resonance (NMR) chemical shifts in biological systems, but the large sizes of these systems requires approximations in the chemical models used to represent them. Achieving good convergence in the predicted chemical shieldings is necessary before one can unravel how other complex structural and dynamical factors affect the NMR measurements. Here, we investigate how to balance trade-offs between using a better basis set or a larger cluster model for predicting the chemical shieldings of the substrates in two representative examples of protein-substrate systems involving different domains in tryptophan synthase: the N-(4′-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F9) ligand which binds in the α active site, and the 2-aminophenol (2AP) quinonoid intermediate formed in the β active site. We first demonstrate that a chemically intuitive three-layer, locally dense basis model that uses a large basis on the substrate, a medium triple-zeta basis to describe its hydrogen-bonding partners and/or surrounding van derWaals cavity, and a crude basis set for more distant atoms provides chemical shieldings in good agreement with much more expensive large basis calculations. Second, long-range quantum mechanical interactions are important, and one can accurately estimate them as a small-basis correction to larger-basis calculations on a smaller cluster. The combination of these approaches enables one to perform density functional theory NMR chemical shift calculations in protein systems that are well-converged with respect to both basis set and cluster size. PMID:25993979

  4. ABS_PDv14.0 beta

    Energy Science and Technology Software Center (ESTSC)

    2008-07-30

    ABS_PD stands for Adaptive Background Suppression and Peak Detection algorithm. This algorithm uses an iterative process to simultaneously estimate the background and identify peaks in signals where the bakground is slowly varying relative to the peaks to be detected. This is sone on an automated manner. The algorithm can be applied to many physical processes where the desired signal is superimposed on a background. The main advantage of this algorithm is that the background canmore » be variable and doesn't have to be known. All that is required is that a set of basis function,s capa ble of representing the background, needs to be defined. Depending on the signals to be processed, the basis functions can be as simple as low order polynomials. The current algorithm has polynomials built-in and allows for additional basis functions to be defined by the user. An additional advantage is that the algorithm does not use any derivatives in detection of peaks and thus allows for analysis of noisy data. The noise level in the data is automatically taken into account when setting thresholds for peak detection.« less

  5. GAUSSIAN 76: An ab initio Molecular Orbital Program

    DOE R&D Accomplishments Database

    Binkley, J. S.; Whiteside, R.; Hariharan, P. C.; Seeger, R.; Hehre, W. J.; Lathan, W. A.; Newton, M. D.; Ditchfield, R.; Pople, J. A.

    1978-01-01

    Gaussian 76 is a general-purpose computer program for ab initio Hartree-Fock molecular orbital calculations. It can handle basis sets involving s, p and d-type Gaussian functions. Certain standard sets (STO-3G, 4-31G, 6-31G*, etc.) are stored internally for easy use. Closed shell (RHF) or unrestricted open shell (UHF) wave functions can be obtained. Facilities are provided for geometry optimization to potential minima and for limited potential surface scans.

  6. The study of molecular spectroscopy by ab initio methods

    NASA Technical Reports Server (NTRS)

    Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.

    1991-01-01

    This review illustrates the potential of theory for solving spectroscopic problems. The accuracy of approximate techniques for including electron correlation have been calibrated by comparison with full configuration-interaction calculations. Examples of the application of ab initio calculations to vibrational, rotational, and electronic spectroscopy are given. It is shown that the state-averaged, complete active space self-consistent field, multireference configuration-interaction procedure provides a good approach for treating several electronic states accurately in a common molecular orbital basis.

  7. AB 1725: A Comprehensive Analysis.

    ERIC Educational Resources Information Center

    California Community Colleges, Sacramento. Board of Governors.

    A summary and analysis is provided of California Assembly Bill (AB) 1725, a reform bill that provides new direction and support for the state's community colleges. The analysis addresses each of the eight sections of the bill: (1) mission, highlighting reforms related to mission statements, transfer core curriculum, remedial limits, articulation…

  8. Properties of AdeABC and AdeIJK Efflux Systems of Acinetobacter baumannii Compared with Those of the AcrAB-TolC System of Escherichia coli

    PubMed Central

    Sugawara, Etsuko

    2014-01-01

    Acinetobacter baumannii contains RND-family efflux systems AdeABC and AdeIJK, which pump out a wide range of antimicrobial compounds, as judged from the MIC changes occurring upon deletion of the responsible genes. However, these studies may miss changes because of the high backgrounds generated by the remaining pumps and by β-lactamases, and it is unclear how the activities of these pumps compare quantitatively with those of the well-studied AcrAB-TolC system of Escherichia coli. We expressed adeABC and adeIJK of A. baumannii, as well as E. coli acrAB, in an E. coli host from which acrAB was deleted. The A. baumannii pumps were functional in E. coli, and the MIC changes that were observed largely confirmed the substrate range already reported, with important differences. Thus, the AdeABC system pumped out all β-lactams, an activity that was often missed in deletion studies. When the expression level of the pump genes was adjusted to a similar level for a comparison with AcrAB-TolC, we found that both A. baumannii efflux systems pumped out a wide range of compounds, but AdeABC was less effective than AcrAB-TolC in the extrusion of lipophilic β-lactams, novobiocin, and ethidium bromide, although it was more effective at tetracycline efflux. AdeIJK was remarkably more effective than a similar level of AcrAB-TolC in the efflux of β-lactams, novobiocin, and ethidium bromide, although it was less so in the efflux of erythromycin. These results thus allow us to compare these efflux systems on a quantitative basis, if we can assume that the heterologous systems are fully functional in the E. coli host. PMID:25246403

  9. AB193. Rectourethral fistula

    PubMed Central

    Jiang, Hai

    2016-01-01

    Objective To investigate the treatment of rectourethral fistula. Methods Eleven cases of male patients with rectourethral fistula were treated in our department from 2011 to 2015. Age 16–66 years old. Causes: three cases of patients with congenital closed anus, four cases of traumatic pelvic fracture with urethral distraction and rectum injury, four cases after radical prostatectomy. The size of the fistula was 0.5–1.5 cm. In addition to the leakage of urine in the large fistula, urine mixed with stool samples. Three patients with congenital closed anal postoperative patients with posterior or anterior median sagittal approach for resection of the fistula, hierarchical closed urethral and rectal wall defect, at least three layer (between the urethral and rectal suture layer), indwelling catheter for 3–4 weeks, no cystostomy. Sigmoid colostomy underwent prior to the surgery. Of which six cases were repaired by perineal approach, one case by abdominal perineal approach, one case by abdominal repair. According to size of fistula and the surrounding scar decide whether or not to adopt tissue interposition, this group of five cases with gracilis muscle flap, one case with bulbocavernosus muscle flap interposed between the rectum and urethra; one case repaired by sigmoid colon pull-through procedure. Post-operation indwelling catheterization for 3–4 weeks with cystostomy. Results A total of 10 patients were successful, and no leakage of urine was found after removal of the catheter. One patient improved, occasionally a small amount of drops of urine voiding from anus. Reoperation was successful after 6 months. Recovered enteric continuity 3–6 months post-operation. Conclusions The median sagittal approach provide good exposure for the repair of congenital rectourethral fistula; perineal approach is a good choice for patients caused by trauma or surgery; complete resection of scar around the fistula, tension-free anastomosis, tissue interposition and sigmoid

  10. Representative Amino Acid Side-Chain Interactions in Protein-DNA Complexes: A Comparison of Highly Accurate Correlated Ab Initio Quantum Mechanical Calculations and Efficient Approaches for Applications to Large Systems.

    PubMed

    Hostaš, Jiří; Jakubec, Dávid; Laskowski, Roman A; Gnanasekaran, Ramachandran; Řezáč, Jan; Vondrášek, Jiří; Hobza, Pavel

    2015-09-01

    Representative pairs of amino acid side chains and nucleic acid bases extracted from available high-quality structures of protein-DNA complexes were analyzed using a range of computational methods. CCSD(T)/CBS interaction energies were calculated for the chosen 272 pairs. These reference interaction energies were used to test the MP2.5/CBS, MP2.X/CBS, MP2-F12, DFT-D3, PM6, and Amber force field methods. Method MP2.5 provided excellent agreement with reference data (root-mean-square error (RMSE) of 0.11 kcal/mol), which is more than 1 order of magnitude faster than the CCSD(T) method. When MP2-F12 and MP2.5 were combined, the results were within reasonable accuracy (0.20 kcal/mol), with a computational savings of almost 2 orders of magnitude. Therefore, this method is a promising tool for accurate calculations of interaction energies in protein-DNA motifs of up to ∼100 atoms, for which CCSD(T)/CBS benchmark calculations are not feasible. B3-LYP-D3 calculated with def2-TZVPP and def2-QZVP basis sets yielded sufficiently good results with a reasonably small RMSE. This method provided better results for neutral systems, whereas positively charged species exhibited the worst agreement with the benchmark data. The Amber force field yielded unbalanced results-performing well for systems containing nonpolar amino acids but severely underestimating interaction energies for charged complexes. The semiempirical PM6 method with corrections for hydrogen bonding and dispersion energy (PM6-D3H4) exhibited considerably smaller error than the Amber force field, which makes it an effective tool for modeling extended protein-ligand complexes (of up to 10,000 atoms). PMID:26575904

  11. An improved ab initio structure for fluorine peroxide (FOOF)

    NASA Astrophysics Data System (ADS)

    Mack, Hans-Georg; Oberhammer, Heinz

    1988-03-01

    Ab initio calculations with the 6-31G* and Dunning (9s5p/4s2p) basis sets augmented with p and d functions at various levels of theory (RHF, MP2, MP3, and MP4) were carried out on F 2O 2. The best result was obtained at the MP2 level with the Dunning basis plus one set of d functions on fluorine and two sets of d functions on oxygen. These calculations reproduce the experimental bond lengths to within 0.01 Å and the angles to within the experimental uncertainties.

  12. Ab initio study of the molecular structure and vibrational spectrum of nitric acid and its protonated forms

    NASA Technical Reports Server (NTRS)

    Lee, Timothy J.; Rice, Julia E.

    1992-01-01

    The equilibrium structures, harmonic vibrational frequencies, IR intensities, and relative energetics of HNO3 and its protonated form H2NO3+ were investigated using double-zeta plus polarization and triple-zeta plus polarization basis sets in conjunction with high-level ab initio methods. The latter include second-order Moller-Plesset perturbation theory, the single and double excitation coupled cluster (CCSD) methods, a perturbational estimate of the effects of connected triple excitations (CCSD(T)), and the self-consistent field. To determine accurate energy differences CCSD(T) energies were computed using large atomic natural orbital basis sets. Four different isomers of H2NO3+ were considered. The lowest energy form of protonated nitric acid was found to correspond to a complex between H2O and NO2+, which is consistent with earlier theoretical and experimental studies.

  13. Extension of a regularizing algorithm for the determination of equilibrium geometry and force field of free molecules from joint use of electron diffraction, molecular spectroscopy and ab initio data on systems with large-amplitude oscillatory motion

    NASA Astrophysics Data System (ADS)

    Kochikov, I. V.; Tarasov, Yu. I.; Spiridonov, V. P.; Kuramshina, G. M.; Yagola, A. G.; Saakjan, A. S.; Popik, M. V.; Samdal, S.

    1999-08-01

    The previously developed integrated algorithm for the joint treatment of gas-phase electron diffraction and vibrational spectroscopic data is extended to include systems with large-amplitude oscillatory motion. In addition, the treatment is augmented by the inclusion of microwave rotational constants. As in the previous work, the analysis of data from experimental sources is guided by quantum mechanical molecular geometry and force field optimization results. The computed force field matrix can be corrected empirically with the aid of suitable scale factors. Centrifugal distortion corrections to interatomic distances are included. The standard deviations of the parameters determined and the correlation coefficients can now be estimated. The principal design of the developed computer program is outlined, and some methodological problems associated with diffraction analysis of molecules with large-amplitude motion are discussed. To provide an example of a problem susceptible to attack by the present method an account is made of the re-analysis of diffraction data for 4-fluorobenzaldehyde collected earlier on the Balzers apparatus in Oslo.

  14. An ab initioCalculation of the Dipole Moment Surfaces and the Vibrational Transition Moments of the H 2Te Molecule

    NASA Astrophysics Data System (ADS)

    Kozin, Igor N.; Jensen, Per; Li, Yan; Buenker, Robert J.; Hirsch, Gerhard; Klee, Stefan

    1997-01-01

    The present work reports an ab initioMRD-CI calculation of the dipole moment surfaces for the electronic ground state of the H 2Te molecule. Using the ab initioresults, we calculate the vibrational transition moments, and we simulate the far-infrared spectrum of H 2Te by means of the MORBID program system. We obtain the equilibrium value of the dipole moment from the ab initiocalculation as 0.377 Debye based on our initial theoretical treatment which was employed over a wide range of molecular geometries. However, the use of an improved AO basis at the equilibrium geometry of H 2Te lowers this result to 0.298 Debye. The comparison of our simulated far-infrared spectrum with the experimental spectrum suggests that this value is too large, and that the correct value is certainly larger than 0.19 Debye and very probably smaller than 0.26 Debye. From the ab initiodata, we predict many vibrational transition moments for H 2Te, D 2Te, and HDTe. We hope that these results will be of assistance in the interpretation of the rotation-vibration spectrum of these molecules.

  15. Antigen identification and characterization of lung cancer specific monoclonal antibodies produced by mAb proteomics.

    PubMed

    Wang, Dongdong; Hincapie, Marina; Guergova-Kuras, Mariana; Kadas, Janos; Takacs, Laszlo; Karger, Barry L

    2010-04-01

    A mass spectrometric (MS)-based strategy for antigen (Ag) identification and characterization of globally produced monoclonal antibodies (mAbs) is described. Mice were immunized with a mixture of native glycoproteins, isolated from the pooled plasma of patients with nonsmall cell lung cancer (NSCLC), to generate a library of IgG-secreting hybridomas. Prior to immunization, the pooled NSCLC plasma was subjected to 3-sequential steps of affinity fractionation, including high abundant plasma protein depletion, glycoprotein enrichment, and polyclonal antibody affinity chromatography normalization. In this paper, to demonstrate the high quality of the globally produced mAbs, we selected 3 mAbs of high differentiating power against a matched, pooled normal plasma sample. After production of large quantities of the mAbs from ascites fluids, Ag identification was achieved by immunoaffinity purification, SDS-PAGE, Western blotting, and MS analysis of in-gel digest products. One antigen was found to be complement factor H, and the other two were mapped to different subunits of haptoglobin (Hpt). The 2 Hpt mAbs were characterized in detail to assess the quality of the mAbs produced by the global strategy. The affinity of one of the mAbs to the Hpt native tetramer form was found to have a K(D) of roughly 10(-9) M and to be 2 orders of magnitude lower than the reduced form, demonstrating the power of the mAb proteomics technology in generating mAbs to the natural form of the proteins in blood. The binding of this mAb to the beta-chain of haptoglobin was also dependent on glycosylation on this chain. The characterization of mAbs in this work reveals that the global mAb proteomics process can generate high-quality lung cancer specific mAbs capable of recognizing proteins in their native state. PMID:20146545

  16. Antigen Identification and Characterization of Lung Cancer Specific Monoclonal Antibodies Produced by mAb Proteomics

    PubMed Central

    Wang, Dongdong; Hincapie, Marina; Guergova-Kuras, Mariana; Kadas, Janos; Takacs, Laszlo; Karger, Barry L.

    2010-01-01

    A mass spectrometric (MS)-based strategy for antigen (Ag) identification and characterization of globally produced monoclonal antibodies (mAbs) is described. Mice were immunized with a mixture of native glycoproteins, isolated from the pooled plasma of patients with non-small cell lung cancer (NSCLC), to generate a library of IgG-secreting hybridomas. Prior to immunization, the pooled NSCLC plasma was subjected to 3-sequential steps of affinity fractionation, including high abundant plasma protein depletion, glycoprotein enrichment and polyclonal antibody affinity chromatography normalization. In this paper, in order to demonstrate the high quality of the globally produced mAbs, we selected 3 mAbs of high differentiating power against a matched, pooled normal plasma sample. After production of large quantities of the mAbs from ascites fluids, Ag identification was achieved by immunoaffinity purification, SDS-PAGE, Western blotting and MS analysis of in-gel digest products. One antigen was found to be complement factor H, and the other two were mapped to different subunits of haptoglobin (Hpt). The 2 Hpt mAbs were characterized in detail in order to assess the quality of the mAbs produced by the global strategy. The affinity of one of the mAbs to the Hpt native tetramer form was found to have a KD of roughly 10−9 M and to be 2 orders of magnitude lower than the reduced form, demonstrating the power of the mAb proteomics technology in generating mAbs to the natural form of the proteins in blood. The binding of this mAb to the β-chain of haptoglobin was also dependent on glycosylation on this chain. The characterization of mAbs in this work reveals that the global mAb proteomics process can generate high-quality lung cancer specific mAbs capable of recognizing proteins in their native state. PMID:20146545

  17. An operational multispectral scanner for bathymetric surveys - The ABS NORDA scanner

    NASA Technical Reports Server (NTRS)

    Haimbach, Stephen P.; Joy, Richard T.; Hickman, G. Daniel

    1987-01-01

    The Naval Ocean Research and Development Activity (NORDA) is developing the Airborne Bathymetric Survey (ABS) system, which will take shallow water depth soundings from a Navy P-3 aircraft. The system combines active and passive sensors to obtain optical measurements of water depth. The ABS NORDA Scanner is the systems passive multispectral scanner whose design goal is to provide 100 percent coverage of the seafloor, to depths of 20 m in average coastal waters. The ABS NORDA Scanner hardware and operational environment is discussed in detail. The optical model providing the basis for depth extraction is reviewed and the proposed data processing routine discussed.

  18. 7-(2-Thienyl)-7-Deazaadenosine (AB61), a New Potent Nucleoside Cytostatic with a Complex Mode of Action.

    PubMed

    Perlíková, Pavla; Rylová, Gabriela; Nauš, Petr; Elbert, Tomáš; Tloušťová, Eva; Bourderioux, Aurelie; Slavětínská, Lenka Poštová; Motyka, Kamil; Doležal, Dalibor; Znojek, Pawel; Nová, Alice; Harvanová, Monika; Džubák, Petr; Šiller, Michal; Hlaváč, Jan; Hajdúch, Marián; Hocek, Michal

    2016-05-01

    7-(2-Thienyl)-7-deazaadenosine (AB61) showed nanomolar cytotoxic activities against various cancer cell lines but only mild (micromolar) activities against normal fibroblasts. The selectivity of AB61 was found to be due to inefficient phosphorylation of AB61 in normal fibroblasts. The phosphorylation of AB61 in the leukemic CCRF-CEM cell line proceeds well and it was shown that AB61 is incorporated into both DNA and RNA, preferentially as a ribonucleotide. It was further confirmed that a triphosphate of AB61 is a substrate for both RNA and DNA polymerases in enzymatic assays. Gene expression analysis suggests that AB61 affects DNA damage pathways and protein translation/folding machinery. Indeed, formation of large 53BP1 foci was observed in nuclei of AB61-treated U2OS-GFP-53BP1 cells indicating DNA damage. Random incorporation of AB61 into RNA blocked its translation in an in vitro assay and reduction of reporter protein expression was also observed in mice after 4-hour treatment with AB61. AB61 also significantly reduced tumor volume in mice bearing SK-OV-3, BT-549, and HT-29 xenografts. The results indicate that AB61 is a promising compound with unique mechanism of action and deserves further development as an anticancer agent. Mol Cancer Ther; 15(5); 922-37. ©2016 AACR. PMID:26819331

  19. What Are the Ground State Structures of C20 and C24? An Explicitly Correlated Ab Initio Approach.

    PubMed

    Manna, Debashree; Martin, Jan M L

    2016-01-14

    A new benchmark study has been performed for six isomers of C20 and four isomers of C24 using explicitly correlated methods, together with coupled cluster theory with large basis sets and DFT with advanced functionals. The relative energy trends obtained are extremely sensitive to the methods used. Combining our best CCSD(T)-MP2 difference with our best MP2 basis set limit, the dehydrocorannulene bowl is found to be the most stable for C20, followed by the cage at about 8 kcal/mol, and the ring at about 46 kcal/mol. For C24, the D3d cage is found to be the most stable isomer, followed at only a few kilocalories per mole by dehydrocoronene, and at larger separations by then octahedral cage and the ring, respectively. This makes C24 the smallest classical fullerene. The estimated residual basis set error of the estimated CCSD(T) basis set limit is conservatively expected to be ±1 kcal/mol. In general, DFT exhibits large errors for relative energies with RMSD values in the 8-34 kcal/mol range. However, among the DFT functionals, the DSD-PBEP86-D3BJ double hybrid comes close to our best ab initio results, while the ωB97X-V range-separated hybrid is in semiquantitative agreement. PMID:26654916

  20. Applications of Basis Light-Front Quantization to QED

    NASA Astrophysics Data System (ADS)

    Vary, James P.; Zhao, Xingbo; Ilderton, Anton; Honkanen, Heli; Maris, Pieter; Brodsky, Stanley J.

    2014-06-01

    Hamiltonian light-front quantum field theory provides a framework for calculating both static and dynamic properties of strongly interacting relativistic systems. Invariant masses, correlated parton amplitudes and time-dependent scattering amplitudes, possibly with strong external time-dependent fields, represent a few of the important applications. By choosing the light-front gauge and adopting an orthonormal basis function representation, we obtain a large, sparse, Hamiltonian matrix eigenvalue problem for mass eigenstates that we solve by adapting ab initio no-core methods of nuclear many-body theory. In the continuum limit, the infinite matrix limit, we recover full covariance. Guided by the symmetries of light-front quantized theory, we adopt a two-dimensional harmonic oscillator basis for transverse modes that corresponds with eigensolutions of the soft-wall anti-de Sitter/quantum chromodynamics (AdS/QCD) model obtained from light-front holography. We outline our approach and present results for non-linear Compton scattering, evaluated non-perturbatively, where a strong and time-dependent laser field accelerates the electron and produces states of higher invariant mass i.e. final states with photon emission.

  1. Guiding ab initio calculations by alchemical derivatives.

    PubMed

    to Baben, M; Achenbach, J O; von Lilienfeld, O A

    2016-03-14

    We assess the concept of alchemical transformations for predicting how a further and not-tested change in composition would change materials properties. This might help to guide ab initio calculations through multidimensional property-composition spaces. Equilibrium volumes, bulk moduli, and relative lattice stability of fcc and bcc 4d transition metals Zr, Nb, Mo, Tc, Ru, Rh, Pd, and Ag are calculated using density functional theory. Alchemical derivatives predict qualitative trends in lattice stability while equilibrium volumes and bulk moduli are predicted with less than 9% and 28% deviation, respectively. Predicted changes in equilibrium volume and bulk moduli for binary and ternary mixtures of Rh-Pd-Ag are in qualitative agreement even for predicted bulk modulus changes as large as +100% or -50%. Based on these results, it is suggested that alchemical transformations could be meaningful for enhanced sampling in the context of virtual high-throughput materials screening projects. PMID:26979677

  2. Guiding ab initio calculations by alchemical derivatives

    NASA Astrophysics Data System (ADS)

    to Baben, M.; Achenbach, J. O.; von Lilienfeld, O. A.

    2016-03-01

    We assess the concept of alchemical transformations for predicting how a further and not-tested change in composition would change materials properties. This might help to guide ab initio calculations through multidimensional property-composition spaces. Equilibrium volumes, bulk moduli, and relative lattice stability of fcc and bcc 4d transition metals Zr, Nb, Mo, Tc, Ru, Rh, Pd, and Ag are calculated using density functional theory. Alchemical derivatives predict qualitative trends in lattice stability while equilibrium volumes and bulk moduli are predicted with less than 9% and 28% deviation, respectively. Predicted changes in equilibrium volume and bulk moduli for binary and ternary mixtures of Rh-Pd-Ag are in qualitative agreement even for predicted bulk modulus changes as large as +100% or -50%. Based on these results, it is suggested that alchemical transformations could be meaningful for enhanced sampling in the context of virtual high-throughput materials screening projects.

  3. Basis selection in LOBPCG

    NASA Astrophysics Data System (ADS)

    Hetmaniuk, U.; Lehoucq, R.

    2006-10-01

    The purpose of our paper is to discuss basis selection for Knyazev's locally optimal block preconditioned conjugate gradient (LOBPCG) method. An inappropriate choice of basis can lead to ill-conditioned Gram matrices in the Rayleigh-Ritz analysis that can delay convergence or produce inaccurate eigenpairs. We demonstrate that the choice of basis is not merely related to computing in finite precision arithmetic. We propose a representation that maintains orthogonality of the basis vectors and so has excellent numerical properties.

  4. Explicitly correlated atomic orbital basis second order Møller-Plesset theory

    NASA Astrophysics Data System (ADS)

    Hollman, David S.; Wilke, Jeremiah J.; Schaefer, Henry F.

    2013-02-01

    The scope of problems treatable by ab initio wavefunction methods has expanded greatly through the application of local approximations. In particular, atomic orbital (AO) based wavefunction methods have emerged as powerful techniques for exploiting sparsity and have been applied to biomolecules as large as 1707 atoms [S. A. Maurer, D. S. Lambrecht, D. Flaig, and C. Ochsenfeld, J. Chem. Phys. 136, 144107 (2012)], 10.1063/1.3693908. Correlated wavefunction methods, however, converge notoriously slowly to the basis set limit and, excepting the use of large basis sets, will suffer from a severe basis set incompleteness error (BSIE). The use of larger basis sets is prohibitively expensive for AO basis methods since, for example, second-order Møller-Plesset perturbation theory (MP2) scales linearly with the number of atoms, but still scales as O(N^5) in the number of functions per atom. Explicitly correlated F12 methods have been shown to drastically reduce BSIE for even modestly sized basis sets. In this work, we therefore explore an atomic orbital based formulation of explicitly correlated MP2-F12 theory. We present working equations for the new method, which produce results identical to the widely used molecular orbital (MO) version of MP2-F12 without resorting to a delocalized MO basis. We conclude with a discussion of several possible approaches to a priori screening of contraction terms in our method and the prospects for a linear scaling implementation of AO-MP2-F12. The discussion includes concrete examples involving noble gas dimers and linear alkane chains.

  5. Phonocatalysis. An ab initio simulation experiment

    NASA Astrophysics Data System (ADS)

    Kim, Kwangnam; Kaviany, Massoud

    2016-06-01

    Using simulations, we postulate and show that heterocatalysis on large-bandgap semiconductors can be controlled by substrate phonons, i.e., phonocatalysis. With ab initio calculations, including molecular dynamic simulations, the chemisorbed dissociation of XeF6 on h-BN surface leads to formation of XeF4 and two surface F/h-BN bonds. The reaction pathway and energies are evaluated, and the sorption and reaction emitted/absorbed phonons are identified through spectral analysis of the surface atomic motion. Due to large bandgap, the atomic vibration (phonon) energy transfer channels dominate and among them is the match between the F/h-BN covalent bond stretching and the optical phonons. We show that the chemisorbed dissociation (the pathway activation ascent) requires absorption of large-energy optical phonons. Then using progressively heavier isotopes of B and N atoms, we show that limiting these high-energy optical phonons inhibits the chemisorbed dissociation, i.e., controllable phonocatalysis.

  6. The structure and binding energies of the van der Waals complexes of Ar and N2 with phenol and its cation, studied by high level ab initio and density functional theory calculations

    NASA Astrophysics Data System (ADS)

    Vincent, Mark A.; Hillier, Ian H.; Morgado, Claudio A.; Burton, Neil A.; Shan, Xiao

    2008-01-01

    We have investigated, using both ab initio and density functional theory methods, the minimum energy structures and corresponding binding energies of the van der Waals complexes between phenol and argon or the nitrogen molecule, and the corresponding complexes involving the phenol cation. Structures were obtained at the MP2 level using a large basis, and the corresponding energies were corrected for basis set superposition error (BSSE), higher order electron correlation effects, and for basis set size. The structures of the global minima were further refined for the effects of BSSE and the corresponding binding energies were evaluated. For each neutral species, we find only a single true minimum, π bonded for argon and OH bonded for nitrogen. For both cationic species, we find that the OH-bonded complex is preferred over other minima which we have identified as having Ar or N2 between exogeneous atoms. The ab initio calculations are generally in excellent agreement with experimental binding energies and rotational constants. We find that the B3LYP functional is particularly poor at describing these complexes, while a density functional theory (DFT) method with an empirical correction for dispersive interactions (DFT-D) is very successful, as are some of the new functionals proposed by Zhao and Truhlar [J. Phys. Chem. A 109, 5656 (2005); J. Chem. Theory Comput. 2, 1009 (2006); Phys. Chem. Chem. Phys. 7, 2701 (2005); J. Phys. Chem. A 108, 6908 (2004)]. Both the ab initio and DFT-D methods accurately predict the intermolecular vibrational modes.

  7. Five ab initio potential energy and dipole moment surfaces for hydrated NaCl and NaF. I. Two-body interactions

    NASA Astrophysics Data System (ADS)

    Wang, Yimin; Bowman, Joel M.; Kamarchik, Eugene

    2016-03-01

    We report full-dimensional, ab initio-based potentials and dipole moment surfaces for NaCl, NaF, Na+H2O, F-H2O, and Cl-H2O. The NaCl and NaF potentials are diabatic ones that dissociate to ions. These are obtained using spline fits to CCSD(T)/aug-cc-pV5Z energies. In addition, non-linear least square fits using the Born-Mayer-Huggins potential are presented, providing accurate parameters based strictly on the current ab initio energies. The long-range behavior of the NaCl and NaF potentials is shown to go, as expected, accurately to the point-charge Coulomb interaction. The three ion-H2O potentials are permutationally invariant fits to roughly 20 000 coupled cluster CCSD(T) energies (awCVTZ basis for Na+ and aVTZ basis for Cl- and F-), over a large range of distances and H2O intramolecular configurations. These potentials are switched accurately in the long range to the analytical ion-dipole interactions, to improve computational efficiency. Dipole moment surfaces are fits to MP2 data; for the ion-ion cases, these are well described in the intermediate- and long-range by the simple point-charge expression. The performance of these new fits is examined by direct comparison to additional ab initio energies and dipole moments along various cuts. Equilibrium structures, harmonic frequencies, and electronic dissociation energies are also reported and compared to direct ab initio results. These indicate the high fidelity of the new PESs.

  8. Five ab initio potential energy and dipole moment surfaces for hydrated NaCl and NaF. I. Two-body interactions.

    PubMed

    Wang, Yimin; Bowman, Joel M; Kamarchik, Eugene

    2016-03-21

    We report full-dimensional, ab initio-based potentials and dipole moment surfaces for NaCl, NaF, Na(+)H2O, F(-)H2O, and Cl(-)H2O. The NaCl and NaF potentials are diabatic ones that dissociate to ions. These are obtained using spline fits to CCSD(T)/aug-cc-pV5Z energies. In addition, non-linear least square fits using the Born-Mayer-Huggins potential are presented, providing accurate parameters based strictly on the current ab initio energies. The long-range behavior of the NaCl and NaF potentials is shown to go, as expected, accurately to the point-charge Coulomb interaction. The three ion-H2O potentials are permutationally invariant fits to roughly 20,000 coupled cluster CCSD(T) energies (awCVTZ basis for Na(+) and aVTZ basis for Cl(-) and F(-)), over a large range of distances and H2O intramolecular configurations. These potentials are switched accurately in the long range to the analytical ion-dipole interactions, to improve computational efficiency. Dipole moment surfaces are fits to MP2 data; for the ion-ion cases, these are well described in the intermediate- and long-range by the simple point-charge expression. The performance of these new fits is examined by direct comparison to additional ab initio energies and dipole moments along various cuts. Equilibrium structures, harmonic frequencies, and electronic dissociation energies are also reported and compared to direct ab initio results. These indicate the high fidelity of the new PESs. PMID:27004880

  9. AB INITIO AND CALPHAD THERMODYNAMICS OF MATERIALS

    SciTech Connect

    Turchi, P A

    2004-04-14

    Ab initio electronic structure methods can supplement CALPHAD in two major ways for subsequent applications to stability in complex alloys. The first one is rather immediate and concerns the direct input of ab initio energetics in CALPHAD databases. The other way, more involved, is the assessment of ab initio thermodynamics {acute a} la CALPHAD. It will be shown how these results can be used within CALPHAD to predict the equilibrium properties of multi-component alloys.

  10. Magnetic topology and prominence patterns on AB Doradus

    NASA Astrophysics Data System (ADS)

    Donati, J.-F.; Collier Cameron, A.; Hussain, G. A. J.; Semel, M.

    1999-01-01

    We report new Zeeman--Doppler imaging observations of the rapidly rotating young K0 dwarf AB Doradus, obtained with the Anglo-Australian Telescope in 1996 December. From such observations, simultaneous brightness and magnetic images of the stellar photosphere of AB Dor were reconstructed at three different epochs. The magnetic topology of AB Dor is found to be very complex, with at least 12 different radial field regions of opposite polarities located all around the star. Significant azimuthal field fluxes are also detected in the form of one negative polarity region close to the equator, a series of positive polarity patches at intermediate latitudes and an almost complete ring of negative polarity encircling the rotational pole at high latitudes. In particular, the azimuthal polarities we reconstruct are in very good agreement with those obtained by Donati & Cameron, confirming that this field component is directly related to the dynamo-generated large-scale toroidal magnetic structure. The triple polarity latitudinal pattern observed for this structure in the upper hemisphere of AB Dor indicates that the degree of the underlying large-scale poloidal structure in an axisymmetric spherical harmonics expansion is equal to or greater than five. It also strengthens the idea that the dynamo processes operating in AB Dor feature a non-solar component distributed throughout the convective zone. From the subtle distortion of successive brightness images, we can also confirm the surface differential rotation first measured on this star by Donati & Cameron in both sense and magnitude, with a pole rotating more slowly than the equator by about one part in 220. Finally, the rotation periods we measure for four prominences (from the recurrence rate of their spectral signatures in Balmer lines) confirm the presumption that such clouds are anchored at intermediate to high latitudes. The intrinsic variability of these prominences is not associated with any abrupt changes in the

  11. A procedure for computing accurate ab initio quartic force fields: Application to HO2+ and H2O

    NASA Astrophysics Data System (ADS)

    Huang, Xinchuan; Lee, Timothy J.

    2008-07-01

    A procedure for the calculation of molecular quartic force fields (QFFs) is proposed and investigated. The goal is to generate highly accurate ab initio QFFs that include many of the so-called ``small'' effects that are necessary to achieve high accuracy. The small effects investigated in the present study include correlation of the core electrons (core correlation), extrapolation to the one-particle basis set limit, correction for scalar relativistic contributions, correction for higher-order correlation effects, and inclusion of diffuse functions in the one-particle basis set. The procedure is flexible enough to allow for some effects to be computed directly, while others may be added as corrections. A single grid of points is used and is centered about an initial reference geometry that is designed to be as close as possible to the final ab initio equilibrium structure (with all effects included). It is shown that the least-squares fit of the QFF is not compromised by the added corrections, and the balance between elimination of contamination from higher-order force constants while retaining energy differences large enough to yield meaningful quartic force constants is essentially unchanged from the standard procedures we have used for many years. The initial QFF determined from the least-squares fit is transformed to the exact minimum in order to eliminate gradient terms and allow for the use of second-order perturbation theory for evaluation of spectroscopic constants. It is shown that this step has essentially no effect on the quality of the QFF largely because the initial reference structure is, by design, very close to the final ab initio equilibrium structure. The procedure is used to compute an accurate, purely ab initio QFF for the H2O molecule, which is used as a benchmark test case. The procedure is then applied to the ground and first excited electronic states of the HO2+ molecular cation. Fundamental vibrational frequencies and spectroscopic

  12. The Structural Basis of Antibody-Antigen Recognition

    PubMed Central

    Sela-Culang, Inbal; Kunik, Vered; Ofran, Yanay

    2013-01-01

    The function of antibodies (Abs) involves specific binding to antigens (Ags) and activation of other components of the immune system to fight pathogens. The six hypervariable loops within the variable domains of Abs, commonly termed complementarity determining regions (CDRs), are widely assumed to be responsible for Ag recognition, while the constant domains are believed to mediate effector activation. Recent studies and analyses of the growing number of available Ab structures, indicate that this clear functional separation between the two regions may be an oversimplification. Some positions within the CDRs have been shown to never participate in Ag binding and some off-CDRs residues often contribute critically to the interaction with the Ag. Moreover, there is now growing evidence for non-local and even allosteric effects in Ab-Ag interaction in which Ag binding affects the constant region and vice versa. This review summarizes and discusses the structural basis of Ag recognition, elaborating on the contribution of different structural determinants of the Ab to Ag binding and recognition. We discuss the CDRs, the different approaches for their identification and their relationship to the Ag interface. We also review what is currently known about the contribution of non-CDRs regions to Ag recognition, namely the framework regions (FRs) and the constant domains. The suggested mechanisms by which these regions contribute to Ag binding are discussed. On the Ag side of the interaction, we discuss attempts to predict B-cell epitopes and the suggested idea to incorporate Ab information into B-cell epitope prediction schemes. Beyond improving the understanding of immunity, characterization of the functional role of different parts of the Ab molecule may help in Ab engineering, design of CDR-derived peptides, and epitope prediction. PMID:24115948

  13. Multiple time step integrators in ab initio molecular dynamics

    SciTech Connect

    Luehr, Nathan; Martínez, Todd J.; Markland, Thomas E.

    2014-02-28

    Multiple time-scale algorithms exploit the natural separation of time-scales in chemical systems to greatly accelerate the efficiency of molecular dynamics simulations. Although the utility of these methods in systems where the interactions are described by empirical potentials is now well established, their application to ab initio molecular dynamics calculations has been limited by difficulties associated with splitting the ab initio potential into fast and slowly varying components. Here we present two schemes that enable efficient time-scale separation in ab initio calculations: one based on fragment decomposition and the other on range separation of the Coulomb operator in the electronic Hamiltonian. We demonstrate for both water clusters and a solvated hydroxide ion that multiple time-scale molecular dynamics allows for outer time steps of 2.5 fs, which are as large as those obtained when such schemes are applied to empirical potentials, while still allowing for bonds to be broken and reformed throughout the dynamics. This permits computational speedups of up to 4.4x, compared to standard Born-Oppenheimer ab initio molecular dynamics with a 0.5 fs time step, while maintaining the same energy conservation and accuracy.

  14. Ab-Initio Shell Model with a Core

    SciTech Connect

    Lisetskiy, A F; Barrett, B R; Kruse, M; Navratil, P; Stetcu, I; Vary, J P

    2008-06-04

    We construct effective 2- and 3-body Hamiltonians for the p-shell by performing 12{h_bar}{Omega} ab initio no-core shell model (NCSM) calculations for A=6 and 7 nuclei and explicitly projecting the many-body Hamiltonians onto the 0{h_bar}{Omega} space. We then separate these effective Hamiltonians into 0-, 1- and 2-body contributions (also 3-body for A=7) and analyze the systematic behavior of these different parts as a function of the mass number A and size of the NCSM basis space. The role of effective 3- and higher-body interactions for A > 6 is investigated and discussed.

  15. Ab initio Study of He Stability in hcp-Ti

    SciTech Connect

    Dai, Yunya; Yang, Li; Peng, SM; Long, XG; Gao, Fei; Zu, Xiaotao T.

    2010-12-20

    The stability of He in hcp-Ti was studied using ab initio method based on density functional theory. The results indicate that a single He atom prefers to occupy the tetrahedral site rather than the octahedral site. The interaction of He defects with Ti atoms has been used to explain the relative stabilities of He point defects in hcp-Ti. The relative stability of He defects in hcp-Ti is useful for He clustering and bubble nucleation in metal tritides, which provides the basis for development of improved atomistic models.

  16. Ultrasound Biomicroscopy Comparison of Ab Interno and Ab Externo Intraocular Lens Scleral Fixation.

    PubMed

    Horiguchi, Lie; Garcia, Patricia Novita; Malavazzi, Gustavo Ricci; Allemann, Norma; Gomes, Rachel L R

    2016-01-01

    Purpose. To compare ab interno and ab externo scleral fixation of posterior chamber intraocular lenses (PCIOL) using ultrasound biomicroscopy (UBM). Methods. Randomized patients underwent ab externo or ab interno scleral fixation of a PCIOL. Ultrasound biomicroscopy was performed 3 to 6 months postoperatively, to determine PCIOL centration, IOL distance to the iris at 12, 3, 6, and 9 hours, and haptics placement in relation to the ciliary sulcus. Results. Fifteen patients were enrolled in the study. The ab externo technique was used in 7 eyes (46.6%) and the ab interno in 8 eyes (53.3%). In the ab externo technique, 14 haptics were located: 4 (28.57%) in the ciliary sulcus; 2 (14.28%) anterior to the sulcus; and 8 (57.14%) posterior to the sulcus, 6 in the ciliary body and 2 posterior to the ciliary body. In the ab interno group, 4 haptics (25.0%) were in the ciliary sulcus, 2 (12.50%) anterior to the sulcus, and 10 (75.0%) posterior to the sulcus, 4 in the ciliary body and 6 posterior to the ciliary body. Conclusions. Ab externo and ab interno scleral fixation techniques presented similar results in haptic placement. Ab externo technique presented higher vertical tilt when compared to the ab interno. PMID:27293878

  17. Ultrasound Biomicroscopy Comparison of Ab Interno and Ab Externo Intraocular Lens Scleral Fixation

    PubMed Central

    Horiguchi, Lie; Garcia, Patricia Novita; Malavazzi, Gustavo Ricci; Allemann, Norma

    2016-01-01

    Purpose. To compare ab interno and ab externo scleral fixation of posterior chamber intraocular lenses (PCIOL) using ultrasound biomicroscopy (UBM). Methods. Randomized patients underwent ab externo or ab interno scleral fixation of a PCIOL. Ultrasound biomicroscopy was performed 3 to 6 months postoperatively, to determine PCIOL centration, IOL distance to the iris at 12, 3, 6, and 9 hours, and haptics placement in relation to the ciliary sulcus. Results. Fifteen patients were enrolled in the study. The ab externo technique was used in 7 eyes (46.6%) and the ab interno in 8 eyes (53.3%). In the ab externo technique, 14 haptics were located: 4 (28.57%) in the ciliary sulcus; 2 (14.28%) anterior to the sulcus; and 8 (57.14%) posterior to the sulcus, 6 in the ciliary body and 2 posterior to the ciliary body. In the ab interno group, 4 haptics (25.0%) were in the ciliary sulcus, 2 (12.50%) anterior to the sulcus, and 10 (75.0%) posterior to the sulcus, 4 in the ciliary body and 6 posterior to the ciliary body. Conclusions. Ab externo and ab interno scleral fixation techniques presented similar results in haptic placement. Ab externo technique presented higher vertical tilt when compared to the ab interno. PMID:27293878

  18. The equilibrium geometry of the SC3H radical: an ab initio study

    NASA Astrophysics Data System (ADS)

    Flores, J. R.; Gómez, F. J.

    The SC3 H radical is known by experiment to have a linear equilibrium structure, but even rather high-level ab initio computations give a bent equilibrium geometry. A theoretical study of the SCCCH radical has been carried out in order to analyse the influence of several factors in the computed equilibrium structure. Quadratic configuration interaction QCISD(T) and restricted coupled cluster RCCSD(T) computations have been performed in combination with large basis sets. Spin-orbit effects have been taken into account through the Breit-Pauli Hamiltonian using multi-configuration SCF and configuration interaction wavefunctions. Our final results indicate that the equilibrium structure must be linear, in agreement with the experimental studies [MCCARTHY, M. C., VRTILEK, J. M., GOTTLIEB, C. A., WANG, W., and THADDEUS, P., 1994, Astrophys. J., 431, L127; HIRAHARA, Y., OHSHIMA, Y., and ENDO, Y., 1994, J. chem. Phys., 101, 7342]. Both spin-orbit and electron correlation effects appear to be of comparable importance, but an adequate computation of the correlation energy has been much more difficult and has ultimately required bas.is set extrapolations.

  19. An ab initio investigation of possible intermediates in the reaction of the hydroxyl and hydroperoxyl radicals

    NASA Technical Reports Server (NTRS)

    Jackels, C. F.

    1985-01-01

    Ab initio quantum chemical techniques are used to investigate covalently-bonded and hydrogen-bonded species that may be important intermediates in the reaction of hydroxyl and hydroperoxyl radicals. Stable structures of both types are identified. Basis sets of polarized double zeta quality and large scale configuration interaction wave functions are utilized. Based on electronic energies, the covalently bonded HOOOH species is 26.4 kcal/mol more stable than the OH and HO2 radicals. Similarly, the hydrogen bonded HO---HO2 species has an electronic energy 4.7 kcal/mol below that of the component radicals, after correction is made for the basis set superposition error. The hydrogen bonded form is planar, possesses one relatively normal hydrogen bond, and has the lowest energy 3A' and 1A' states that are essentially degenerate. The 1A" and 3A" excited states produced by rotation of the unpaired OH electron into the molecular plane are very slightly bound.

  20. Atomic Spectral Methods for Ab Initio Molecular Electronic Energy Surfaces: Transitioning From Small-Molecule to Biomolecular-Suitable Approaches.

    PubMed

    Mills, Jeffrey D; Ben-Nun, Michal; Rollin, Kyle; Bromley, Michael W J; Li, Jiabo; Hinde, Robert J; Winstead, Carl L; Sheehy, Jeffrey A; Boatz, Jerry A; Langhoff, Peter W

    2016-08-25

    Continuing attention has addressed incorportation of the electronically dynamical attributes of biomolecules in the largely static first-generation molecular-mechanical force fields commonly employed in molecular-dynamics simulations. We describe here a universal quantum-mechanical approach to calculations of the electronic energy surfaces of both small molecules and large aggregates on a common basis which can include such electronic attributes, and which also seems well-suited to adaptation in ab initio molecular-dynamics applications. In contrast to the more familiar orbital-product-based methodologies employed in traditional small-molecule computational quantum chemistry, the present approach is based on an "ex-post-facto" method in which Hamiltonian matrices are evaluated prior to wave function antisymmetrization, implemented here in the support of a Hilbert space of orthonormal products of many-electron atomic spectral eigenstates familiar from the van der Waals theory of long-range interactions. The general theory in its various forms incorporates the early semiempirical atoms- and diatomics-in-molecules approaches of Moffitt, Ellison, Tully, Kuntz, and others in a comprehensive mathematical setting, and generalizes the developments of Eisenschitz, London, Claverie, and others addressing electron permutation symmetry adaptation issues, completing these early attempts to treat van der Waals and chemical forces on a common basis. Exact expressions are obtained for molecular Hamiltonian matrices and for associated energy eigenvalues as sums of separate atomic and interaction-energy terms, similar in this respect to the forms of classical force fields. The latter representation is seen to also provide a long-missing general definition of the energies of individual atoms and of their interactions within molecules and matter free from subjective additional constraints. A computer code suite is described for calculations of the many-electron atomic eigenspectra and

  1. Hamiltonian light-front field theory in a basis function approach

    SciTech Connect

    Vary, J. P.; Honkanen, H.; Li Jun; Maris, P.; Brodsky, S. J.; Harindranath, A.; Sternberg, P.; Ng, E. G.; Yang, C.

    2010-03-15

    Hamiltonian light-front quantum field theory constitutes a framework for the nonperturbative solution of invariant masses and correlated parton amplitudes of self-bound systems. By choosing the light-front gauge and adopting a basis function representation, a large, sparse, Hamiltonian matrix for mass eigenstates of gauge theories is obtained that is solvable by adapting the ab initio no-core methods of nuclear many-body theory. Full covariance is recovered in the continuum limit, the infinite matrix limit. There is considerable freedom in the choice of the orthonormal and complete set of basis functions with convenience and convergence rates providing key considerations. Here we use a two-dimensional harmonic oscillator basis for transverse modes that corresponds with eigensolutions of the soft-wall anti-de Sitter/quantum chromodynamics (AdS/QCD) model obtained from light-front holography. We outline our approach and present illustrative features of some noninteracting systems in a cavity. We illustrate the first steps toward solving quantum electrodynamics (QED) by obtaining the mass eigenstates of an electron in a cavity in small basis spaces and discuss the computational challenges.

  2. Safety Basis Report

    SciTech Connect

    R.J. Garrett

    2002-01-14

    As part of the internal Integrated Safety Management Assessment verification process, it was determined that there was a lack of documentation that summarizes the safety basis of the current Yucca Mountain Project (YMP) site characterization activities. It was noted that a safety basis would make it possible to establish a technically justifiable graded approach to the implementation of the requirements identified in the Standards/Requirements Identification Document. The Standards/Requirements Identification Documents commit a facility to compliance with specific requirements and, together with the hazard baseline documentation, provide a technical basis for ensuring that the public and workers are protected. This Safety Basis Report has been developed to establish and document the safety basis of the current site characterization activities, establish and document the hazard baseline, and provide the technical basis for identifying structures, systems, and components (SSCs) that perform functions necessary to protect the public, the worker, and the environment from hazards unique to the YMP site characterization activities. This technical basis for identifying SSCs serves as a grading process for the implementation of programs such as Conduct of Operations (DOE Order 5480.19) and the Suspect/Counterfeit Items Program. In addition, this report provides a consolidated summary of the hazards analyses processes developed to support the design, construction, and operation of the YMP site characterization facilities and, therefore, provides a tool for evaluating the safety impacts of changes to the design and operation of the YMP site characterization activities.

  3. Comparative studies of the spectroscopy of CuCl2: DFT versus standard ab initio approaches.

    PubMed

    Ramírez-Solís, A; Poteau, R; Vela, A; Daudey, J P

    2005-04-22

    The X2Pi g-2Sigma g+, X2Pi g-2Delta g, X2Pi g-2Sigma u+, X2Pi g-2Pi u transitions on CuCl2 have been studied using several exchange-correlation functionals from the various types of density functional theory (DFT) approaches like local density approximation (LDA), generalized gradient approximation (GGA), hybrid and meta-GGA. The results are compared with the experience and with those coming from the most sophisticated nondynamic and dynamic electronic correlation treatments using the same relativistic effective core potentials and especially developed basis sets to study the electronic structure of the five lowest states and the corresponding vertical and adiabatic transition energies. The calculated transition energies for three of the hybrid functionals (B3LYP, B97-2, and PBE0) are in very good agreement with the benchmark ab initio results and experimental figures. All of the other functionals largely overestimate the X2Pi g-2Sigma g+ and X2Pi g-2Delta g transition energies, many of them even placing the 2Delta g ligand field state above the charge transfer 2Pi u and 2Sigma u+ states. The relative weight of the Hartree-Fock exchange in the definition of the functional used appears to play a key role in the accurate description of the LambdaSSigma density defined by the orientation of the 3d hole (sigma, pi, or delta) on Cu in the field of both chlorine atoms, but no simple connection of this weight with the quality of the spectra has been found. Mulliken charges and spin densities are carefully analyzed; a possible link between the extent of spin density on the metal for the X2Pi g state and the performance of the various functionals was observed, suggesting that those that lead to the largest values (close to 0.65) are the ones that best reproduce these four transitions. Most functionals lead to a remarkably low ionicity for the three ligand field states even for the best performing functionals, compared to the complete active space (SCF) (21, 14) ab initio

  4. Comparative studies of the spectroscopy of CuCl2: DFT versus standard ab initio approaches

    NASA Astrophysics Data System (ADS)

    Ramírez-Solís, A.; Poteau, R.; Vela, A.; Daudey, J. P.

    2005-04-01

    The XΠg2-Σg +2, XΠg2-Δg2, XΠg2-Σu +2, XΠg2-Πu2 transitions on CuCl2 have been studied using several exchange-correlation functionals from the various types of density functional theory (DFT) approaches like local density approximation (LDA), generalized gradient approximation (GGA), hybrid and meta-GGA. The results are compared with the experience and with those coming from the most sophisticated nondynamic and dynamic electronic correlation treatments using the same relativistic effective core potentials and especially developed basis sets to study the electronic structure of the five lowest states and the corresponding vertical and adiabatic transition energies. The calculated transition energies for three of the hybrid functionals (B3LYP, B97-2, and PBE0) are in very good agreement with the benchmark ab initio results and experimental figures. All of the other functionals largely overestimate the XΠg2-Σg +2 and XΠg2-Δg2 transition energies, many of them even placing the Δg2 ligand field state above the charge transfer Πu2 and Σu +2 states. The relative weight of the Hartree-Fock exchange in the definition of the functional used appears to play a key role in the accurate description of the ΛSΣ density defined by the orientation of the 3d hole (σ, π, or δ) on Cu in the field of both chlorine atoms, but no simple connection of this weight with the quality of the spectra has been found. Mulliken charges and spin densities are carefully analyzed; a possible link between the extent of spin density on the metal for the XΠg2 state and the performance of the various functionals was observed, suggesting that those that lead to the largest values (close to 0.65) are the ones that best reproduce these four transitions. Most functionals lead to a remarkably low ionicity for the three ligand field states even for the best performing functionals, compared to the complete active space (SCF) (21, 14) ab initio values. These findings show that not only large

  5. Ab initio study of magnesium and magnesium hydride nanoclusters and nanocrystals: examining optimal structures and compositions for efficient hydrogen storage.

    PubMed

    Koukaras, Emmanuel N; Zdetsis, Aristides D; Sigalas, Michael M

    2012-09-26

    On the basis of the attractive possibility of efficient hydrogen storage in light metal hydrides, we have examined a large variety of Mg(n)H(m) nanoclusters and (MgH(2))(n) nanocrystals (n = 2-216, m = 2-436) using high level coupled cluster, CCSD(T), ab initio methods, and judicially chosen density functional calculations of comparable quality and (near chemical) accuracy. Our calculated desorption energies as a function of size and percentage of hydrogen have pinpointed optimal regions of sizes and concentrations of hydrogen which are in full agreement with recent experimental findings. Furthermore, our results reproduce the experimental desorption energy of 75.5 kJ/mol for the infinite system with remarkable accuracy (76.5 ± 1.5 kJ/mol). PMID:22920498

  6. Ab initio investigations of A-site doping on the structure and electric polarization of HoMnO3

    NASA Astrophysics Data System (ADS)

    S, Sathya Sheela; C, Kanagaraj; Natesan, Baskaran

    2015-06-01

    We have investigated the effect of A-site doping on the structure and electric polarization of orthorhombic HoMnO3 using ab initio density functional theory calculations. We find that the substitution of rare earth ions, such as Lu, Y and La in place of Ho in orthorhombic HoMnO3 modifies the local structure around Mn ions drastically, and leads to the formation of two distinct Mn sites Mn(0) and Mn(1). As a result, large variance between Mn(0)O6 and Mn(1)O6 octahedral distortions arises. This variance in the octahedral distortions drives the disparate hopping of electrons between the eg orbitals enhancing the electronic polarization with increasing rare earth ion radius. The largest polarization of 7 µC/cm2 is obtained for La doped HoMnO3. This increase in polarization has been explained on the basis of radius mismatch induced local structural effects.

  7. Structural basis for the antibody neutralization of Herpes simplex virus

    SciTech Connect

    Lee, Cheng-Chung; Lin, Li-Ling; Chan, Woan-Eng; Ko, Tzu-Ping; Lai, Jiann-Shiun; Wang, Andrew H.-J.

    2013-10-01

    The gD–E317-Fab complex crystal revealed the conformational epitope of human mAb E317 on HSV gD, providing a molecular basis for understanding the viral neutralization mechanism. Glycoprotein D (gD) of Herpes simplex virus (HSV) binds to a host cell surface receptor, which is required to trigger membrane fusion for virion entry into the host cell. gD has become a validated anti-HSV target for therapeutic antibody development. The highly inhibitory human monoclonal antibody E317 (mAb E317) was previously raised against HSV gD for viral neutralization. To understand the structural basis of antibody neutralization, crystals of the gD ectodomain bound to the E317 Fab domain were obtained. The structure of the complex reveals that E317 interacts with gD mainly through the heavy chain, which covers a large area for epitope recognition on gD, with a flexible N-terminal and C-terminal conformation. The epitope core structure maps to the external surface of gD, corresponding to the binding sites of two receptors, herpesvirus entry mediator (HVEM) and nectin-1, which mediate HSV infection. E317 directly recognizes the gD–nectin-1 interface and occludes the HVEM contact site of gD to block its binding to either receptor. The binding of E317 to gD also prohibits the formation of the N-terminal hairpin of gD for HVEM recognition. The major E317-binding site on gD overlaps with either the nectin-1-binding residues or the neutralizing antigenic sites identified thus far (Tyr38, Asp215, Arg222 and Phe223). The epitopes of gD for E317 binding are highly conserved between two types of human herpesvirus (HSV-1 and HSV-2). This study enables the virus-neutralizing epitopes to be correlated with the receptor-binding regions. The results further strengthen the previously demonstrated therapeutic and diagnostic potential of the E317 antibody.

  8. Ab initio tight-binding Hamiltonian for transition metal dichalcogenides

    NASA Astrophysics Data System (ADS)

    Fang, Shiang; Kuate Defo, Rodrick; Shirodkar, Sharmila N.; Lieu, Simon; Tritsaris, Georgios A.; Kaxiras, Efthimios

    2015-11-01

    We present an accurate ab initio tight-binding Hamiltonian for the transition metal dichalcogenides, MoS2, MoSe2, WS2, WSe2, with a minimal basis (the d orbitals for the metal atoms and p orbitals for the chalcogen atoms) based on a transformation of the Kohn-Sham density functional theory Hamiltonian to a basis of maximally localized Wannier functions. The truncated tight-binding Hamiltonian, with only on-site, first, and partial second neighbor interactions, including spin-orbit coupling, provides a simple physical picture and the symmetry of the main band-structure features. Interlayer interactions between adjacent layers are modeled by transferable hopping terms between the chalcogen p orbitals. The full-range tight-binding Hamiltonian can be reduced to hybrid-orbital k .p effective Hamiltonians near the band extrema that capture important low-energy excitations. These ab initio Hamiltonians can serve as the starting point for applications to interacting many-body physics including optical transitions and Berry curvature of bands, of which we give some examples.

  9. A spinor boson AB chain

    NASA Astrophysics Data System (ADS)

    Cruz Reyes, Greis Julieth; Franco, Roberto; Silva Valencia, Jereson; Universidad Santo Tomas Collaboration; Universidad Nacional de Colombia Collaboration

    Recent research is focused on superlattices arising from optical lattices, which allow a tunable environment. Experimentally bosons present transitions from superfluid to Mott insulator by changing the energy offset in the unit cell [Nat. Commun. 5:5735 (2014)]. Many studies displayed that ground state of spinless boson systems on superlattices present superfluid, Mott insulator and an additional CDW phase created by the energy shift between the sites into the unit cell [Phys. Rev. A 83, 053621 (2011)]. The first confinement methods were magnetic traps, which freezes the spin; with optical lattices the grade of freedom of spin plays an important role. We consider bosons with spin S =1 on a superlattice made by two sites with energy offset per unit cell (AB chain). The Hamiltonian that describes the system is the Bose-Hubbard model with the superlattice potential (W) and the exchange interaction (V) parameters. This model supports CDW, Mott insulator and superfluid phases. For W near to U, with V =0, Mott phase disappears, but for V increasing, a new CDW appears due to the spin interaction, while the half-integer CDW decrease. These results are widely different from spinless boson, where the CDW phases are stables.

  10. Ab initio phonon limited transport

    NASA Astrophysics Data System (ADS)

    Verstraete, Matthieu

    We revisit the thermoelectric (TE) transport properties of two champion materials, PbTe and SnSe, using fully first principles methods. In both cases the performance of the material is due to subtle combinations of structural effects, scattering, and phase space reduction. In PbTe anharmonic effects are completely opposite to the predicted quasiharmonic evolution of phonon frequencies and to frequently (and incorrectly) cited extrapolations of experiments. This stabilizes the material at high T, but also tends to enhance its thermal conductivity, in a non linear manner, above 600 Kelvin. This explains why PbTe is in practice limited to room temperature applications. SnSe has recently been shown to be the most efficient TE material in bulk form. This is mainly due to a strongly enhanced carrier concentration and electrical conductivity, after going through a phase transition from 600 to 800 K. We calculate the transport coefficients as well as the defect concentrations ab initio, showing excellent agreement with experiment, and elucidating the origin of the double phase transition as well as the new charge carriers. AH Romero, EKU Gross, MJ Verstraete, and O Hellman PRB 91, 214310 (2015) O. Hellman, IA Abrikosov, and SI Simak, PRB 84 180301 (2011)

  11. The B AB AR detector

    NASA Astrophysics Data System (ADS)

    Aubert, B.; Bazan, A.; Boucham, A.; Boutigny, D.; De Bonis, I.; Favier, J.; Gaillard, J.-M.; Jeremie, A.; Karyotakis, Y.; Le Flour, T.; Lees, J. P.; Lieunard, S.; Petitpas, P.; Robbe, P.; Tisserand, V.; Zachariadou, K.; Palano, A.; Chen, G. P.; Chen, J. C.; Qi, N. D.; Rong, G.; Wang, P.; Zhu, Y. S.; Eigen, G.; Reinertsen, P. L.; Stugu, B.; Abbott, B.; Abrams, G. S.; Amerman, L.; Borgland, A. W.; Breon, A. B.; Brown, D. N.; Button-Shafer, J.; Clark, A. R.; Dardin, S.; Day, C.; Dow, S. F.; Fan, Q.; Gaponenko, I.; Gill, M. S.; Goozen, F. R.; Gowdy, S. J.; Gritsan, A.; Groysman, Y.; Hernikl, C.; Jacobsen, R. G.; Jared, R. C.; Kadel, R. W.; Kadyk, J.; Karcher, A.; Kerth, L. T.; Kipnis, I.; Kluth, S.; Kral, J. F.; Lafever, R.; LeClerc, C.; Levi, M. E.; Lewis, S. A.; Lionberger, C.; Liu, T.; Long, M.; Luo, L.; Lynch, G.; Luft, P.; Mandelli, E.; Marino, M.; Marks, K.; Matuk, C.; Meyer, A. B.; Minor, R.; Mokhtarani, A.; Momayezi, M.; Nyman, M.; Oddone, P. J.; Ohnemus, J.; Oshatz, D.; Patton, S.; Pedrali-Noy, M.; Perazzo, A.; Peters, C.; Pope, W.; Pripstein, M.; Quarrie, D. R.; Rasson, J. E.; Roe, N. A.; Romosan, A.; Ronan, M. T.; Shelkov, V. G.; Stone, R.; Strother, P. D.; Telnov, A. V.; von der Lippe, H.; Weber, T. F.; Wenzel, W. A.; Zizka, G.; Bright-Thomas, P. G.; Hawkes, C. M.; Kirk, A.; Knowles, D. J.; O'Neale, S. W.; Watson, A. T.; Watson, N. K.; Deppermann, T.; Koch, H.; Krug, J.; Kunze, M.; Lewandowski, B.; Peters, K.; Schmuecker, H.; Steinke, M.; Andress, J. C.; Barlow, N. R.; Bhimji, W.; Chevalier, N.; Clark, P. J.; Cottingham, W. N.; De Groot, N.; Dyce, N.; Foster, B.; Mass, A.; McFall, J. D.; Wallom, D.; Wilson, F. F.; Abe, K.; Hearty, C.; McKenna, J. A.; Thiessen, D.; Camanzi, B.; Harrison, T. J.; McKemey, A. K.; Tinslay, J.; Antohin, E. I.; Blinov, V. E.; Bukin, A. D.; Bukin, D. A.; Buzykaev, A. R.; Dubrovin, M. S.; Golubev, V. B.; Ivanchenko, V. N.; Kolachev, G. M.; Korol, A. A.; Kravchenko, E. A.; Mikhailov, S. F.; Onuchin, A. P.; Salnikov, A. A.; Serednyakov, S. I.; Skovpen, Yu. I.; Telnov, V. I.; Yushkov, A. N.; Booth, J.; Lankford, A. J.; Mandelkern, M.; Pier, S.; Stoker, D. P.; Zioulas, G.; Ahsan, A.; Arisaka, K.; Buchanan, C.; Chun, S.; Faccini, R.; MacFarlane, D. B.; Prell, S. A.; Rahatlou, Sh.; Raven, G.; Sharma, V.; Burke, S.; Callahan, D.; Campagnari, C.; Dahmes, B.; Hale, D.; Hart, P. A.; Kuznetsova, N.; Kyre, S.; Levy, S. L.; Long, O.; Lu, A.; May, J.; Richman, J. D.; Verkerke, W.; Witherell, M.; Yellin, S.; Beringer, J.; DeWitt, J.; Dorfan, D. E.; Eisner, A. M.; Frey, A.; Grillo, A. A.; Grothe, M.; Heusch, C. A.; Johnson, R. P.; Kroeger, W.; Lockman, W. S.; Pulliam, T.; Rowe, W.; Sadrozinski, H.; Schalk, T.; Schmitz, R. E.; Schumm, B. A.; Seiden, A.; Spencer, E. N.; Turri, M.; Walkowiak, W.; Wilder, M.; Williams, D. C.; Chen, E.; Dubois-Felsmann, G. P.; Dvoretskii, A.; Hanson, J. E.; Hitlin, D. G.; Kolomensky, Yu. G.; Metzler, S.; Oyang, J.; Porter, F. C.; Ryd, A.; Samuel, A.; Weaver, M.; Yang, S.; Zhu, R. Y.; Devmal, S.; Geld, T. L.; Jayatilleke, S.; Jayatilleke, S. M.; Mancinelli, G.; Meadows, B. T.; Sokoloff, M. D.; Bloom, P.; Broomer, B.; Erdos, E.; Fahey, S.; Ford, W. T.; Gaede, F.; van Hoek, W. C.; Johnson, D. R.; Michael, A. K.; Nauenberg, U.; Olivas, A.; Park, H.; Rankin, P.; Roy, J.; Sen, S.; Smith, J. G.; Wagner, D. L.; Blouw, J.; Harton, J. L.; Krishnamurthy, M.; Soffer, A.; Toki, W. H.; Warner, D. W.; Wilson, R. J.; Zhang, J.; Brandt, T.; Brose, J.; Dahlinger, G.; Dickopp, M.; Dubitzky, R. S.; Eckstein, P.; Futterschneider, H.; Kocian, M. L.; Krause, R.; Müller-Pfefferkorn, R.; Schubert, K. R.; Schwierz, R.; Spaan, B.; Wilden, L.; Behr, L.; Bernard, D.; Bonneaud, G. R.; Brochard, F.; Cohen-Tanugi, J.; Ferrag, S.; Fouque, G.; Gastaldi, F.; Matricon, P.; Mora de Freitas, P.; Renard, C.; Roussot, E.; T'Jampens, S.; Thiebaux, C.; Vasileiadis, G.; Verderi, M.; Anjomshoaa, A.; Bernet, R.; Di Lodovico, F.; Muheim, F.; Playfer, S.; Swain, J. E.; Falbo, M.; Bozzi, C.; Dittongo, S.; Folegani, M.; Piemontese, L.; Ramusino, A. C.; Treadwell, E.; Anulli, F.; Baldini-Ferroli, R.; Calcaterra, A.; de Sangro, R.; Falciai, D.; Finocchiaro, G.; Patteri, P.; Peruzzi, I. M.; Piccolo, M.; Xie, Y.; Zallo, A.; Bagnasco, S.; Buzzo, A.; Contri, R.; Crosetti, G.; Fabbricatore, P.; Farinon, S.; Lo Vetere, M.; Macri, M.; Minutoli, S.; Monge, M. R.; Musenich, R.; Pallavicini, M.; Parodi, R.; Passaggio, S.; Pastore, F. C.; Patrignani, C.; Pia, M. G.; Priano, C.; Robutti, E.; Santroni, A.; Bartoldus, R.; Dignan, T.; Hamilton, R.; Mallik, U.; Cochran, J.; Crawley, H. B.; Fischer, P. A.; Lamsa, J.; McKay, R.; Meyer, W. T.; Rosenberg, E. I.; Albert, J. N.; Beigbeder, C.; Benkebil, M.; Breton, D.; Cizeron, R.; Du, S.; Grosdidier, G.; Hast, C.; Höcker, A.; Lacker, H. M.; LePeltier, V.; Lutz, A. M.

    2002-02-01

    B AB AR, the detector for the SLAC PEP-II asymmetric e +e - B Factory operating at the ϒ(4 S) resonance, was designed to allow comprehensive studies of CP-violation in B-meson decays. Charged particle tracks are measured in a multi-layer silicon vertex tracker surrounded by a cylindrical wire drift chamber. Electromagnetic showers from electrons and photons are detected in an array of CsI crystals located just inside the solenoidal coil of a superconducting magnet. Muons and neutral hadrons are identified by arrays of resistive plate chambers inserted into gaps in the steel flux return of the magnet. Charged hadrons are identified by d E/d x measurements in the tracking detectors and by a ring-imaging Cherenkov detector surrounding the drift chamber. The trigger, data acquisition and data-monitoring systems, VME- and network-based, are controlled by custom-designed online software. Details of the layout and performance of the detector components and their associated electronics and software are presented.

  12. Recent progress in ab initio density matrix renormalization group methodology

    NASA Astrophysics Data System (ADS)

    Hachmann, Johannes; Dorando, Jonathan J.; Kin-Lic Chan, Garnet

    2008-03-01

    We present some recent developments in the ab initio density matrix renormalization group (DMRG) method for quantum chemical problems, in particular our local, quadratic scaling algorithm [1] for low dimensional systems. This method is particularly suited for the description of strong nondynamic correlation, and allows us to compute numerically exact (FCI) correlated energies for large active spaces, up to one order of magnitude larger then can be done by conventional CASCI techniques. Other features of this method are its inherent multireference nature, compactness, variational results, size-consistency and size-extensivity. In addition we will review the problems (predominantly organic electronic materials) on which we applied the ab initio DMRG: 1) metal-insulator transition in hydrogen chains [1] 2) all-trans polyacetylene [1] 3) acenes [2] 4) polydiacetylenes [3]. References [1] Hachmann, Cardoen, Chan, JCP 125 (2006), 144101. [2] Hachmann, Dorando, Avil'es, Chan, JCP 127 (2007), 134309. [3] unpublished.

  13. Ab Initio Studies of Calcium Carbonate Hydration.

    PubMed

    Lopez-Berganza, Josue A; Diao, Yijue; Pamidighantam, Sudhakar; Espinosa-Marzal, Rosa M

    2015-11-25

    Ab initio simulations of large hydrated calcium carbonate clusters are challenging due to the existence of multiple local energy minima. Extensive conformational searches around hydrated calcium carbonate clusters (CaCO3·nH2O for n = 1-18) were performed to find low-energy hydration structures using an efficient combination of Monte Carlo searches, density-functional tight binding (DFTB+) method, and density-functional theory (DFT) at the B3LYP level, or Møller-Plesset perturbation theory at the MP2 level. This multilevel optimization yields several low-energy structures for hydrated calcium carbonate. Structural and energetics analysis of the hydration of these clusters revealed a first hydration shell composed of 12 water molecules. Bond-length and charge densities were also determined for different cluster sizes. The solvation of calcium carbonate in bulk water was investigated by placing the explicitly solvated CaCO3·nH2O clusters in a polarizable continuum model (PCM). The findings of this study provide new insights into the energetics and structure of hydrated calcium carbonate and contribute to the understanding of mechanisms where calcium carbonate formation or dissolution is of relevance. PMID:26505205

  14. Characterization of Cry34Ab1 and Cry35Ab1 insecticidal crystal proteins expressed in transgenic corn plants and Pseudomonas fluorescens.

    PubMed

    Gao, Yong; Schafer, Barry W; Collins, Randy A; Herman, Rod A; Xu, Xiaoping; Gilbert, Jeffrey R; Ni, Weiting; Langer, Vickie L; Tagliani, Laura A

    2004-12-29

    Cry34Ab1 and Cry35Ab1 proteins, identified from Bacillus thuringiensis strain PS149B1, act together to control corn rootworms. Transgenic corn lines coexpressing the two proteins were developed to protect corn against rootworm damage. Large quantities of the two proteins were needed to conduct studies required for assessing the safety of this transgenic corn crop. Because it was technically infeasible to obtain sufficient quantities of high purity Cry34Ab1 and Cry35Ab1 proteins from the transgenic corn plants, the proteins were produced using a recombinant Pseudomonas fluorescens (Pf) production system. The two proteins from both the transgenic corn and the Pf were purified and characterized. The proteins from each host had the expected molecular mass and were immunoreactive to specific antibodies in enzyme-linked immunosorbent assay and Western blot analysis. Data from N-terminal sequencing, tryptic peptide mass fingerprinting, internal peptide sequencing, and biological activity provided direct evidence that the Cry34Ab1 and Cry35Ab1 proteins produced in Pf and transgenic corn were, respectively, comparable or equivalent molecules. In addition, neither protein had detectable glycosylation regardless of the host. PMID:15612796

  15. Microwave and ab initio studies of rare gas-methane van der Waals complexes

    NASA Astrophysics Data System (ADS)

    Liu, Yaqian; Jäger, Wolfgang

    2004-05-01

    Rotational spectra of the weakly bound Kr-methane van der Waals complex were recorded using a pulsed molecular beam Fourier transform microwave spectrometer in the range from 3.5 to 18 GHz. Spectra of 25 isotopomers of Kr-methane were assigned and analyzed. For isotopomers containing CH4, 13CH4, and CD4, two sets of transitions with K=0 and one with K=1 were recorded, correlating to the j=0, 1, and 2 rotational levels of free methane, respectively (j is the rotational angular momentum quantum number of the methane monomer). For isotopomers containing CH3D and CHD3, two K=0 components were recorded, correlating to the jk=00 and 11 rotational levels of free methane (k corresponds to the projection of j onto the C3 axis of CH3D and CHD3). The obtained spectroscopic results were used to derive van der Waals bond distance R, van der Waals stretching frequency νs, and the corresponding stretching force constant ks. Nuclear spin statistical weights of individual states were obtained from molecular symmetry group analyses and were compared with the observed relative transition intensities. The tentatively assigned j=2 transitions were more intense than predicted from symmetry considerations. This is attributed to a relatively large effective dipole moment of this state, supported by ab initio dipole moment calculations. Ab initio potential energy calculations of Kr-CH4 and Ar-CH4 were done at the coupled cluster level of theory, with single and double excitations and perturbative inclusion of triple excitations, using the aug-cc-pVTZ basis set supplemented with bond functions. The theoretical results show that the angular dynamics of the dimer does not change significantly when the binding partner of methane changes from Ar to Kr. The dipole moment of Ar-CH4 was calculated at various configurations, providing a qualitative explanation for the unsuccessful spectral searches for rotational transitions of Ar-CH4.

  16. Collective rotation from ab initio theory

    NASA Astrophysics Data System (ADS)

    Caprio, M. A.; Maris, P.; Vary, J. P.; Smith, R.

    2015-08-01

    Through ab initio approaches in nuclear theory, we may now seek to quantitatively understand the wealth of nuclear collective phenomena starting from the underlying internucleon interactions. No-core configuration interaction (NCCI) calculations for p-shell nuclei give rise to rotational bands, as evidenced by rotational patterns for excitation energies, electromagnetic moments and electromagnetic transitions. In this review, NCCI calculations of 7-9Be are used to illustrate and explore ab initio rotational structure, and the resulting predictions for rotational band properties are compared with experiment. We highlight the robustness of ab initio rotational predictions across different choices for the internucleon interaction.

  17. AB209630, a long non-coding RNA decreased expression in hypopharyngeal squamous cell carcinoma, influences proliferation, invasion, metastasis, and survival

    PubMed Central

    Li, Wenming; Wang, Juan; Xiang, Xuan; Li, Guojun; Pan, Xinliang; Lei, Dapeng

    2016-01-01

    Long noncoding RNAs (lncRNAs) are associated with the development, progression, and prognosis of human cancers. However, the clinical significance and biological function of lncRNAs in hypopharyngeal squamous cell carcinoma (HSCC) remain largely unknown. We characterized the novel lncRNA AB209630 in vivo and in vitro. First, using qRT-PCR, we evaluated whether AB209630 levels differ between HSCC tissues/cell lines and adjacent normal tissues/cell lines. We then assessed whether AB209630 expression levels stimulate or inhibit proliferation, invasion, apoptosis, and metastasis in vitro. Finally, we investigated whether AB209630 levels in tumor tissues were associated with survival outcomes. Our results demonstrated that AB209630 levels were markedly lower in HSCC tissues and cells than in normal tissues and cells, and increased expression of AB209630 level significantly inhibited growth, metastasis, and invasion and stimulated apoptosis in vitro. In addition, patients with decreased expression of AB209630 had a significantly poorer prognosis than those with high AB209630 expression. These data suggest that increased expression of AB209630 might either stimulate or inhibit biological activities involved in HSCC development, indicating a potential application of AB209630 in future treatment for this disease. This study suggest that AB209630 functions as a tumor suppressor in HSCC, and its decreased expression may help predict a poor prognostic outcome of HSCC. Our future work will focus on the mechanisms of whether and how AB209630 as a tumor suppressor gene is involved in HSCC development. PMID:26895099

  18. Ab Initio Study of Polonium

    SciTech Connect

    Zabidi, Noriza Ahmad; Kassim, Hasan Abu; Shrivastava, Keshav N.

    2008-05-20

    Polonium is the only element with a simple cubic (sc) crystal structure. Atoms in solid polonium sit at the corners of a simple cubic unit cell and no where else. Polonium has a valence electron configuration 6s{sup 2}6p{sup 4} (Z = 84). The low temperature {alpha}-phase transforms into the rhombohedral (trigonal) {beta} structure at {approx}348 K. The sc {alpha}-Po unit cell constant is a = 3.345 A. The beta form of polonium ({beta}-Po) has the lattice parameters, a{sub R} = 3.359 A and a rhombohedral angle 98 deg. 13'. We have performed an ab initio electronic structure calculation by using the density functional theory. We have performed the calculation with and without spin-orbit (SO) coupling by using both the LDA and the GGA for the exchange-correlations. The k-points in a simple cubic BZ are determined by R (0.5, 0.5, 0.5), {gamma} (0, 0, 0), X (0.5, 0, 0), M (0.5, 0.5, 0) and {gamma} (0, 0, 0). Other directions of k-points are {gamma} (0, 0, 0), X (0.5, 0, 0), R (0.5, 0.5, 0.5) and {gamma} (0, 0, 0). The SO splittings of p states at the {gamma} point in the GGA+SO scheme for {alpha}-Po are 0.04 eV and 0.02 eV while for the {beta}-Po these are 0.03 eV and 0.97 eV. We have also calculated the vibrational spectra for the unit cells in both the structures. We find that exchanging of a Po atom by Pb atom produces several more bands and destabilizes the {beta} phase.

  19. Human malignant mesothelioma is recapitulated in immunocompetent BALB/c mice injected with murine AB cells.

    PubMed

    Mezzapelle, Rosanna; Rrapaj, Eltjona; Gatti, Elena; Ceriotti, Chiara; Marchis, Francesco De; Preti, Alessandro; Spinelli, Antonello E; Perani, Laura; Venturini, Massimo; Valtorta, Silvia; Moresco, Rosa Maria; Pecciarini, Lorenza; Doglioni, Claudio; Frenquelli, Michela; Crippa, Luca; Recordati, Camilla; Scanziani, Eugenio; de Vries, Hilda; Berns, Anton; Frapolli, Roberta; Boldorini, Renzo; D'Incalci, Maurizio; Bianchi, Marco E; Crippa, Massimo P

    2016-01-01

    Malignant Mesothelioma is a highly aggressive cancer, which is difficult to diagnose and treat. Here we describe the molecular, cellular and morphological characterization of a syngeneic system consisting of murine AB1, AB12 and AB22 mesothelioma cells injected in immunocompetent BALB/c mice, which allows the study of the interplay of tumor cells with the immune system. Murine mesothelioma cells, like human ones, respond to exogenous High Mobility Group Box 1 protein, a Damage-Associated Molecular Pattern that acts as a chemoattractant for leukocytes and as a proinflammatory mediator. The tumors derived from AB cells are morphologically and histologically similar to human MM tumors, and respond to treatments used for MM patients. Our system largely recapitulates human mesothelioma, and we advocate its use for the study of MM development and treatment. PMID:26961782

  20. Human malignant mesothelioma is recapitulated in immunocompetent BALB/c mice injected with murine AB cells

    PubMed Central

    Mezzapelle, Rosanna; Rrapaj, Eltjona; Gatti, Elena; Ceriotti, Chiara; Marchis, Francesco De; Preti, Alessandro; Spinelli, Antonello E.; Perani, Laura; Venturini, Massimo; Valtorta, Silvia; Moresco, Rosa Maria; Pecciarini, Lorenza; Doglioni, Claudio; Frenquelli, Michela; Crippa, Luca; Recordati, Camilla; Scanziani, Eugenio; de Vries, Hilda; Berns, Anton; Frapolli, Roberta; Boldorini, Renzo; D’Incalci, Maurizio; Bianchi, Marco E.; Crippa, Massimo P.

    2016-01-01

    Malignant Mesothelioma is a highly aggressive cancer, which is difficult to diagnose and treat. Here we describe the molecular, cellular and morphological characterization of a syngeneic system consisting of murine AB1, AB12 and AB22 mesothelioma cells injected in immunocompetent BALB/c mice, which allows the study of the interplay of tumor cells with the immune system. Murine mesothelioma cells, like human ones, respond to exogenous High Mobility Group Box 1 protein, a Damage-Associated Molecular Pattern that acts as a chemoattractant for leukocytes and as a proinflammatory mediator. The tumors derived from AB cells are morphologically and histologically similar to human MM tumors, and respond to treatments used for MM patients. Our system largely recapitulates human mesothelioma, and we advocate its use for the study of MM development and treatment. PMID:26961782

  1. Ab initio SCF calculations on hydrogen bonded cresol isomers

    NASA Astrophysics Data System (ADS)

    Pohl, M.; Kleinermanns, K.

    1988-12-01

    Ab initio GAUSSIAN 80 calculations with two different basis sets (STO-3G and 4 31 G*) were performed on hydrogen bonded cresol isomers for comparison with experimental data from free jet fluorescence excitation spectroscopy. For m-cresol, the calculated barriers for hindered internal rotation of the OH-group and the CH3-group are in good agreement with experiment. The calculations show the trans-linear configuration of p-cresol· B-clusters ( B = H2O, CH3OH) to be more stable than the all-planar configuration. This agrees with CI calculations and microwave spectroscopic investigations of the water dimer. Calculations of both the intermolecular stretch and bend frequencies of p-cresol· B-clusters show little dependence on the all-planar or trans-linear configuration but a strong dependence on the choice of the basis set. With the minimal basis set STO-3G, the vibrational energies are generally too high. The agreement between the calculated vibrational frequencies from the 4 31 G* basis set and the experimental values is fair.

  2. Basis of Articulation.

    ERIC Educational Resources Information Center

    Kelz, Heinrich

    This article intends to shed light on the somewhat nebulous term "basis of articulation," which is used frequently in Eastern European phonetic and linguistic literature but highly neglected in contemporary American literature. In a historical approach, it is shown how the term originated and developed, how it is defined by various authors, and…

  3. Ab initio molecular dynamics of liquid water using embedded-fragment second-order many-body perturbation theory towards its accurate property prediction

    PubMed Central

    Willow, Soohaeng Yoo; Salim, Michael A.; Kim, Kwang S.; Hirata, So

    2015-01-01

    A direct, simultaneous calculation of properties of a liquid using an ab initio electron-correlated theory has long been unthinkable. Here we present structural, dynamical, and response properties of liquid water calculated by ab initio molecular dynamics using the embedded-fragment spin-component-scaled second-order many-body perturbation method with the aug-cc-pVDZ basis set. This level of theory is chosen as it accurately and inexpensively reproduces the water dimer potential energy surface from the coupled-cluster singles, doubles, and noniterative triples with the aug-cc-pVQZ basis set, which is nearly exact. The calculated radial distribution function, self-diffusion coefficient, coordinate number, and dipole moment, as well as the infrared and Raman spectra are in excellent agreement with experimental results. The shapes and widths of the OH stretching bands in the infrared and Raman spectra and their isotropic-anisotropic Raman noncoincidence, which reflect the diverse local hydrogen-bond environment, are also reproduced computationally. The simulation also reveals intriguing dynamic features of the environment, which are difficult to probe experimentally, such as a surprisingly large fluctuation in the coordination number and the detailed mechanism by which the hydrogen donating water molecules move across the first and second shells, thereby causing this fluctuation. PMID:26400690

  4. Ab initio molecular dynamics of liquid water using embedded-fragment second-order many-body perturbation theory towards its accurate property prediction.

    PubMed

    Willow, Soohaeng Yoo; Salim, Michael A; Kim, Kwang S; Hirata, So

    2015-01-01

    A direct, simultaneous calculation of properties of a liquid using an ab initio electron-correlated theory has long been unthinkable. Here we present structural, dynamical, and response properties of liquid water calculated by ab initio molecular dynamics using the embedded-fragment spin-component-scaled second-order many-body perturbation method with the aug-cc-pVDZ basis set. This level of theory is chosen as it accurately and inexpensively reproduces the water dimer potential energy surface from the coupled-cluster singles, doubles, and noniterative triples with the aug-cc-pVQZ basis set, which is nearly exact. The calculated radial distribution function, self-diffusion coefficient, coordinate number, and dipole moment, as well as the infrared and Raman spectra are in excellent agreement with experimental results. The shapes and widths of the OH stretching bands in the infrared and Raman spectra and their isotropic-anisotropic Raman noncoincidence, which reflect the diverse local hydrogen-bond environment, are also reproduced computationally. The simulation also reveals intriguing dynamic features of the environment, which are difficult to probe experimentally, such as a surprisingly large fluctuation in the coordination number and the detailed mechanism by which the hydrogen donating water molecules move across the first and second shells, thereby causing this fluctuation. PMID:26400690

  5. A Complete and Accurate Ab Initio Repeat Finding Algorithm.

    PubMed

    Lian, Shuaibin; Chen, Xinwu; Wang, Peng; Zhang, Xiaoli; Dai, Xianhua

    2016-03-01

    It has become clear that repetitive sequences have played multiple roles in eukaryotic genome evolution including increasing genetic diversity through mutation, changes in gene expression and facilitating generation of novel genes. However, identification of repetitive elements can be difficult in the ab initio manner. Currently, some classical ab initio tools of finding repeats have already presented and compared. The completeness and accuracy of detecting repeats of them are little pool. To this end, we proposed a new ab initio repeat finding tool, named HashRepeatFinder, which is based on hash index and word counting. Furthermore, we assessed the performances of HashRepeatFinder with other two famous tools, such as RepeatScout and Repeatfinder, in human genome data hg19. The results indicated the following three conclusions: (1) The completeness of HashRepeatFinder is the best one among these three compared tools in almost all chromosomes, especially in chr9 (8 times of RepeatScout, 10 times of Repeatfinder); (2) in terms of detecting large repeats, HashRepeatFinder also performed best in all chromosomes, especially in chr3 (24 times of RepeatScout and 250 times of Repeatfinder) and chr19 (12 times of RepeatScout and 60 times of Repeatfinder); (3) in terms of accuracy, HashRepeatFinder can merge the abundant repeats with high accuracy. PMID:26272474

  6. Ab initio and DFT study of the geometric structures and static dipole (hyper)polarizabilities of aromatic anions.

    PubMed

    Castellano, O; Bermúdez, Y; Giffard, M; Mabon, G; Cubillan, N; Sylla, M; Nguyen-Phu, X; Hinchliffe, A; Soscún, H

    2005-11-17

    The geometries and the static dipole (hyper)polarizabilities (alpha, beta, gamma) of a series of aromatic anions were investigated at the ab initio (HF, MP2, and MP4) and density functional theory DFT (B3LYP) levels of theory. The anions chosen for the present study are the benzenethiolate (Ph-S-), benzenecarboxylate (Ph-CO2-), benzenesulfinate (Ph-SO2-), benzenesulfonate (Ph-SO3-), and 1,3-benzenedicarboxylate (1,3-Ph-(CO2)2(2-)). For benzenethiolate anion, additional alpha, beta, and gamma calculations were performed at the coupled cluster CCSD level with MP2 optimized geometries. The standard diffuse and polarized 6-31+G(d,p) basis set was employed in conjunction to the ab initio and DFT methods. Additional HF calculations were performed with the 6-311++G(3d,3p) basis set for all the anions. The correlated electric properties were evaluated numerically within the formalism of finite field. The optimized geometries were analyzed in terms of the few reports about the phenolate and sulfonate ions. The results show that electron correlation effects on the polarizabilities are very important in all the anion series. Was found that Ph-SO2- is highly polarizable in terms of alpha and beta, and the Ph-S- is the highest second hyperpolarizable in the series. The results of alpha were rationalized in terms of the analysis of the polarization of charge based in Mulliken atomic population and the structural features of the optimized geometries of anions, whereas the large differences in the beta and gamma values in the series were respectively interpreted in terms of the bond length alternation BLA and the separation of charge in the aromatic ring by effects of the substitution. These results allowed us to suggest the benzenesulfinate and benzenethiolate anions as promising candidates that should be incorporated in ionic materials for second and third-order nonlinear optical devices. PMID:16833334

  7. Argon Interaction with Gold Surfaces: Ab Initio-Assisted Determination of Pair Ar-Au Potentials for Molecular Dynamics Simulations.

    PubMed

    Grenier, Romain; To, Quy-Dong; de Lara-Castells, María Pilar; Léonard, Céline

    2015-07-01

    Global potentials for the interaction between the Ar atom and gold surfaces are investigated and Ar-Au pair potentials suitable for molecular dynamics simulations are derived. Using a periodic plane-wave representation of the electronic wave function, the nonlocal van-der-Waals vdW-DF2 and vdW-OptB86 approaches have been proved to describe better the interaction. These global interaction potentials have been decomposed to produce pair potentials. Then, the pair potentials have been compared with those derived by combining the dispersionless density functional dlDF for the repulsive part with an effective pairwise dispersion interaction. These repulsive potentials have been obtained from the decomposition of the repulsive interaction between the Ar atom and the Au2 and Au4 clusters and the dispersion coefficients have been evaluated by means of ab initio calculations on the Ar+Au2 complex using symmetry adapted perturbation theory. The pair potentials agree very well with those evaluated through periodic vdW-DF2 calculations. For benchmarking purposes, CCSD(T) calculations have also been performed for the ArAu and Ar+Au2 systems using large basis sets and extrapolations to the complete basis set limit. This work highlights that ab initio calculations using very small surface clusters can be used either as an independent cross-check to compare the performance of state-of-the-art vdW-corrected periodic DFT approaches or, directly, to calculate the pair potentials necessary in further molecular dynamics calculations. PMID:26046588

  8. Anatomical repair of large incisional hernias.

    PubMed Central

    Loh, A.; Rajkumar, J. S.; South, L. M.

    1992-01-01

    We present a method of repair for large incisional hernias using lateral relieving incisions of the anterior rectus sheath. This is a modification of the methods previously described by Young (1), Hunter (2) and Maguire and Young (3). There were no recurrences in the 13 patients reviewed. Other methods of repair for large incisional hernias are discussed. Images Figure 2a,b Figure 3a,b Figure 4 Figure 5 PMID:1567126

  9. Formation of large voids in the amorphous phase-change memory Ge2Sb2Te5 alloy.

    PubMed

    Sun, Zhimei; Zhou, Jian; Blomqvist, Andreas; Johansson, Börje; Ahuja, Rajeev

    2009-02-20

    On the basis of ab initio molecular dynamics simulations, large voids mainly surrounded by Te atoms are observed in molten and amorphous Ge2Sb2Te5, which is due to the clustering of two- and threefold coordinated Te atoms. Furthermore, pressure shows a significant effect on the clustering of the under coordinated Te atoms and hence the formation of large voids. The present results demonstrate that both vacancies and Te play an important role in the fast reversible phase transition process. PMID:19257687

  10. Carbonate-Bridged Lanthanoid Triangles: Single-Molecule Magnet Behavior, Inelastic Neutron Scattering, and Ab Initio Studies.

    PubMed

    Giansiracusa, Marcus J; Vonci, Michele; Van den Heuvel, Willem; Gable, Robert W; Moubaraki, Boujemaa; Murray, Keith S; Yu, Dehong; Mole, Richard A; Soncini, Alessandro; Boskovic, Colette

    2016-06-01

    comprised of contributions from large MJ, giving rise to slow magnetic relaxation. Although no direct evidence for intramolecular RE···RE magnetic coupling is observed in either magnetic or INS studies, on the basis of the ab initio calculations, we find noncollinear magnetic axes in 1-Er that are coplanar with the erbium triangle and radially arranged with respect to the triangle's centroid; thus, we argue that the absence of magnetic coupling in this system arises from dipolar and antiferromagnetic superexchange interactions that cancel each other out. PMID:27203849

  11. Determination of organic priority pollutants and emerging compounds in wastewater and snow samples using multiresidue protocols on the basis of microextraction by packed sorbents coupled to large volume injection gas chromatography-mass spectrometry analysis.

    PubMed

    Prieto, A; Schrader, S; Moeder, M

    2010-09-17

    This paper describes the development and validation of a new procedure for the simultaneous determination of 41 multi-class priority and emerging organic pollutants in water samples using microextraction by packed sorbent (MEPS) followed by large volume injection-gas chromatography-mass spectrometry (LVI-GC-MS). Apart from method parameter optimization the influence of humic acids as matrix components on the extraction efficiency of MEPS procedure was also evaluated. The list of target compounds includes polycyclic aromatic hydrocarbons (PAHs), polychlorinated biphenyls (PCBs), phthalate esters (PEs), nonylphenols (NPs), bisphenol A (BPA) and selected steroid hormones. The performance of the new at-line microextraction-LVI-GC-MS protocol was compared to standard solid-phase extraction (SPE) and LVI-GC-MS analysis. LODs for 100 mL samples (SPE) ranged from 0.2 to 736 ng L(-1) were obtained. LODs for 800 microL of sample (MEPS) were between 0.2 and 266 ng L(-1). In the case of MEPS methodology even a sample volume of only 800 microL allowed to detect the target compounds. These results demonstrate the high sensitivity of both procedures which permitted to obtain good recoveries (>75%) for all cases. The precision of the methods, calculated as relative standard deviation (RSD) was below 21% for all compounds and both methodologies. Finally, the developed methods were applied to the determination of target analytes in various samples, including snow and wastewater. PMID:20719318

  12. Symmetry-guided large-scale shell-model theory

    NASA Astrophysics Data System (ADS)

    Launey, Kristina D.; Dytrych, Tomas; Draayer, Jerry P.

    2016-07-01

    In this review, we present a symmetry-guided strategy that utilizes exact as well as partial symmetries for enabling a deeper understanding of and advancing ab initio studies for determining the microscopic structure of atomic nuclei. These symmetries expose physically relevant degrees of freedom that, for large-scale calculations with QCD-inspired interactions, allow the model space size to be reduced through a very structured selection of the basis states to physically relevant subspaces. This can guide explorations of simple patterns in nuclei and how they emerge from first principles, as well as extensions of the theory beyond current limitations toward heavier nuclei and larger model spaces. This is illustrated for the ab initio symmetry-adapted no-core shell model (SA-NCSM) and two significant underlying symmetries, the symplectic Sp(3 , R) group and its deformation-related SU(3) subgroup. We review the broad scope of nuclei, where these symmetries have been found to play a key role-from the light p-shell systems, such as 6Li, 8B, 8Be, 12C, and 16O, and sd-shell nuclei exemplified by 20Ne, based on first-principle explorations; through the Hoyle state in 12C and enhanced collectivity in intermediate-mass nuclei, within a no-core shell-model perspective; up to strongly deformed species of the rare-earth and actinide regions, as investigated in earlier studies. A complementary picture, driven by symmetries dual to Sp(3 , R) , is also discussed. We briefly review symmetry-guided techniques that prove useful in various nuclear-theory models, such as Elliott model, ab initio SA-NCSM, symplectic model, pseudo- SU(3) and pseudo-symplectic models, ab initio hyperspherical harmonics method, ab initio lattice effective field theory, exact pairing-plus-shell model approaches, and cluster models, including the resonating-group method. Important implications of these approaches that have deepened our understanding of emergent phenomena in nuclei, such as enhanced

  13. Water body mapping method with HJ-1A/B satellite imagery

    NASA Astrophysics Data System (ADS)

    Lu, Shanlong; Wu, Bingfang; Yan, Nana; Wang, Hao

    2011-06-01

    This paper proposes an integrated water body mapping method with HJ-1A/B satellite imagery, the CCD (charge coupled device) data of the Chinese environmental satellites that were launched on September 6th, 2008. It combines the difference between NDVI and NDWI (NDVI-NDWI) with SLOPE and near-infrared (NIR) band. The NDVI-NDWI index is used to enhance the contrast between water bodies and the surrounding surface features; the topographic SLOPE is used to eliminate the mountain shadow; and the NIR band is used to reduce the effects of artificial construction land. The objectives are evaluating the potential of the HJ-1A/B imagery on water body monitoring, and proposing ideally mapping method. The test study results indicated that the NDVI-NDWI index is superior to the single index of NDVI and NDWI to enhance the contrast between water bodies and the rest of the features. On the basis of the accurately mapped water bodies in the HJ-1A/B CCD images of the study area, we conclude that the HJ-1A/B multi-spectral satellite images is an ideal data source for high spatial and temporal resolution water bodies monitoring. And the integrated water body mapping method is suitable for the applications of HJ-1A/B multi-spectral satellite images in this field.

  14. Optimized energy landscape exploration using the ab initio based activation-relaxation technique.

    PubMed

    Machado-Charry, Eduardo; Béland, Laurent Karim; Caliste, Damien; Genovese, Luigi; Deutsch, Thierry; Mousseau, Normand; Pochet, Pascal

    2011-07-21

    Unbiased open-ended methods for finding transition states are powerful tools to understand diffusion and relaxation mechanisms associated with defect diffusion, growth processes, and catalysis. They have been little used, however, in conjunction with ab initio packages as these algorithms demanded large computational effort to generate even a single event. Here, we revisit the activation-relaxation technique (ART nouveau) and introduce a two-step convergence to the saddle point, combining the previously used Lanczós algorithm with the direct inversion in interactive subspace scheme. This combination makes it possible to generate events (from an initial minimum through a saddle point up to a final minimum) in a systematic fashion with a net 300-700 force evaluations per successful event. ART nouveau is coupled with BigDFT, a Kohn-Sham density functional theory (DFT) electronic structure code using a wavelet basis set with excellent efficiency on parallel computation, and applied to study the potential energy surface of C(20) clusters, vacancy diffusion in bulk silicon, and reconstruction of the 4H-SiC surface. PMID:21786982

  15. ``Phantom'' Modes in Ab Initio Tunneling Calculations: Implications for Theoretical Materials Optimization, Tunneling, and Transport

    NASA Astrophysics Data System (ADS)

    Barabash, Sergey V.; Pramanik, Dipankar

    2015-03-01

    Development of low-leakage dielectrics for semiconductor industry, together with many other areas of academic and industrial research, increasingly rely upon ab initio tunneling and transport calculations. Complex band structure (CBS) is a powerful formalism to establish the nature of tunneling modes, providing both a deeper understanding and a guided optimization of materials, with practical applications ranging from screening candidate dielectrics for lowest ``ultimate leakage'' to identifying charge-neutrality levels and Fermi level pinning. We demonstrate that CBS is prone to a particular type of spurious ``phantom'' solution, previously deemed true but irrelevant because of a very fast decay. We demonstrate that (i) in complex materials, phantom modes may exhibit very slow decay (appearing as leading tunneling terms implying qualitative and huge quantitative errors), (ii) the phantom modes are spurious, (iii) unlike the pseudopotential ``ghost'' states, phantoms are an apparently unavoidable artifact of large numerical basis sets, (iv) a presumed increase in computational accuracy increases the number of phantoms, effectively corrupting the CBS results despite the higher accuracy achieved in resolving the true CBS modes and the real band structure, and (v) the phantom modes cannot be easily separated from the true CBS modes. We discuss implications for direct transport calculations. The strategy for dealing with the phantom states is discussed in the context of optimizing high-quality high- κ dielectric materials for decreased tunneling leakage.

  16. Low-energy states of manganese-oxo corrole and corrolazine: multiconfiguration reference ab initio calculations.

    PubMed

    Zhao, Hailiang; Pierloot, Kristine; Langner, Ernie H G; Swarts, Jannie C; Conradie, Jeanet; Ghosh, Abhik

    2012-04-01

    Manganese(V)-oxo corrole and corrolazine have been studied with ab initio multiconfiguration reference methods (CASPT2 and RASPT2) and large atomic natural orbital (ANO) basis sets. The calculations confirm the expected singlet d(δ)(2) ground states for both complexes and rule out excited states within 0.5 eV of the ground states. The lowest excited states are a pair of Mn(V) triplet states with d(δ)(1)d(π)(1) configurations 0.5-0.75 eV above the ground state. Manganese(IV)-oxo macrocycle radical states are much higher in energy, ≥1.0 eV relative to the ground state. The macrocyclic ligands in the ground states of the complexes are thus unambiguously 'innocent'. The approximate similarity of the spin state energetics of the corrole and corrolazine complexes suggests that the latter macrocycle on its own does not afford any special stabilization for the Mn(V)O center. The remarkable stability of an Mn(V)O octaarylcorrolazine thus appears to be ascribable to the steric protection afforded by the β-aryl groups. PMID:22432719

  17. Efficient Ab initio Modeling of Random Multicomponent Alloys

    NASA Astrophysics Data System (ADS)

    Jiang, Chao; Uberuaga, Blas P.

    2016-03-01

    We present in this Letter a novel small set of ordered structures (SSOS) method that allows extremely efficient ab initio modeling of random multicomponent alloys. Using inverse II-III spinel oxides and equiatomic quinary bcc (so-called high entropy) alloys as examples, we demonstrate that a SSOS can achieve the same accuracy as a large supercell or a well-converged cluster expansion, but with significantly reduced computational cost. In particular, because of this efficiency, a large number of quinary alloy compositions can be quickly screened, leading to the identification of several new possible high-entropy alloy chemistries. The SSOS method developed here can be broadly useful for the rapid computational design of multicomponent materials, especially those with a large number of alloying elements, a challenging problem for other approaches.

  18. Efficient Ab initio Modeling of Random Multicomponent Alloys.

    PubMed

    Jiang, Chao; Uberuaga, Blas P

    2016-03-11

    We present in this Letter a novel small set of ordered structures (SSOS) method that allows extremely efficient ab initio modeling of random multicomponent alloys. Using inverse II-III spinel oxides and equiatomic quinary bcc (so-called high entropy) alloys as examples, we demonstrate that a SSOS can achieve the same accuracy as a large supercell or a well-converged cluster expansion, but with significantly reduced computational cost. In particular, because of this efficiency, a large number of quinary alloy compositions can be quickly screened, leading to the identification of several new possible high-entropy alloy chemistries. The SSOS method developed here can be broadly useful for the rapid computational design of multicomponent materials, especially those with a large number of alloying elements, a challenging problem for other approaches. PMID:27015491

  19. Radioactive Waste Management Basis

    SciTech Connect

    Perkins, B K

    2009-06-03

    The purpose of this Radioactive Waste Management Basis is to describe the systematic approach for planning, executing, and evaluating the management of radioactive waste at LLNL. The implementation of this document will ensure that waste management activities at LLNL are conducted in compliance with the requirements of DOE Order 435.1, Radioactive Waste Management, and the Implementation Guide for DOE Manual 435.1-1, Radioactive Waste Management Manual. Technical justification is provided where methods for meeting the requirements of DOE Order 435.1 deviate from the DOE Manual 435.1-1 and Implementation Guide.

  20. TOPICAL REVIEW: Ab initio symplectic no-core shell model

    NASA Astrophysics Data System (ADS)

    Dytrych, T.; Sviratcheva, K. D.; Draayer, J. P.; Bahri, C.; Vary, J. P.

    2008-12-01

    The no-core shell model (NCSM) is a prominent ab initio method that yields a good description of the low-lying states in few-nucleon systems as well as in more complex p-shell nuclei. Nevertheless, its applicability is limited by the rapid growth of the many-body basis with larger model spaces and increasing number of nucleons. The symplectic no-core shell model (Sp-NCSM) aspires to extend the scope of the NCSM beyond the p-shell region by augmenting the conventional spherical harmonic oscillator basis with the physically relevant symplectic \\SpR{3} symmetry-adapted configurations of the symplectic shell model that describe naturally the monopole-quadrupole vibrational and rotational modes, and also partially incorporate α-cluster correlations. In this review, the models underpinning the Sp-NCSM approach, namely, the NCSM, the Elliott SU(3) model and the symplectic shell model, are discussed. Following this, a prescription for constructing translationally invariant symplectic configurations in the spherical harmonic oscillator basis is given. This prescription is utilized to unveil the extent to which symplectic configurations enter into low-lying states in 12C and 16O nuclei calculated within the framework of the NCSM with the JISP16 realistic nucleon-nucleon interaction. The outcomes of this proof-of-principle study are presented in detail.

  1. Ab interno trabeculectomy: patient selection and perspectives.

    PubMed

    Vinod, Kateki; Gedde, Steven J

    2016-01-01

    Ab interno trabeculectomy is one among several recently introduced minimally invasive glaucoma surgeries that avoid a conjunctival incision and full-thickness sclerostomy involved in traditional glaucoma surgery. Ablation of the trabecular meshwork and inner wall of Schlemm's canal is performed in an arcuate fashion via a clear corneal incision, alone or in combination with phacoemulsification cataract surgery. Intraocular pressure reduction following ab interno trabeculectomy is limited by resistance in distal outflow pathways and generally stabilizes in the mid-to-high teens. Relief of medication burden has been demonstrated by some studies. A very low rate of complications, most commonly transient hyphema and intraocular pressure elevations in the immediate postoperative period, have been reported. However, available data are derived from small retrospective and prospective case series. Randomized, controlled trials are needed to better elucidate the potential merits of ab interno trabeculectomy in the combined setting versus phacoemulsification cataract surgery alone and to compare it with other minimally invasive glaucoma surgeries. PMID:27574396

  2. Gravitational Spacecraft Anomalies as well as the at Present Relatively large Uncertainty of Newton's Gravity Constant are Explained on the Basis of Force-Effects Due to a so-far Unknown Form of Space-like Matter

    NASA Astrophysics Data System (ADS)

    Volkamer, Klaus

    2009-09-01

    Recently, a so-far unkown form of invisible, space-like and field-like form of matter was detected with real, weighable mass content. The observed quanta of field-like matter exhibit in laboratory experiments Planck mass mp = ℏ.c/G = 21.77 μg. They show either a positive or negative sign, and can be understood as candidates for dark matter (field-like quanta with positive sign) and dark energy (field-like quanta with negative sign). Due to the observed gravitational as well as ``topological'' (i.e. form-specific at phase borders) interaction of space-like matter, celestial objects at various scales build up quantized fields of these forms of matter around their centers' of gravity of normal matter which reach far beyond the observable surface of such objects, respectively. From the description of he space-like matter fields of the Sun and the Earth a quantitative explanation of the reported NASA spacecraft anomalies of Pioneer 10 and 11 as well as of NEAR-Shoemaker at its fly-by maneuver on 01/23/1998 at Earth is given. The structure of such quantized subtle matter fields of stars allow in addition a description of the formation of ``normal'' and ``inverse'' planetary nebulae at the end of a star's life.-So-far unknown physical force-effects between subtle matter fields, bound, due to the topological interaction, at various metals were observed. The results allow an explanation of the at present relatively large uncertainty ΔG/G in the determination of Newton's constant of gravity. Devices for a more precise and reproducible determination of G (with an accuracy comparable to that of Planck's quantum of action, for example) should be made of beryllium or aluminium/berylUum alloy to eliminate subtle matter effects or such effects should be taken into account by a comprehensive theoretical modeling of their material and shape (and also to some extend time and place) depending force actions. Thus, both anomalies, the reported NASA spacecraft acceleration anomalies

  3. Ab initio computational study of beta-cellobiose conformers using B3LYP/6-311++G**.

    PubMed

    Strati, Gina L; Willett, Julious L; Momany, Frank A

    2002-11-01

    The molecular structure of 27 conformers of beta-cellobiose were studied in vacuo through gradient geometry optimization using B3LYP density functionals and the 6-311++G** basis set. The conformationally dependent geometry changes and energies were explored as well as the hydrogen-bonding network. The lowest electronic energy structures found were not those suggested from available crystallographic and NMR solution data, where the glycosidic dihedral angles fall in the region (phi, psi) approximately (40 degrees, -20 degrees ). Rather, 'flipped' conformations in which the dihedral angles are in the range (phi, psi) approximately (180 degrees, 0 degrees ) are energetically more stable by approximately 2.5 kcal/mol over the 'experimentally accepted' structure. Further, when the vibrational free energy, deltaG, obtained from the calculated frequencies, is compared throughout the series, structures with (phi, psi) in the experimentally observed range still have higher free energy ( approximately 2.0 kcal/mol) than 'flipped' forms. The range of bridging dihedral angles of the 'normal' conformers, resulting from the variance in the phi dihedral is larger than that found in the 'flipped' forms. Due to this large flat energy surface for the normal conformations, we surmise that the summation of populations of these conformations will favor the 'normal' conformations, although evidence suggests that polar solvent effects may play the dominant role in providing stability for the 'normal' forms. Even though some empirical studies previously found the 'flipped' conformations to be lowest in energy, these studies have been generally discredited because they were in disagreement with experimental results. Most of the DFT/ab initio conformations reported here have not been reported previously in the ab initio literature, in part because the use of less rigorous theoretical methods, i.e. smaller basis sets, have given results in general agreement with experimental data, that is

  4. Characterization of mAb dimers reveals predominant dimer forms common in therapeutic mAbs.

    PubMed

    Plath, Friederike; Ringler, Philippe; Graff-Meyer, Alexandra; Stahlberg, Henning; Lauer, Matthias E; Rufer, Arne C; Graewert, Melissa A; Svergun, Dmitri; Gellermann, Gerald; Finkler, Christof; Stracke, Jan O; Koulov, Atanas; Schnaible, Volker

    2016-07-01

    The formation of undesired high molecular weight species such as dimers is an important quality attribute for therapeutic monoclonal antibody formulations. Therefore, the thorough understanding of mAb dimerization and the detailed characterization mAb dimers is of great interest for future pharmaceutical development of therapeutic antibodies. In this work, we focused on the analyses of different mAb dimers regarding size, surface properties, chemical identity, overall structure and localization of possible dimerization sites. Dimer fractions of different mAbs were isolated to a satisfactory purity from bulk material and revealed 2 predominant overall structures, namely elongated and compact dimer forms. The elongated dimers displayed one dimerization site involving the tip of the Fab domain. Depending on the stress applied, these elongated dimers are connected either covalently or non-covalently. In contrast, the compact dimers exhibited non-covalent association. Several interaction points were detected for the compact dimers involving the hinge region or the base of the Fab domain. These results indicate that mAb dimer fractions are rather complex and may contain more than one kind of dimer. Nevertheless, the overall appearance of mAb dimers suggests the existence of 2 predominant dimeric structures, elongated and compact, which are commonly present in preparations of therapeutic mAbs. PMID:27031922

  5. Operator evolution for ab initio theory of light nuclei

    NASA Astrophysics Data System (ADS)

    Schuster, Micah; Quaglioni, Sofia; Johnson, Calvin; Jurgenson, Eric; Navrátil, Petr

    2014-09-01

    The past two decades have seen a revolution in ab initio calculations of nuclear properties. One key element has been the development of a rigorous effective interaction theory, applying unitary transformations to soften the nuclear Hamiltonian and hence accelerate the convergence as a function of the model space size. For consistency, however, one ought to apply the same transformation to other operators when calculating transitions and mean values from the eigenstates of the renormalized Hamiltonian. Working in a translationally invariant harmonic oscillator basis for the two- and three-nucleon systems, we evolve the Hamiltonian, square radius, and total dipole strength operators by the similarity renormalization group (SRG). The inclusion of up to three-body matrix elements in the 4He nucleus all but completely restores the invariance of the expectation values under the transformation. We also consider a Gaussian operator with adjustable range; short ranges have the largest absolute renormalization when including two- and three-body induced terms, while at long ranges the induced three-body contribution takes on increased relative importance. The past two decades have seen a revolution in ab initio calculations of nuclear properties. One key element has been the development of a rigorous effective interaction theory, applying unitary transformations to soften the nuclear Hamiltonian and hence accelerate the convergence as a function of the model space size. For consistency, however, one ought to apply the same transformation to other operators when calculating transitions and mean values from the eigenstates of the renormalized Hamiltonian. Working in a translationally invariant harmonic oscillator basis for the two- and three-nucleon systems, we evolve the Hamiltonian, square radius, and total dipole strength operators by the similarity renormalization group (SRG). The inclusion of up to three-body matrix elements in the 4He nucleus all but completely restores

  6. WEST Physics Basis

    NASA Astrophysics Data System (ADS)

    Bourdelle, C.; Artaud, J. F.; Basiuk, V.; Bécoulet, M.; Brémond, S.; Bucalossi, J.; Bufferand, H.; Ciraolo, G.; Colas, L.; Corre, Y.; Courtois, X.; Decker, J.; Delpech, L.; Devynck, P.; Dif-Pradalier, G.; Doerner, R. P.; Douai, D.; Dumont, R.; Ekedahl, A.; Fedorczak, N.; Fenzi, C.; Firdaouss, M.; Garcia, J.; Ghendrih, P.; Gil, C.; Giruzzi, G.; Goniche, M.; Grisolia, C.; Grosman, A.; Guilhem, D.; Guirlet, R.; Gunn, J.; Hennequin, P.; Hillairet, J.; Hoang, T.; Imbeaux, F.; Ivanova-Stanik, I.; Joffrin, E.; Kallenbach, A.; Linke, J.; Loarer, T.; Lotte, P.; Maget, P.; Marandet, Y.; Mayoral, M. L.; Meyer, O.; Missirlian, M.; Mollard, P.; Monier-Garbet, P.; Moreau, P.; Nardon, E.; Pégourié, B.; Peysson, Y.; Sabot, R.; Saint-Laurent, F.; Schneider, M.; Travère, J. M.; Tsitrone, E.; Vartanian, S.; Vermare, L.; Yoshida, M.; Zagorski, R.; Contributors, JET

    2015-06-01

    With WEST (Tungsten Environment in Steady State Tokamak) (Bucalossi et al 2014 Fusion Eng. Des. 89 907-12), the Tore Supra facility and team expertise (Dumont et al 2014 Plasma Phys. Control. Fusion 56 075020) is used to pave the way towards ITER divertor procurement and operation. It consists in implementing a divertor configuration and installing ITER-like actively cooled tungsten monoblocks in the Tore Supra tokamak, taking full benefit of its unique long-pulse capability. WEST is a user facility platform, open to all ITER partners. This paper describes the physics basis of WEST: the estimated heat flux on the divertor target, the planned heating schemes, the expected behaviour of the L-H threshold and of the pedestal and the potential W sources. A series of operating scenarios has been modelled, showing that ITER-relevant heat fluxes on the divertor can be achieved in WEST long pulse H-mode plasmas.

  7. Genetic basis of hyperlysinemia

    PubMed Central

    2013-01-01

    Background Hyperlysinemia is an autosomal recessive inborn error of L-lysine degradation. To date only one causal mutation in the AASS gene encoding α-aminoadipic semialdehyde synthase has been reported. We aimed to better define the genetic basis of hyperlysinemia. Methods We collected the clinical, biochemical and molecular data in a cohort of 8 hyperlysinemia patients with distinct neurological features. Results We found novel causal mutations in AASS in all affected individuals, including 4 missense mutations, 2 deletions and 1 duplication. In two patients originating from one family, the hyperlysinemia was caused by a contiguous gene deletion syndrome affecting AASS and PTPRZ1. Conclusions Hyperlysinemia is caused by mutations in AASS. As hyperlysinemia is generally considered a benign metabolic variant, the more severe neurological disease course in two patients with a contiguous deletion syndrome may be explained by the additional loss of PTPRZ1. Our findings illustrate the importance of detailed biochemical and genetic studies in any hyperlysinemia patient. PMID:23570448

  8. Genetic basis of autoimmunity

    PubMed Central

    Marson, Alexander; Housley, William J.; Hafler, David A.

    2015-01-01

    Autoimmune diseases affect up to approximately 10% of the population. While rare Mendelian autoimmunity syndromes can result from monogenic mutations disrupting essential mechanisms of central and peripheral tolerance, more common human autoimmune diseases are complex disorders that arise from the interaction between polygenic risk factors and environmental factors. Although the risk attributable to most individual nucleotide variants is modest, genome-wide association studies (GWAS) have the potential to provide an unbiased view of biological pathways that drive human autoimmune diseases. Interpretation of GWAS requires integration of multiple genomic datasets including dense genotyping, cis-regulatory maps of primary immune cells, and genotyped studies of gene expression in relevant cell types and cellular conditions. Improved understanding of the genetic basis of autoimmunity may lead to a more sophisticated understanding of underlying cellular phenotypes and, eventually, novel diagnostics and targeted therapies. PMID:26030227

  9. A coupled channel study of HN2 unimolecular decay based on a global ab initio potential surface

    NASA Technical Reports Server (NTRS)

    Koizumi, Hiroyasu; Schatz, George C.; Walch, Stephen P.

    1991-01-01

    The unimolecular decay lifetimes of several vibrational states of HN2 are determined on the basis of an accurate coupled channel dynamics study using a global analytical potential surface. The surface reproduces the ab initio points with an rms error of 0.08 kcal/mol for energies below 20 kcal/mol. Modifications to the potential that describe the effect of improving the basis set in the ab initio calculations are provided. Converged coupled channel calculations are performed for the ground rotational state of HN2 to determine the lifetimes of the lowest ten vibrational states. Only the ground vibrational state (000) and first excited bend (001) are found to have lifetimes longer than 1 ps. The lifetimes of these states are estimated at 3 x 10 to the -9th and 2 x 10 to the -10th s, respectively. Variation of these results with quality of the ab initio calculations is not more than a factor of 5.

  10. Structural and molecular basis for Ebola virus neutralization by protective human antibodies.

    PubMed

    Misasi, John; Gilman, Morgan S A; Kanekiyo, Masaru; Gui, Miao; Cagigi, Alberto; Mulangu, Sabue; Corti, Davide; Ledgerwood, Julie E; Lanzavecchia, Antonio; Cunningham, James; Muyembe-Tamfun, Jean Jacques; Baxa, Ulrich; Graham, Barney S; Xiang, Ye; Sullivan, Nancy J; McLellan, Jason S

    2016-03-18

    Ebola virus causes hemorrhagic fever with a high case fatality rate for which there is no approved therapy. Two human monoclonal antibodies, mAb100 and mAb114, in combination, protect nonhuman primates against all signs of Ebola virus disease, including viremia. Here, we demonstrate that mAb100 recognizes the base of the Ebola virus glycoprotein (GP) trimer, occludes access to the cathepsin-cleavage loop, and prevents the proteolytic cleavage of GP that is required for virus entry. We show that mAb114 interacts with the glycan cap and inner chalice of GP, remains associated after proteolytic removal of the glycan cap, and inhibits binding of cleaved GP to its receptor. These results define the basis of neutralization for two protective antibodies and may facilitate development of therapies and vaccines. PMID:26917592

  11. Ab initio calculations in a uniform magnetic field using periodic supercells

    SciTech Connect

    Cai, W; Galli, G

    2003-10-21

    We present a formulation of ab initio electronic structure calculations in a finite magnetic field, which retains the simplicity and efficiency of techniques widely used in first principles molecular dynamics simulations, based on plane-wave basis sets and Fourier transforms. In addition we discuss results obtained with this method for the energy spectrum of interacting electrons in quantum wells, and for the electronic properties of dense fluid deuterium in a uniform magnetic field.

  12. Characterization of the Cephalosporium acremonium pcbAB gene encoding alpha-aminoadipyl-cysteinyl-valine synthetase, a large multidomain peptide synthetase: linkage to the pcbC gene as a cluster of early cephalosporin biosynthetic genes and evidence of multiple functional domains.

    PubMed Central

    Gutiérrez, S; Díez, B; Montenegro, E; Martín, J F

    1991-01-01

    A 24-kb region of Cephalosporium acremonium C10 DNA was cloned by hybridization with the pcbAB and pcbC genes of Penicillium chrysogenum. A 3.2-kb BamHI fragment of this region complemented the mutation in the structural pcbC gene of the C. acremonium N2 mutant, resulting in cephalosporin production. A functional alpha-aminoadipyl-cysteinyl-valine (ACV) synthetase was encoded by a 15.6-kb EcoRI-BamHI DNA fragment, as shown by complementation of an ACV synthetase-deficient mutant of P. chrysogenum. Two transcripts of 1.15 and 11.4 kb were found by Northern (RNA blot) hybridization with probes internal to the pcbC and pcbAB genes, respectively. An open reading frame of 11,136 bp was located upstream of the pcbC gene that matched the 11.4-kb transcript initiation and termination regions. It encoded a protein of 3,712 amino acids with a deduced Mr of 414,791. The nucleotide sequence of the gene showed 62.9% similarity to the pcbAB gene encoding the ACV synthetase of P. chrysogenum; 54.9% of the amino acids were identical in both ACV synthetases. Three highly repetitive regions occur in the deduced amino acid sequence of C. acremonium ACV synthetase. Each is similar to the three repetitive domains in the deduced sequence of P. chrysogenum ACV synthetase and also to the amino acid sequence of gramicidin synthetase I and tyrocidine synthetase I of Bacillus brevis. These regions probably correspond to amino acid activating domains in the ACV synthetase protein. In addition, a thioesterase domain was present in the ACV synthetases of both fungi. A similarity has been found between the domains existing in multienzyme nonribosomal peptide synthetases and polyketide and fatty acid synthetases. The pcbAB gene is linked to the pcbC gene, forming a cluster of early cephalosporin-biosynthetic genes. Images PMID:1706706

  13. Accurate ab Initio Spin Densities

    PubMed Central

    2012-01-01

    We present an approach for the calculation of spin density distributions for molecules that require very large active spaces for a qualitatively correct description of their electronic structure. Our approach is based on the density-matrix renormalization group (DMRG) algorithm to calculate the spin density matrix elements as a basic quantity for the spatially resolved spin density distribution. The spin density matrix elements are directly determined from the second-quantized elementary operators optimized by the DMRG algorithm. As an analytic convergence criterion for the spin density distribution, we employ our recently developed sampling-reconstruction scheme [J. Chem. Phys.2011, 134, 224101] to build an accurate complete-active-space configuration-interaction (CASCI) wave function from the optimized matrix product states. The spin density matrix elements can then also be determined as an expectation value employing the reconstructed wave function expansion. Furthermore, the explicit reconstruction of a CASCI-type wave function provides insight into chemically interesting features of the molecule under study such as the distribution of α and β electrons in terms of Slater determinants, CI coefficients, and natural orbitals. The methodology is applied to an iron nitrosyl complex which we have identified as a challenging system for standard approaches [J. Chem. Theory Comput.2011, 7, 2740]. PMID:22707921

  14. On the Strain Rate Sensitivity of Abs and Abs Plus Fused Deposition Modeling Parts

    NASA Astrophysics Data System (ADS)

    Vairis, A.; Petousis, M.; Vidakis, N.; Savvakis, K.

    2016-06-01

    In this work the effect of strain rate on the tensile strength of fused deposition modeling parts built with Acrylonitrile-butadiene-styrene (ABS) and ABS plus material is presented. ASTM D638-02a specimens were built with ABS and ABS plus and they were tested on a Schenck Trebel Co. tensile test machine at three different test speeds, equal, lower, and higher to the test speed required by the ASTM D638-02a standard. The experimental tensile strength results were compared and evaluated. The fracture surfaces of selected specimens were examined with a scanning electron microscope, to determine failure mode of the filament strands. It was found that, as the test speed increases, specimens develop higher tensile strength and have higher elastic modulus. Specimens tested in the highest speed of the experiment had on average about 10% higher elastic modulus and developed on average about 11% higher tensile strength.

  15. Rapid digestion of Cry34Ab1 and Cry35Ab1 in simulated gastric fluid.

    PubMed

    Herman, Rod A; Schafer, Barry W; Korjagin, Valerie A; Ernest, April D

    2003-11-01

    Two genes were identified in Bacillus thuringiensis Berliner (Bt) that code for the proteins that comprise a Cry34Ab1/Cry35Ab1 binary insecticidal crystal protein. Maize, Zea mays L., plants have been transformed to express the Cry34Ab1/Cry35Ab1 proteins, and as a result, these plants are resistant to attack by western corn rootworm, Diabrotica virgifera virgifera LeConte, a major pest in the Midwestern corn-growing area of the U.S.A. As part of the safety assessment for the proteins, digestibility studies were conducted. Digestion experiments with both proteins demonstrated rapid degradation in simulated gastric fluid, comparable to other registered plant-incorporated protectants. Quantitative and qualitative approaches for determining digestibility are illustrated. PMID:14582981

  16. High-throughput mAb expression and purification platform based on transient CHO.

    PubMed

    Barnard, Gavin C; Hougland, Maria D; Rajendra, Yashas

    2015-01-01

    A high-cell-density transient transfection system was recently developed in our laboratory based on a CHO-GS-KO cell line. This method yields monoclonal antibody titers up to 350 mg/L from a simple 7-day process, in volumes ranging from 2 mL to 2 L. By performing transfections in 24-deep-well plates, a large number of mAbs can be expressed simultaneously. We coupled this new high-throughput transfection process to a semiautomated protein A purification process. Using a Biomek FX(p) liquid handling robot, up to 72 unique mAbs can be simultaneously purified. Our primary goal was to obtain >0.25 mg of purified mAb at a concentration of >0.5 mg/mL, without any concentration or buffer-exchange steps. We optimized both the batch-binding and the batch elution steps. The length of the batch-binding step was important to minimize mAb losses in the flowthrough fraction. The elution step proved to be challenging to simultaneously maximize protein recovery and protein concentration. We designed a variable volume elution strategy based on the average supernatant titer. Finally, we present two case studies. In the first study, we produced 56 affinity maturation mAb variants at an average yield of 0.33 ± 0.05 mg (average concentration of 0.65 ± 0.10 mg/mL). In a second study, we produced 42 unique mAbs, from an early-stage discovery effort, at an average yield of 0.79 ± 0.31 mg (average concentration of 1.59 ± 0.63 mg/mL). The combination of parallel high-yielding transient transfection and semiautomated high-throughput protein A purification represents a valuable mAb drug discovery tool. PMID:25403790

  17. Basis set dependence using DFT/B3LYP calculations to model the Raman spectrum of thymine.

    PubMed

    Bielecki, Jakub; Lipiec, Ewelina

    2016-02-01

    Raman spectroscopy (including surface enhanced Raman spectroscopy (SERS) and tip enhanced Raman spectroscopy (TERS)) is a highly promising experimental method for investigations of biomolecule damage induced by ionizing radiation. However, proper interpretation of changes in experimental spectra for complex systems is often difficult or impossible, thus Raman spectra calculations based on density functional theory (DFT) provide an invaluable tool as an additional layer of understanding of underlying processes. There are many works that address the problem of basis set dependence for energy and bond length consideration, nevertheless there is still lack of consistent research on basis set influence on Raman spectra intensities for biomolecules. This study fills this gap by investigating of the influence of basis set choice for the interpretation of Raman spectra of the thymine molecule calculated using the DFT/B3LYP framework and comparing these results with experimental spectra. Among 19 selected Pople's basis sets, the best agreement was achieved using 6-31[Formula: see text](d,p), 6-31[Formula: see text](d,p) and 6-11[Formula: see text]G(d,p) sets. Adding diffuse functions or polarized functions for small basis set or use of a medium or large basis set without diffuse or polarized functions is not sufficient to reproduce Raman intensities correctly. The introduction of the diffuse functions ([Formula: see text]) on hydrogen atoms is not necessary for gas phase calculations. This work serves as a benchmark for further research on the interaction of ionizing radiation with DNA molecules by means of ab initio calculations and Raman spectroscopy. Moreover, this work provides a set of new scaling factors for Raman spectra calculation in the framework of DFT/B3LYP method. PMID:26508426

  18. FTIR, Raman spectra and ab initio calculations of 2-mercaptobenzothiazole.

    PubMed

    Rai, Amareshwar K; Singh, Rachana; Singh, K N; Singh, V B

    2006-02-01

    FTIR and Raman spectra of a rubber vulcanization accelerator, 2-mercaptobenzothiazole (MBT), were recorded in the solid phase. The harmonic vibrational wavenumbers, for both the toutomeric forms of MBT, as well as for its dimeric complex, have been calculated, using ab initio RHF and density functional B3LYP methods invoking different basis sets upto RHF/6-31G** and B3LYP/6-31G** and the results were compared with the experimental values. Conformational studies have been also carried out regarding its toutomeric monomer forms and its dimer form. With all the basis sets the thione form of MBT (II) is predicted to be more stable than thiol form (I) and dimeric conformation (III) is predicted to be more stable with monomeric conformations (I) and (II). Vibrational assignments have been made, and it has been found that the calculated normal mode frequencies of dimeric conformation (III) are required for the analysis of IR and Raman bands of the MBT. The predicted shift in NH- stretching vibration towards the lower wave number side with the B3LYP/6-31G** calculations for the most stable dimer form (III), is in better agreement with experimental results. The intermolecular sulfur-nitrogen distance in N-H...S hydrogen bond was found to be 3.35 angstroms from these calculations, is also in agreement to the experimental value. PMID:16098806

  19. Ab Initio Simulation of the Photoelectron Spectrum for Methoxy Radical

    NASA Astrophysics Data System (ADS)

    Cheng, Lan; Weichman, Marissa L.; Kim, Jongjin B.; Ichino, Takatoshi; Neumark, Daniel; Stanton, John F.

    2015-06-01

    A theoretical simulation of the photoelectron spectrum for the ground state of methoxy radical is reported based on the quasidiabatic model Hamiltonian originally proposed by Köppel, Domcke, and Cederbaum. The parameters in the model Hamiltonian have been obtained from ab initio coupled-cluster calculations. The linear and quadratic force constants have been calculated using equation-of-motion coupled-cluster ionization potential method with the singles, doubles, and triples (EOMIP-CCSDT) truncation scheme together with atomic natural orbital basis sets of triple-zeta quality (ANO1). The cubic and quartic force constants have been obtained from EOMIP-CCSD calculations with ANO basis sets of double-zeta quality (ANO0), and the spin-orbit coupling constant has been computed at the EOMIP-CCSD/pCVTZ level. The nuclear Schroedinger equation has been solved using the Lanzcos algorithm to obtain vibronic energy levels as well as the corresponding intensities. The simulated spectrum compares favorably with the recent high-resolution slow electron velocity-map imaging experiment for vibronic levels up to 2000 cm-1.

  20. Ab initio study of the bonding in diatomic nickel

    SciTech Connect

    Noell, J.O.; Newton, M.D.; Hay, P.J.; Martin, R.L.; Bobrowicz, F.W.

    1980-09-01

    Hartree--Fock, GVB, and configuration interaction calculations were performed for diatomic nickel using an ab initio effective core potential. A basis set specifically optimized for the /sup 3/D state of atomic nickel is found to be far superior to the more common basis obtained from the /sup 3/F atomic state. Correlation effects are found to be significant in determining the bond energy. In particular, the two electrons of the s--s bond must be appropriately correlated. In addition, correlation effects which one would interpret as being principally intra-atomic in character are found to have a marked effect on the molecular properties. The theoretically predicted bond dissociation energy (D/sub e/) of 43.4 kcal/mol is significantly lower than the experimental estimate of 55 +- 5 kcal/mol. However, molecular partition functions calculated using the present results indicate that the experimental value should be revised downward to a value of approx.46 +- 5 kcal/mol, in good agreement with our calculations. An interatomic distance of 4.27 bohr is computed and compared with experimental estimates. Spectroscopic parameters for dipole-allowed transitions from the ground state were determined from SCF and GVB calculations and discussed in relation to the experimentally observed visible and ultraviolet spectra attributed to Ni/sub 2/.

  1. Accurate ab initio vibrational energies of methyl chloride.

    PubMed

    Owens, Alec; Yurchenko, Sergei N; Yachmenev, Andrey; Tennyson, Jonathan; Thiel, Walter

    2015-06-28

    Two new nine-dimensional potential energy surfaces (PESs) have been generated using high-level ab initio theory for the two main isotopologues of methyl chloride, CH3 (35)Cl and CH3 (37)Cl. The respective PESs, CBS-35( HL), and CBS-37( HL), are based on explicitly correlated coupled cluster calculations with extrapolation to the complete basis set (CBS) limit, and incorporate a range of higher-level (HL) additive energy corrections to account for core-valence electron correlation, higher-order coupled cluster terms, scalar relativistic effects, and diagonal Born-Oppenheimer corrections. Variational calculations of the vibrational energy levels were performed using the computer program TROVE, whose functionality has been extended to handle molecules of the form XY 3Z. Fully converged energies were obtained by means of a complete vibrational basis set extrapolation. The CBS-35( HL) and CBS-37( HL) PESs reproduce the fundamental term values with root-mean-square errors of 0.75 and 1.00 cm(-1), respectively. An analysis of the combined effect of the HL corrections and CBS extrapolation on the vibrational wavenumbers indicates that both are needed to compute accurate theoretical results for methyl chloride. We believe that it would be extremely challenging to go beyond the accuracy currently achieved for CH3Cl without empirical refinement of the respective PESs. PMID:26133427

  2. Accurate ab initio vibrational energies of methyl chloride

    NASA Astrophysics Data System (ADS)

    Owens, Alec; Yurchenko, Sergei N.; Yachmenev, Andrey; Tennyson, Jonathan; Thiel, Walter

    2015-06-01

    Two new nine-dimensional potential energy surfaces (PESs) have been generated using high-level ab initio theory for the two main isotopologues of methyl chloride, CH335Cl and CH337Cl. The respective PESs, CBS-35 HL, and CBS-37 HL, are based on explicitly correlated coupled cluster calculations with extrapolation to the complete basis set (CBS) limit, and incorporate a range of higher-level (HL) additive energy corrections to account for core-valence electron correlation, higher-order coupled cluster terms, scalar relativistic effects, and diagonal Born-Oppenheimer corrections. Variational calculations of the vibrational energy levels were performed using the computer program TROVE, whose functionality has been extended to handle molecules of the form XY 3Z. Fully converged energies were obtained by means of a complete vibrational basis set extrapolation. The CBS-35 HL and CBS-37 HL PESs reproduce the fundamental term values with root-mean-square errors of 0.75 and 1.00 cm-1, respectively. An analysis of the combined effect of the HL corrections and CBS extrapolation on the vibrational wavenumbers indicates that both are needed to compute accurate theoretical results for methyl chloride. We believe that it would be extremely challenging to go beyond the accuracy currently achieved for CH3Cl without empirical refinement of the respective PESs.

  3. Ab initio calculation of infrared intensities for hydrogen peroxide

    NASA Technical Reports Server (NTRS)

    Rogers, J. D.; Hillman, J. J.

    1982-01-01

    Results of an ab initio SCF quantum mechanical study are used to derive estimates for the infrared intensities of the fundamental vibrations of hydrogen peroxide. Atomic polar tensors (APTs) were calculated on the basis of a 4-31G basis set, and used to derive absolute intensities for the vibrational transitions. Comparison of the APTs calculated for H2O2 with those previously obtained for H2O and CH3OH, and of the absolute intensities derived from the H2O2 APTs with those derived from APTs transferred from H2O and CH3OH, reveals the sets of values to differ by no more than a factor of two, supporting the validity of the theoretical calculation. Values of the infrared intensities obtained correspond to A1 = 14.5 km/mol, A2 = 0.91 km/mol, A3 = 0.058 km/mol, A4 = 123 km/mol, A5 = 46.2 km/mol, and A6 = 101 km/mol. Charge, charge flux and overlap contributions to the dipole moment derivatives are also computed.

  4. Accurate ab initio vibrational energies of methyl chloride

    SciTech Connect

    Owens, Alec; Yurchenko, Sergei N.; Yachmenev, Andrey; Tennyson, Jonathan; Thiel, Walter

    2015-06-28

    Two new nine-dimensional potential energy surfaces (PESs) have been generated using high-level ab initio theory for the two main isotopologues of methyl chloride, CH{sub 3}{sup 35}Cl and CH{sub 3}{sup 37}Cl. The respective PESs, CBS-35{sup  HL}, and CBS-37{sup  HL}, are based on explicitly correlated coupled cluster calculations with extrapolation to the complete basis set (CBS) limit, and incorporate a range of higher-level (HL) additive energy corrections to account for core-valence electron correlation, higher-order coupled cluster terms, scalar relativistic effects, and diagonal Born-Oppenheimer corrections. Variational calculations of the vibrational energy levels were performed using the computer program TROVE, whose functionality has been extended to handle molecules of the form XY {sub 3}Z. Fully converged energies were obtained by means of a complete vibrational basis set extrapolation. The CBS-35{sup  HL} and CBS-37{sup  HL} PESs reproduce the fundamental term values with root-mean-square errors of 0.75 and 1.00 cm{sup −1}, respectively. An analysis of the combined effect of the HL corrections and CBS extrapolation on the vibrational wavenumbers indicates that both are needed to compute accurate theoretical results for methyl chloride. We believe that it would be extremely challenging to go beyond the accuracy currently achieved for CH{sub 3}Cl without empirical refinement of the respective PESs.

  5. Basis Selection for Wavelet Regression

    NASA Technical Reports Server (NTRS)

    Wheeler, Kevin R.; Lau, Sonie (Technical Monitor)

    1998-01-01

    A wavelet basis selection procedure is presented for wavelet regression. Both the basis and the threshold are selected using cross-validation. The method includes the capability of incorporating prior knowledge on the smoothness (or shape of the basis functions) into the basis selection procedure. The results of the method are demonstrated on sampled functions widely used in the wavelet regression literature. The results of the method are contrasted with other published methods.

  6. An Efficient Time-Stepping Scheme for Ab Initio Molecular Dynamics Simulations

    NASA Astrophysics Data System (ADS)

    Tsuchida, Eiji

    2016-08-01

    In ab initio molecular dynamics simulations of real-world problems, the simple Verlet method is still widely used for integrating the equations of motion, while more efficient algorithms are routinely used in classical molecular dynamics. We show that if the Verlet method is used in conjunction with pre- and postprocessing, the accuracy of the time integration is significantly improved with only a small computational overhead. We also propose several extensions of the algorithm required for use in ab initio molecular dynamics. The validity of the processed Verlet method is demonstrated in several examples including ab initio molecular dynamics simulations of liquid water. The structural properties obtained from the processed Verlet method are found to be sufficiently accurate even for large time steps close to the stability limit. This approach results in a 2× performance gain over the standard Verlet method for a given accuracy. We also show how to generate a canonical ensemble within this approach.

  7. Interaction-induced dipoles of hydrogen molecules colliding with helium atoms: a new ab initio dipole surface for high-temperature applications.

    PubMed

    Li, Xiaoping; Mandal, Anirban; Miliordos, Evangelos; Hunt, Katharine L C

    2012-01-28

    We report new ab initio results for the interaction-induced dipole moments Δμ of hydrogen molecules colliding with helium atoms. These results are needed in order to calculate collision-induced absorption spectra at high temperatures; applications include modeling the radiative profiles of very cool white dwarf stars, with temperatures from 3500 K to 9000 K. We have evaluated the dipoles based on finite-field calculations, with coupled cluster methods in MOLPRO 2006 and aug-cc-pV5Z (spdfg) basis sets for both the H and He centers. We have obtained values of Δμ for eight H(2) bond lengths ranging from 0.942 a.u. to 2.801 a.u., for 15 intermolecular separations R ranging from 2.0 a.u. to 10.0 a.u., and for 19 different relative orientations. In general, our values agree well with earlier ab initio results, for the geometrical configurations that are treated in common, but we have determined more points on the collision-induced dipole surface by an order of magnitude. These results make it possible to calculate transition probabilities for molecules in excited vibrational states, overtones, and rotational transitions with ΔJ > 4. We have cast our results in the symmetry-adapted form needed for absorption line shape calculations, by expressing Δμ as a series in the spherical harmonics of the orientation angles of the intermolecular vector and of a unit vector along the H(2) bond axis. The expansion coefficients depend on the H(2) bond length and the intermolecular distance R. For large separations R, we show that the ab initio values of the leading coefficients converge to the predictions from perturbation theory, including both classical multipole polarization and dispersion effects. PMID:22299884

  8. Accurate Bond Energies of Hydrocarbons from Complete Basis Set Extrapolated Multi-Reference Singles and Doubles Configuration Interaction

    SciTech Connect

    Oyeyemi, Victor B.; Pavone, Michele; Carter, Emily A.

    2011-11-03

    Quantum chemistry has become one of the most reliable tools for characterizing the thermochemical underpinnings of reactions, such as bond dissociation energies (BDEs). The accurate prediction of these particular properties (BDEs) are challenging for ab initio methods based on perturbative corrections or coupled cluster expansions of the single-determinant Hartree-Fock wave function: the processes of bond breaking and forming are inherently multi-configurational and require an accurate description of non-dynamical electron correlation. To this end, we present a systematic ab initio approach for computing BDEs that is based on three components: (1) multi-reference single and double excitation configuration interaction (MRSDCI) for the electronic energies; (2) a two-parameter scheme for extrapolating MRSDCI energies to the complete basis set limit; and (3) DFT-B3LYP calculations of minimumenergy structures and vibrational frequencies to account for zero point energy and thermal corrections. We validated our methodology against a set of reliable experimental BDE values of C*C and C*H bonds of hydrocarbons. The goal of chemical accuracy is achieved, on average, without applying any empirical corrections to the MRSDCI electronic energies. We then use this composite scheme to make predictions of BDEs in a large number of hydrocarbon molecules for which there are no experimental data, so as to provide needed thermochemical estimates for fuel molecules.

  9. Collective rotation from ab initio theory

    NASA Astrophysics Data System (ADS)

    Caprio, Mark A.; Maris, Pieter; Vary, James P.

    2015-10-01

    The challenge of ab initio nuclear theory is to quantitatively predict the complex and highly-correlated behavior of the nuclear many-body system, starting from the underlying internucleon interactions. We may now seek to understand the wealth of nuclear collective phenomena through ab initio approaches. No-core configuration interaction (NCCI) calculations for p-shell nuclei give rise to rotational bands, as evidenced by rotational patterns for excitation energies, electromagnetic moments, and electromagnetic transitions. In this talk, the intrinsic structure of these bands is discussed, and the predicted rotational bands are compared to experiment. Supported by the US DOE under Award Nos. DE-FG02-95ER-40934, DESC0008485 (SciDAC/NUCLEI), and DE-FG02-87ER40371 and the US NSF under Award No. 0904782. Computational resources provided by NERSC (US DOE Contract No. DE-AC02-05CH11231).

  10. Ab initio infrared and Raman spectra

    NASA Technical Reports Server (NTRS)

    Fredkin, D. R.; White, S. R.; Wilson, K. R.; Komornicki, A.

    1983-01-01

    It is pointed out that with increased computer power and improved computational techniques, such as the gradients developed in recent years, it is becoming practical to compute spectra ab initio, from the fundamental constants of nature, for systems of increasing complexity. The present investigation has the objective to explore several possible ab initio approaches to spectra, giving particular attention to infrared and nonresonance Raman. Two approaches are discussed. The sequential approach, in which first the electronic part and then later the nuclear part of the Born-Oppenheimer approximation is solved, is appropriate for small systems. The simultaneous approach, in which the electronic and nuclear parts are solved at the same time, is more appropriate for many-atom systems. A review of the newer quantum gradient techniques is provided, and the infrared and Raman spectral band contours for the water molecule are computed.

  11. Ab initio vel ex eventu. II

    NASA Astrophysics Data System (ADS)

    Thiessen, P. A.; Treder, H.-J.

    Jedes initium wird durch experimenta crucis zum eventus. Jedes theoretisch interpretierbare ex-eventu-Resultat führt auf ein neues Initium. Gerade dies ist die gemeinsame Aussage von Atomistik, Quantenmechanik und Relativitätstheorie.Translated AbstractAb initio vel ex eventu. IIEvery initium becomes an eventus by experimenta crucis. Every theoretically interpretable ex-eventu result leads to a new initium. Right this is the joint assertion of atomism, quantum mechanics, and relativity.

  12. Interaction of vermiculite with Katamin AB

    SciTech Connect

    Veksler, V.I.; Makashev, Yu.A.; Mikhailova, V.N.; Aleksandrova, E.M.; Vinogradova, I.A.

    1986-12-10

    Vermiculite can interact with quaternary ammonium salts to form organmineral compounds in the ion-exchange adsorption of these salts from waste water. This paper investigates this interaction between vermiculite and Katamin AB, a commercial cationic surfactant consisting of an aqueous solution of alkylbenyldimethylammonium chloride homologs. It is shown that ion exchange and molecular sorption occur simultaneously during the interaction and that vermiculite is a highly effective adsorbent in process water purification.

  13. Detection and analysis of tupaia hepatocytes via mAbs against tupaia serum albumin

    PubMed Central

    Liu, Xuan; Yuan, Lunzhi; Yuan, Quan; Zhang, Yali; Wu, Kun; Zhang, Tianying; Wu, Yong; Hou, Wangheng; Wang, Tengyun; Liu, Pingguo; Shih, James Wai Kuo; Cheng, Tong; Xia, Ningshao

    2015-01-01

    On the basis of its close phylogenetic relationship with primates, the development of Tupaia belangeri as an infection animal model and drug metabolism model could provide a new option for preclinical studies, especially in hepatitis virus research. As a replacement for primary human hepatocytes (PHHs), primary tupaia hepatocytes (PTHs) have been widely used. Similar to human serum albumin, tupaia serum albumin (TSA) is the most common liver synthesis protein and is an important biomarker for PTHs and liver function. However, no detection or quantitative method for TSA has been reported. In this study, mouse monoclonal antibodies (mAbs) 4G5 and 9H3 against TSA were developed to recognize PTHs, and they did not show cross-reactivity with serum albumin from common experimental animals, such as the mouse, rat, cow, rabbit, goat, monkey, and chicken. The two mAbs also exhibited good performance in fluorescence activated cell sorting (FACS) analysis and immunofluorescence (IF) detection of PTHs. A chemiluminescent enzyme immune assay method using the two mAbs, with a linear range from 96.89 pg/ml to 49,609.38 pg/ml, was developed for the quantitative detection of TSA. The mAbs and the CLEIA method provide useful tools for research on TSA and PTHs. PMID:26597317

  14. Detection and analysis of tupaia hepatocytes via mAbs against tupaia serum albumin.

    PubMed

    Liu, Xuan; Yuan, Lunzhi; Yuan, Quan; Zhang, Yali; Wu, Kun; Zhang, Tianying; Wu, Yong; Hou, Wangheng; Wang, Tengyun; Liu, Pingguo; Shih, James Wai Kuo; Cheng, Tong; Xia, Ningshao

    2016-05-20

    On the basis of its close phylogenetic relationship with primates, the development of Tupaia belangeri as an infection animal model and drug metabolism model could provide a new option for preclinical studies, especially in hepatitis virus research. As a replacement for primary human hepatocytes (PHHs), primary tupaia hepatocytes (PTHs) have been widely used. Similar to human serum albumin, tupaia serum albumin (TSA) is the most common liver synthesis protein and is an important biomarker for PTHs and liver function. However, no detection or quantitative method for TSA has been reported. In this study, mouse monoclonal antibodies (mAbs) 4G5 and 9H3 against TSA were developed to recognize PTHs, and they did not show cross-reactivity with serum albumin from common experimental animals, such as the mouse, rat, cow, rabbit, goat, monkey, and chicken. The two mAbs also exhibited good performance in fluorescence activated cell sorting (FACS) analysis and immunofluorescence (IF) detection of PTHs. A chemiluminescent enzyme immune assay method using the two mAbs, with a linear range from 96.89 pg/ml to 49,609.38 pg/ml, was developed for the quantitative detection of TSA. The mAbs and the CLEIA method provide useful tools for research on TSA and PTHs. PMID:26597317

  15. CVRQD ab initio ground-state adiabatic potential energy surfaces for the water molecule.

    PubMed

    Barletta, Paolo; Shirin, Sergei V; Zobov, Nikolai F; Polyansky, Oleg L; Tennyson, Jonathan; Valeev, Edward F; Császár, Attila G

    2006-11-28

    The high accuracy ab initio adiabatic potential energy surfaces (PESs) of the ground electronic state of the water molecule, determined originally by Polyansky et al. [Science 299, 539 (2003)] and called CVRQD, are extended and carefully characterized and analyzed. The CVRQD potential energy surfaces are obtained from extrapolation to the complete basis set of nearly full configuration interaction valence-only electronic structure computations, augmented by core, relativistic, quantum electrodynamics, and diagonal Born-Oppenheimer corrections. We also report ab initio calculations of several quantities characterizing the CVRQD PESs, including equilibrium and vibrationally averaged (0 K) structures, harmonic and anharmonic force fields, harmonic vibrational frequencies, vibrational fundamentals, and zero-point energies. They can be considered as the best ab initio estimates of these quantities available today. Results of first-principles computations on the rovibrational energy levels of several isotopologues of the water molecule are also presented, based on the CVRQD PESs and the use of variational nuclear motion calculations employing an exact kinetic energy operator given in orthogonal internal coordinates. The variational nuclear motion calculations also include a simplified treatment of nonadiabatic effects. This sophisticated procedure to compute rovibrational energy levels reproduces all the known rovibrational levels of the water isotopologues considered, H(2) (16)O, H(2) (17)O, H(2) (18)O, and D(2) (16)O, to better than 1 cm(-1) on average. Finally, prospects for further improvement of the ground-state adiabatic ab initio PESs of water are discussed. PMID:17144700

  16. Ab interno trabeculectomy: patient selection and perspectives

    PubMed Central

    Vinod, Kateki; Gedde, Steven J

    2016-01-01

    Ab interno trabeculectomy is one among several recently introduced minimally invasive glaucoma surgeries that avoid a conjunctival incision and full-thickness sclerostomy involved in traditional glaucoma surgery. Ablation of the trabecular meshwork and inner wall of Schlemm’s canal is performed in an arcuate fashion via a clear corneal incision, alone or in combination with phacoemulsification cataract surgery. Intraocular pressure reduction following ab interno trabeculectomy is limited by resistance in distal outflow pathways and generally stabilizes in the mid-to-high teens. Relief of medication burden has been demonstrated by some studies. A very low rate of complications, most commonly transient hyphema and intraocular pressure elevations in the immediate postoperative period, have been reported. However, available data are derived from small retrospective and prospective case series. Randomized, controlled trials are needed to better elucidate the potential merits of ab interno trabeculectomy in the combined setting versus phacoemulsification cataract surgery alone and to compare it with other minimally invasive glaucoma surgeries. PMID:27574396

  17. Ab-inition melting curve of titanium

    NASA Astrophysics Data System (ADS)

    Stutzmann, Vincent; Bouchet, Johann; Bottin, Francois

    2014-03-01

    Thermodynamical properties of titanium are of great interest for aerospace and aviation industries and many studies are done in order to understand its behaviour under pressure (P) and temperature (T) : phase transitions at low T, melting curve at high T and P. In this work we compute the first ab-initio melting curve of titanium. This one is obtained with the Abinit package using DFT, in the GGA approximation, and in the framework of the projector augmented wave method (PAW). At first, we perform ground state calculations and study the five allotropic phases of titanium. Two PAW atomic data are generated with two different cutoff radius. The larger one gives results near previews ab-initio calculations, whereas the smaller one gives results near all electron calculation. Using the second PAW atomic data and performing ab-initio molecular dynamic simulations, we then compute the melting curve of titanium with three different methods. Results show relevance of our calculations, but also discrepencies with experimental data.

  18. Ab initio Bogoliubov coupled cluster theory

    NASA Astrophysics Data System (ADS)

    Signoracci, Angelo; Hagen, Gaute; Duguet, Thomas

    2014-09-01

    Coupled cluster (CC) theory has become a standard method in nuclear theory for realistic ab initio calculations of medium mass nuclei, but remains limited by its requirement of a Slater determinant reference state which reasonably approximates the nuclear system of interest. Extensions of the method, such as equation-of-motion CC, permit the calculation of nuclei with one or two nucleons added or removed from a doubly magic core, yet still only a few dozen nuclei are accessible with modern computational restrictions. In order to extend the applicability of ab initio methods to open-shell systems, the superfluid nature of nuclei must be taken into account. By utilizing Bogoliubov algebra and employing spontaneous symmetry breaking with respect to particle number conservation, superfluid systems can be treated by a single reference state. An ab initio theory to include correlations on top of a Bogoliubov reference state has been developed in the guise of standard CC theory. The formalism and first results of this Bogoliubov coupled cluster theory will be presented to demonstrate the applicability of the method.

  19. Genetic basis of cohesinopathies

    PubMed Central

    Barbero, José L

    2013-01-01

    Cohesin is a ring-form multifunctional protein complex, which was discovered during a search for molecules that keep sister chromatids together during segregation of chromosomes during cell division. In the past decade, a large number of results have also demonstrated a need for the cohesin complex in other crucial events in the life cycle of the cell, including DNA duplication, heterochromatin formation, DNA double-strand break repair, and control of gene expression. The dynamics of the cohesin ring are modulated by a number of accessory and regulatory proteins, known as cohesin cofactors. Loss of function of the cohesin complex is incompatible with life; however, mutations in the genes encoding for cohesin subunits and/or cohesin cofactors, which have very little or a null effect on chromosome segregation, represent a newly recognized class of human genetic disorders known as cohesinopathies. A number of genetic, biochemical, and clinical approaches, and importantly, animal models, can help us to determine the underlying mechanisms for these human diseases. PMID:23882154

  20. Ab initio ground and the first excited adiabatic and quasidiabatic potential energy surfaces of H + + CO system

    NASA Astrophysics Data System (ADS)

    George, D. X. F.; Kumar, Sanjay

    2010-08-01

    Ab initio global adiabatic as well as quasidiabatic potential energy surfaces for the ground and the first excited electronic states of the H + + CO system have been computed as a function of the Jacobi coordinates ( R, r, γ) using Dunning's cc-pVTZ basis set at the internally contracted multi-reference (single and double) configuration interaction level of accuracy. In addition, nonadiabatic coupling matrix elements arising from radial motion, mixing angle and coupling potential have been computed using the ab initio procedure [Simah et al. (1999) [66

  1. Highly anisotropic thermal conductivity of arsenene: An ab initio study

    NASA Astrophysics Data System (ADS)

    Zeraati, Majid; Vaez Allaei, S. Mehdi; Abdolhosseini Sarsari, I.; Pourfath, Mahdi; Donadio, Davide

    2016-02-01

    Elemental two-dimensional (2D) materials exhibit intriguing heat transport and phononic properties. Here we have investigated the lattice thermal conductivity of newly proposed arsenene, the 2D honeycomb structure of arsenic, using ab initio calculations. Solving the Boltzmann transport equation for phonons, we predict a highly anisotropic thermal conductivity of 30.4 and 7.8 W/mK along the zigzag and armchair directions, respectively, at room temperature. Our calculations reveal that phonons with mean free paths between 20 nm and 1 μ m provide the main contribution to the large thermal conductivity in the zigzag direction; mean free paths of phonons contributing to heat transport in the armchair directions range between 20 and 100 nm. The obtained anisotropic thermal conductivity and feasibility of synthesis, in addition to high electron mobility reported elsewhere, make arsenene a promising material for nanoelectronic applications and thermal management.

  2. SEEDS Polarimetric Imagery of the AB Aur Protoplanetary Disk

    NASA Astrophysics Data System (ADS)

    Wisniewski, John P.; Fukagawa, M.; Grady, C.; Hashimoto, J.; Hodapp, K.; Kudo, T.; Munetake, M.; Okamoto, Y.; Tamura, M.; SEEDS Team

    2011-01-01

    The Strategic Exploration of Exoplanets and Disks with Subaru (SEEDS) is a large survey which will be observing roughly 200 protoplanetary and debris disk systems over the next five years using the HiCIAO coronagraph + AO188 system on the Subaru telescope. We present new J-band polarimetric differential imagery of the proto-type Herbig Ae star, AB Aurigae, which diagnoses scattered light from the system between 20 - 540 AU at a resolution of roughly 8 AU. We discuss the morphology we observe in the outer disk region in the context of previous observations of the system, and compare/contrast the morphology in the inner disk region with recent H-band imagery of the system made with HiCIAO (Hashimoto et al 2010). This work was supported in part by NSF grants AST 0802230 and AST 1009314 and the AAS' Chretien International Research Grant.

  3. Ab initio calculations of grain boundaries in bcc metals

    NASA Astrophysics Data System (ADS)

    Scheiber, Daniel; Pippan, Reinhard; Puschnig, Peter; Romaner, Lorenz

    2016-03-01

    In this study, we compute grain boundary (GB) properties for a large set of GBs in bcc transition metals with a special focus on W, Mo and Fe using ab initio density functional theory (DFT) and semi-empirical second nearest neighbour modified embedded atom method (2NN-MEAM) potentials. The GB properties include GB energies, surface energies, GB excess volume and work of separation, which we analyse and then compare to experimental data. We find that the used 2NN-MEAM potentials can predict general trends of GB properties, but do not always reproduce the GB ground state structure and energy found with DFT. In particular, our results explain the experimental finding that W and Mo prefer intergranular fracture, while other bcc metals prefer transgranular cleavage.

  4. An improved choice of oscillator basis for banana shaped nuclides

    SciTech Connect

    Chasman, R.R.

    1994-03-01

    The question of the appropriate choice of oscillator basis functions for studying exotic nuclear shapes is raised. Difficulties with the conventional choice of oscillator basis states are noted for shapes having a large banana component. A prescription for an improved oscillator basis to study these shapes is given. It can be applied in a more general context. New calculations with this improved basis are presented for the banana deformation mode. The change of basis gives results that improve the prospects of finding states in the banana minimum for many isotopes of Tl, Pb and Bi.

  5. Anharmonic force field and vibrational dynamics of CH2F2 up to 5000 cm(-1) studied by Fourier transform infrared spectroscopy and state-of-the-art ab initio calculations.

    PubMed

    Tasinato, Nicola; Regini, Giorgia; Stoppa, Paolo; Pietropolli Charmet, Andrea; Gambi, Alberto

    2012-06-01

    Difluoromethane (CH(2)F(2), HFC-32) is a molecule used in refrigerant mixtures as a replacement of the more environmentally hazardous, ozone depleting, chlorofluorocarbons. On the other hand, presenting strong vibration-rotation bands in the 9 μm atmospheric window, it is a greenhouse gas which contributes to global warming. In the present work, the vibrational and ro-vibrational properties of CH(2)F(2), providing basic data for its atmospheric modeling, are studied in detail by coupling medium resolution Fourier transform infrared spectroscopy to high-level electronic structure ab initio calculations. Experimentally a full quantum assignment and accurate integrated absorption cross sections are obtained up to 5000 cm(-1). Ab initio calculations are carried out by using CCSD(T) theory and large basis sets of either the correlation consistent or atomic natural orbital hierarchies. By using vibrational perturbation theory to second order a complete set of vibrational and ro-vibrational parameters is derived from the ab initio quartic anharmonic force fields, which well compares with the spectroscopic constants retrieved experimentally. An excellent agreement between theory and experiment is achieved for vibrational energy levels and integrated absorption cross sections: transition frequencies up to four quanta of vibrational excitation are reproduced with a root mean square deviation (RMSD) of 7 cm(-1) while intensities are predicted within few km mol(-1) from the experiment. Basis set performances and core correlation effects are discussed throughout the paper. Particular attention is focused in the understanding of the anharmonic couplings which rule the vibrational dynamics of the |ν(1)>, |2ν(8)>, |2ν(2)> three levels interacting system. The reliability of the potential energy and dipole moment surfaces in reproducing the vibrational eigenvalues and intensities as well as in modeling the vibrational and ro-vibrational mixings over the whole 400-5000 cm(-1

  6. Undoing Gender Through Legislation and Schooling: the Case of AB 537 and AB 394 IN California, USA

    NASA Astrophysics Data System (ADS)

    Knotts, Greg

    2009-11-01

    This article investigates California laws AB 537: The Student Safety and Violence Prevention Act of 2000, and the recently enacted AB 394: Safe Place to Learn Act. Both demand that gender identity and sexual orientation be added to the lexicon of anti-harassment protection in public education. However, despite these progressive measures, schools have an unconscious acceptance of heteronormativity and gendered norms, which undermines both the spirit and language of these laws. This paper examines how California schools can both change standard practices and realise the transformative social change that laws like AB 537 and AB 394 can instigate. I assert that the systemic implementation of these laws, through the adoption, enforcement and evaluation of existing AB 537 Task Force Recommendations, is necessary for their success. My second assertion is that AB 537 and AB 394 have the potential to change and reconstitute gender-based and heteronormative standards at school sites.

  7. Molecular Basis of Differential B-Pentamer Stability of Shiga Toxins 1 and 2

    SciTech Connect

    Conrady, Deborah G.; Flagler, Michael J.; Friedmann, David R.; Vander Wielen, Bradley D.; Kovall, Rhett A.; Weiss, Alison A.; Herr, Andrew B.

    2012-06-27

    Escherichia coli strain O157:H7 is a major cause of food poisoning that can result in severe diarrhea and, in some cases, renal failure. The pathogenesis of E. coli O157:H7 is in large part due to the production of Shiga toxin (Stx), an AB{sub 5} toxin that consists of a ribosomal RNA-cleaving A-subunit surrounded by a pentamer of receptor-binding B subunits. There are two major isoforms, Stx1 and Stx2, which differ dramatically in potency despite having 57% sequence identity. Animal studies and epidemiological studies show Stx2 is associated with more severe disease. Although the molecular basis of this difference is unknown, data suggest it is associated with the B-subunit. Mass spectrometry studies have suggested differential B-pentamer stability between Stx1 and Stx2. We have examined the relative stability of the B-pentamers in solution. Analytical ultracentrifugation using purified B-subunits demonstrates that Stx2B, the more deadly isoform, shows decreased pentamer stability compared to Stx1B (EC{sub 50} = 2.3 {micro}M vs. EC{sub 50} = 0.043 {micro}M for Stx1B). X-ray crystal structures of Stx1B and Stx2B identified a glutamine in Stx2 (versus leucine in Stx1) within the otherwise strongly hydrophobic interface between B-subunits. Interchanging these residues switches the stability phenotype of the B-pentamers of Stx1 and Stx2, as demonstrated by analytical ultracentrifugation and circular dichroism. These studies demonstrate a profound difference in stability of the B-pentamers in Stx1 and Stx2, illustrate the mechanistic basis for this differential stability, and provide novel reagents to test the basis for differential pathogenicity of these toxins.

  8. Ab initio charge-carrier mobility model for amorphous molecular semiconductors

    NASA Astrophysics Data System (ADS)

    Massé, Andrea; Friederich, Pascal; Symalla, Franz; Liu, Feilong; Nitsche, Robert; Coehoorn, Reinder; Wenzel, Wolfgang; Bobbert, Peter A.

    2016-05-01

    Accurate charge-carrier mobility models of amorphous organic molecular semiconductors are essential to describe the electrical properties of devices based on these materials. The disordered nature of these semiconductors leads to percolative charge transport with a large characteristic length scale, posing a challenge to the development of such models from ab initio simulations. Here, we develop an ab initio mobility model using a four-step procedure. First, the amorphous morphology together with its energy disorder and intermolecular charge-transfer integrals are obtained from ab initio simulations in a small box. Next, the ab initio information is used to set up a stochastic model for the morphology and transfer integrals. This stochastic model is then employed to generate a large simulation box with modeled morphology and transfer integrals, which can fully capture the percolative charge transport. Finally, the charge-carrier mobility in this simulation box is calculated by solving a master equation, yielding a mobility function depending on temperature, carrier concentration, and electric field. We demonstrate the procedure for hole transport in two important molecular semiconductors, α -NPD and TCTA. In contrast to a previous study, we conclude that spatial correlations in the energy disorder are unimportant for α -NPD. We apply our mobility model to two types of hole-only α -NPD devices and find that the experimental temperature-dependent current density-voltage characteristics of all devices can be well described by only slightly decreasing the simulated energy disorder strength.

  9. 40 CFR 174.506 - Bacillus thuringiensis Cry34Ab1 and Cry35Ab1 proteins in corn; exemption from the requirement of...

    Code of Federal Regulations, 2013 CFR

    2013-07-01

    ... Cry35Ab1 proteins in corn; exemption from the requirement of a tolerance. 174.506 Section 174.506... thuringiensis Cry34Ab1 and Cry35Ab1 proteins in corn; exemption from the requirement of a tolerance. Residues of Bacillus thuringiensis Cry34Ab1 and Cry35Ab1 proteins in corn are exempted from the requirement of...

  10. 40 CFR 174.506 - Bacillus thuringiensis Cry34Ab1 and Cry35Ab1 proteins in corn; exemption from the requirement of...

    Code of Federal Regulations, 2012 CFR

    2012-07-01

    ... Cry35Ab1 proteins in corn; exemption from the requirement of a tolerance. 174.506 Section 174.506... thuringiensis Cry34Ab1 and Cry35Ab1 proteins in corn; exemption from the requirement of a tolerance. Residues of Bacillus thuringiensis Cry34Ab1 and Cry35Ab1 proteins in corn are exempted from the requirement of...