Microscopic and macroscopic polarization within a combined quantum mechanics and molecular mechanics model

A polarizable quantum mechanics and molecular mechanics model has been extended to account for the difference between the macroscopic electric field and the actual electric field felt by the solute molecule. This enables the calculation of effective microscopic properties which can be related to macroscopic susceptibilities directly comparable with experimental results. By seperating the discrete local field into two distinct contribution we define two different microscopic properties, the so-called solute and effective properties. The solute properties account for the pure solvent effects, i.e., effects even when the macroscopic electric field is zero, and the effective properties account for both the pure solvent effects and the effect from the induced dipoles in the solvent due to the macroscopic electric field. We present results for the linear and nonlinear polarizabilities of water and acetonitrile both in the gas phase and in the liquid phase. For all the properties we find that the pure solvent effect increases the properties whereas the induced electric field decreases the properties. Furthermore, we present results for the refractive index, third-harmonic generation (THG), and electric field induced second-harmonic generation (EFISH) for liquid water and acetonitrile. We find in general good agreement between the calculated and experimental results for the refractive index and the THG susceptibility. For the EFISH susceptibility, however, the difference between experiment and theory is larger since the orientational effect arising from the static electric field is not accurately described

Chemodiversity and molecular plasticity: recognition processes as explored by property spaces.

In the last few years, a need to account for molecular flexibility in drug-design methodologies has emerged, even if the dynamic behavior of molecular properties is seldom made explicit. For a flexible molecule, it is indeed possible to compute different values for a given conformation-dependent property and the ensemble of such values defines a property space that can be used to describe its molecular variability; a most representative case is the lipophilicity space. In this review, a number of applications of lipophilicity space and other property spaces are presented, showing that this concept can be fruitfully exploited: to investigate the constraints exerted by media of different levels of structural organization, to examine processes of molecular recognition and binding at an atomic level, to derive informative descriptors to be included in quantitative structure--activity relationships and to analyze protein simulations extracting the relevant information. Much molecular information is neglected in the descriptors used by medicinal chemists, while the concept of property space can fill this gap by accounting for the often-disregarded dynamic behavior of both small ligands and biomacromolecules. Property space also introduces some innovative concepts such as molecular sensitivity and plasticity, which appear best suited to explore the ability of a molecule to adapt itself to the environment variously modulating its property and conformational profiles. Globally, such concepts can enhance our understanding of biological phenomena providing fruitful descriptors in drug-design and pharmaceutical sciences.

On the effect of molecular structure and composition on the photofragmentation of biomolecules

The research on the interaction between radiation and biomolecules pro-vides valuable information for both radiobiology and molecular physics. While radiobiology is interested in the damage inflicted on the molecule upon irradiation, molecular physics exploits these studies to obtain infor-mation about the physical properties of the molecule and the quantum me-chanical processes involved in the interaction. This thesis work investigated how a small change in the structure or composition of a biomolecule changes the response of the molecule to ioniz-ing radiation. Altogether eight different biomolecules were studied: nucleo-sides uridine, 5-methyluridine and thymidine; amino acids alanine, cysteine and serine; and halogenated acetic acids chloro- and bromoacetic acids. The effect of ionizing radiation on these molecules was studied on molecular level, investigating the samples in gas phase. Synchrotron radiation of VUV or soft x-ray range was used to ionize sample molecules, and the subsequent fragmentation processes were investigated with ion mass spectroscopy and ion-ion-electron coincidence spectroscopy. The comparison between the three nucleosides revealed that adding or removing a single functional group can affect not only the bonds from which the molecule ruptures upon ionization but also the charge localiza-tion in the formed fragments. Studies on amino acids and halogenated acetic acids indicated that one simple substitution in the molecule can dramatical-ly change the extent of fragmentation. This thesis work also demonstrates that in order to steer the radiation-induced fragmentation of the molecules, it is not always necessary to alter the amount of energy deposited on the molecules but selecting a suitable substitution may suffice.

Atomic and ionic radii of superheavy elements

Atomic and ionic radii are presented for the elements E104-E120 and E156-E172. It is shown that a number of effects correlated with the large relativistic contraction of orbitals with low angular momentum leads to smaller atoms for higher atomic numbers. It is expected that Cs is the largest atom in nature.

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

Selectivity of bis(calix[4]diquinone) ionophores towards metal ions in solvent dimethylsulfoxide: a molecular mechanics and molecular dynamics study

Molecular modelling studies have been carried out on two bis(calix[4]diqu(inone) ionophores, each created from two (calix[4]diquinone)arenes bridged at their bottom rims via alkyl chains (CH2)(n), 1: n = 3, 2; n = 4, in order to understand the reported selectivity of these ligands towards different sized metal ions such as Na+, K+, Rb+, and Cs+ in dmso solution. Conformational. analyses have been carried out which show that in the lowest energy conformations of the two macrocycles, the individual calix[4]diquinones exhibit a combination of partial cone, 1,3-alternate and cone conformations. The interactions of these alkali metals with the macrocycles have been studied in the gas phase and in a periodic box of solvent dmso by molecular mechanics and molecular dynamics calculations. Molecular mechanics calculations have been carried out on the mode of entry of the ions into the macrocycles and suggest that this is likely to occur from the side of the central cavity, rather than through the main axis of the calix[4]diquinones. There are energy barriers of ca. 19 kcal mol(-1) for this entry path in the gas phase, but in solution no energy barrier is found. Molecular dynamics simulations show that in both 1 and 2, though particularly in the latter macrocycle, one or two solvent molecules are bonded to the metal throughout the course of the simulation, often to the exclusion, of one or more of the ether oxygen atoms. By contrast the carbonyl oxygen atoms remain bonded to the metal atoms throughout with bond lengths that remain significantly less than those to the ether oxygen atoms. Free energy perturbation studies have been carried out in dmso and indicate that for 1, the selectivity follows the order Rb+ approximate to K+ > Cs+ >> Na+, which is partially in agreement with the experimental results. The energy differences are small and indeed the ratio between stability constants found for Cs+ and K+ complexes is only 0.60, showing that 1 has only a slight preference for K+. For the larger receptor 2, which is better suited to metal complexation, the binding affinity follows the pattern Cs+ >> Rb+ >> K+ >> Na+, with energy differences of 5.75, 2.61, 2.78 kcal mol(-1) which is perfectly consistent with experimental results.

Synthesis, mechanism of formation, and molecular orbital calculations of arylamidoximes

A simple and easy synthesis of ten arylamidoximes from arylnitriles and hydroxylamine is described. The formation of the arylamides has been observed to a much lesser extent in the present work. A new mechanism for the formation of arylamidoximes, as well as arylamides, from arylnitriles and hydroxylamine is suggested. Quantum mechanical calculations have been carried out to support this mechanism. The enthalpy of formation in conjunction with atomic charges of the reactants and intermediates helped to understand more about the generation of the products.

Bifurcations and bistability in cavity-assisted photoassociation of Bose-Einstein-condensed molecules

We study the photoassociation of Bose-Einstein condensed atoms into molecules using an optical cavity field. The driven cavity field introduces a dynamical degree of freedom into the photoassociation process, whose role in determining the stationary behavior has not previously been considered. The semiclassical stationary solutions for the atom and molecules as well as the intracavity field are found and their stability and scaling properties are determined in terms of experimentally controllable parameters including driving amplitude of the cavity and the nonlinear interactions between atoms and molecules. For weak cavity driving, we find a bifurcation in the atom and molecule number occurs that signals a transition from a stable steady state to nonlinear Rabi oscillations. For a strongly driven cavity, there exists bistability in the atom and molecule number.

Quantum control of electromagnetically induced transparency dispersion via atomic tunneling in a double-well Bose-Einstein condensate

Electromagnetically induced transparency (EIT) is an important tool for controlling light propagation and nonlinear wave mixing in atomic gases with potential applications ranging from quantum computing to table top tests of general relativity. Here we consider EIT in an atomic Bose-Einstein condensate (BEC) trapped in a double-well potential. A weak probe laser propagates through one of the wells and interacts with atoms in a three-level Lambda configuration. The well through which the probe propagates is dressed by a strong control laser with Rabi frequency Omega(mu), as in standard EIT systems. Tunneling between the wells at the frequency g provides a coherent coupling between identical electronic states in the two wells, which leads to the formation of interwell dressed states. The macroscopic interwell coherence of the BEC wave function results in the formation of two ultranarrow absorption resonances for the probe field that are inside of the ordinary EIT transparency window. We show that these new resonances can be interpreted in terms of the interwell dressed states and the formation of a type of dark state involving the control laser and the interwell tunneling. To either side of these ultranarrow resonances there is normal dispersion with very large slope controlled by g. We discuss prospects for observing these ultranarrow resonances and the corresponding regions of high dispersion experimentally.

Wettability, optical properties and molecular structure of plasma polymerized diethylene glycol dimethyl ether

Modern industry has frequently employed ethylene glycol ethers as monomers in plasma polymerization process to produce different types of coatings. In this work we used a stainless steel plasma reactor to grow thin polymeric films from low pressure RF excited plasma of diethylene glycol dimethyl ether. Plasmas were generated at 5W RF power in the range of 16 Pa to 60 Pa. The molecular structure of plasma polymerized films and their optical properties were analyzed by Fourier Transform Infrared Spectroscopy (FTIR) and UltravioletVisible Spectroscopy, respectively. The IR spectra show C-H stretching at 3000-2900 cm(-1), C=O stretching at 1730-1650 cm(-1), C-H bending at 1440-1380 cm(-1), C-O and C-O-C stretching at 1200-1000 cm(-1). The refraction index was around 1.5 and the optical gap calculated from absorption coefficient presented value near 3.8 eV. Water contact angle of the films ranged from 40 degrees to 35 degrees with corresponding surface energy from 66 to 73x10(-7) J. Because of its favorable optical and hydrophilic characteristics these films can be used in ophthalmic industries as glass lenses coatings.

Quantum coherent tunneling between two atomic-molecular Bose-Einstein condensates

We study the quantum coherent tunneling dynamics of two weakly coupled atomic-molecular Bose-Einstein condensates (AMBEC). A weak link is supposed to be provided by a double-well trap. The regions of parameters where the macroscopic quantum localization of the relative atomic population occurs are revealed. The different dynamical regimes are found depending on the value of nonlinearity, namely, coupled oscillations of population imbalance of atomic and molecular condensate, including irregular oscillations regions, and macroscopic quantum self trapping regimes. Quantum means and quadrature variances are calculated for population of atomic and molecular condensates and the possibility of quadrature squeezing is shown via stochastic simulations within P-positive phase space representation method. Linear tunnel coupling between two AMBEC leads to correlations in quantum statistics.

LARGE ATOMIC and NUCLEAR THREE-BODY SYSTEMS: SCATTERING and BINDING

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Surface morphology and molecular organization of lignins in Langmuir-Blodgett films

Atomic force microscopy (AFM) and Fourier transform infrared spectroscopy (FTIR) are used to investigate molecular organization in Langmuir-Blodgett (LB) films of two kinds of lignins. The lignins were extracted from sugar cane bagasse using distinct extraction processes and are referred to here as ethanol lignin (EL) and saccharification lignin (SAC). AFM images show that LB films from EL have a flat surface in comparison with those from SAC. For the latter, ellipsoidal aggregates are seen oriented perpendicularly to the substrate. This result is confirmed by a combination of transmission and reflection FTIR measurements, which also point to lignin aggregates preferentially oriented perpendicularly to the substrate. For LB films from EL, on the other hand, aggregates are preferentially oriented parallel to the substrate, again consistent with the flat surface observed in AFM data. The vibrational spectroscopy data for cast films from both lignins show random molecular organization, as one should expect.