Structure optimization combining soft-core interaction functions, the diffusion equation method, and molecular dynamics

Smoothing the potential energy surface for structure optimization is a general and commonly applied strategy. We propose a combination of soft-core potential energy functions and a variation of the diffusion equation method to smooth potential energy surfaces, which is applicable to complex systems such as protein structures; The performance of the method was demonstrated by comparison with simulated annealing using the refinement of the undecapeptide Cyclosporin A as a test case. Simulations were repeated many times using different initial conditions and structures since the methods are heuristic and results are only meaningful in a statistical sense.

Rates of electronic energy transfer in conformationally flexible bichromophoric macrocyclic complexes: a combined experimental and molecular modeling study

Electronic energy transfer (EET) rate constants between a naphthalene donor and anthracene acceptor in [ZnL4a](ClO4)(2) and [ZnL4b](ClO4)(2) were determined by time-resolved fluorescence where L-4a and L-4b are the trans and cis isomers of 6-((anthracen-9-yl-methyl)amino)-6,13-dimethyl-13-((naphthalen-1-yl-methyl)amino)-1,4,8,11-tetraazacyclotetradecane, respectively. These isomers differ in the relative disposition of the appended chromophores with respect to the macrocyclic plane. The trans isomer has an energy transfer rate constant (k(EET)) of 8.7 x 10(8) s(-1), whereas that of the cis isomer is significantly faster (2.3 x 10(9) s(-1)). Molecular modeling was used to determine the likely distribution of conformations in CH3CN solution for these complexes in an attempt to identify any distance or orientation dependency that may account for the differing rate constants observed. The calculated conformational distributions together with analysis by H-1 NMR for the [ZnL4a](2+) trans complex in the common trans-III N-based isomer gave a calculated Forster rate constant close to that observed experimentally. For the [ZnL4b](2+) cis complex, the experimentally determined rate constant may be attributed to a combination of trans-Ill and trans-I N-based isomeric forms of the complex in solution.

Atomic diversity, molecular diversity, and chemical diversity: the concept of chemodiversity.

This minireview is meant as an introduction to the following paper. To this end, it presents the general background against which the joint paper should be understood. The first objective of the present paper is thus to clarify some concepts and related terminology, drawing a clear distinction between i) atomic diversity (i.e., atomic-property space), ii) molecular or macromolecular diversity (i.e., molecular- or macromolecular-property spaces), and iii) chemical diversity (i.e., chemical-diversity space). The first refers to the various electronic states an atom can occupy. The second encompasses the conformational and property spaces of a given (macro)molecule. The third pertains to the diversity in structure and properties exhibited by a library or a supramolecular assembly of different chemical compounds. The ground is thus laid for the content of the joint paper, which pertains to case ii, to be placed in its broader chemodiversity context. The second objective of this paper is to point to the concepts of chemodiversity and biodiversity as forming a continuum. Chemodiversity is indeed the material substratum of organisms. In other words, chemodiversity is the material condition for life to emerge and exist. Increasing our knowledge of chemodiversity is thus a condition for a better understanding of life as a process.

Molecular Physics of Elementary Processes relevant to Hypersonics: atom-molecule, molecule-molecule and atom-surface processes.

In the present chapter some prototype gas and gas-surface processes occurring within the hypersonic flow layer surrounding spacecrafts at planetary entry are discussed. The discussion is based on microscopic dynamical calculations of the detailed cross sections and rate coefficients performed using classical mechanics treatments for atoms, molecules and surfaces. Such treatment allows the evaluation of the efficiency of thermal processes (both at equilibrium and nonequilibrium distributions) based on state-to-state and state specific calculations properly averaged over the population of the initial states. The dependence of the efficiency of the considered processes on the initial partitioning of energy among the various degrees of freedom is discussed.

Relativistic effects in physics and chemistry of element 105. [Part] IV.

Results of relativistic (Dirac-Slater and Dirac-Fock) and nonrelativistic (Hartree-Fock-Slater) atomic and molecular calculations have been compared for the group 5 elements Nb, Ta, and Ha and their compounds MCl_5, to elucidate the influence of relativistic effects on their properties especially in going from the 5d element Ta to the 6d element Ha. The analysis of the radial distribution of the valence electrons of the metals for electronic configurations obtained as a result of the molecular calculations and their overlap with ligands show opposite trends in behavior for ns_1/2, np_l/2, and (n -1 )d_5/2 orbitals for Ta and Ha in the relativistic and nonrelativistic cases. Relativistic contraction and energetic stabilization of the ns_1/2 and np_l/2 wave functions and expansion and destabilization of the (n-1)d_5/2 orbitals make hahnium pentahalide more covalent than tantalum pentahalide and increase the bond strength. The nonrelativistic treatment of the wave functions results in an increase in ionicity of the MCl_5 molecules in going from Nb to Ha making element Ha an analog of V. Different trends for the relativistic and nonrelativistic cases are also found for ionization potentials, electronic affinities, and energies of charge-transfer transitions as well as the stability of the maximum oxidation state.

A novel set of ss-N-biaryl ether sulfonamide hydroxamates as potential MMPs inhibitors: Molecular dynamics simulations and molecular properties evaluation

Matrix metalloproteinases (MMPs) constitute a family of zinc-dependent proteases involved in the extracellular matrix degradation. MMP-2 and MMP9 are overexpressed in several human cancer types, including melanoma, thus the development of new compounds to inhibit MMPs' activity is desirable. Molecular dynamic simulation and molecular properties calculations were performed on a set of novel beta-N-biaryl ether sulfonamide-based hydroxamates, reported as MMP-2 and MMP-9 inhibitors, for providing data to develop an exploratory analysis. Thermodynamic, electronic, and steric descriptors have significantly discriminated highly active from moderately and less active inhibitors of MMP-2 whereas apparent partition coefficient at pH 1.5 was also significant for the MMP-9 data set. Compound 47 was considered an outlier in all analysis, indicating the presence of a bulky substituent group in R3 is crucial to this set of inhibitors for the establishment of molecular interactions with the S1 subsite of both enzymes, but there is a limit. (C) 2012 Wiley Periodicals, Inc.

Theoretical investigation of atomic and electronic structures of Ga(2)O(3)(ZnO)(6)

Transparent conducting oxides (TCO) are widely used in technological applications ranging from photovoltaics to thin-film transparent field-effect transistors. In this work we report a first-principles investigation, based on density-functional theory, of the atomic and electronic properties of Ga(2)O(3)(ZnO)(6) (GZO(6)), which is a promising candidate to be used as host oxide for wide band gap TCO applications. We identify a low-energy configuration for the coherent distribution of the Ga and Zn atoms in the cation positions within the experimentally reported orthorhombic GZO(6) structure. Four Ga atoms are located in four-fold sites, while the remaining 12 Ga atoms in the unit cell form four shared Ga agglomerates (a motif of four atoms). The Zn atoms are distributed in the remaining cation sites with effective coordination numbers from 3.90 to 4.50. Furthermore, we identify the natural formation of twin-boundaries in GZO(6), which can explain the zigzag modulations observed experimentally by high-resolution transmission electron microscopy in GZO(n) (n=9). Due to the intrinsic twin-boundary formation, polarity inversion in the ZnO tetrahedrons is present which is facilitated by the formation of the Ga agglomerates. Our analysis shows that the formation of fourfold Ga sites and Ga agglomerates are stabilized by the electronic octet rule, while the distribution of Ga atoms and the formation of the twin-boundary help alleviate excess strain. Finally we identify that the electronic properties of GZO(6) are essentially determined by the electronic properties of ZnO, i.e., there are slight changes in the band gap and optical absorption properties.

Novel solitons in parametric amplifiers and atom lasers

We review recent developments in quantum and classical soliton theory, leading to the possibility of observing both classical and quantum parametric solitons in higher-dimensional environments. In particular, we consider the theory of three bosonic fields interacting via both parametric (cubic) and quartic couplings. In the case of photonic fields in a nonlinear optical medium this corresponds to the process of sum frequency generation (via chi((2)) nonlinearity) modified by the chi((3)) nonlinearity. Potential applications include an ultrafast photonic AND-gate. The simplest quantum solitons or energy eigenstates (bound-state solutions) of the interacting field Hamiltonian are obtained exactly in three space dimensions. They have a point-like structure-even though the corresponding classical theory is nonsingular. We show that the solutions can be regularized with the imposition of a momentum cut-off on the nonlinear couplings. The case of three-dimensional matter-wave solitons in coupled atomic/molecular Bose-Einstein condensates is discussed.

The carbohydrate-binding specificity and molecular modelling of Canavalia maritima and Dioclea grandiflora lectins

The carbohydrate-binding specificity of lectins from the seeds of Canavalia maritima and Dioclea grandiflora was studied by hapten-inhibition of haemagglutination using various sugars and sugar derivatives as inhibitors, including N-acetylneuraminic acid and N-acetylmuramic acid. Despite some discrepancies, both lectins exhibited a very similar carbohydrate-binding specificity as previously reported for other lectins from Diocleinae (tribe Phaseoleae, sub-tribe Diocleinae). Accordingly, both lectins exhibited almost identical hydropathic profiles and their three-dimensional models built up from the atomic coordinates of ConA looked very similar. However, docking experiments of glucose and mannose in their monosaccharide-binding sites, by comparison with the ConA-mannose complex used as a model, revealed conformational changes in side chains of the amino acid residues involved in the binding of monosaccharides. These results fully agree with crystallographic data showing that binding of specific ligands to ConA requires conformational chances of its monosaccharide-binding site.

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

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

The mapping of the local contributions of Fermi and Coulomb correlation into intracule and extracule density distributions

The contributions of the correlated and uncorrelated components of the electron-pair density to atomic and molecular intracule I(r) and extracule E(R) densities and its Laplacian functions ∇2I(r) and ∇2E(R) are analyzed at the Hartree-Fock (HF) and configuration interaction (CI) levels of theory. The topologies of the uncorrelated components of these functions can be rationalized in terms of the corresponding one-electron densities. In contrast, by analyzing the correlated components of I(r) and E(R), namely, IC(r) and EC(R), the effect of electron Fermi and Coulomb correlation can be assessed at the HF and CI levels of theory. Moreover, the contribution of Coulomb correlation can be isolated by means of difference maps between IC(r) and EC(R) distributions calculated at the two levels of theory. As application examples, the He, Ne, and Ar atomic series, the C2-2, N2, O2+2 molecular series, and the C2H4 molecule have been investigated. For these atoms and molecules, it is found that Fermi correlation accounts for the main characteristics of IC(r) and EC(R), with Coulomb correlation increasing slightly the locality of these functions at the CI level of theory. Furthermore, IC(r), EC(R), and the associated Laplacian functions, reveal the short-ranged nature and high isotropy of Fermi and Coulomb correlation in atoms and molecules