Magnetically induced currents and magnetic response properties of molecules

The magnetically induced currents in organic monoring and multiring molecules, in Möbius shaped molecules and in inorganic all-metal molecules have been investigated by means of the Gauge-including magnetically induced currents (GIMIC) method. With the GIMIC method, the ring-current strengths and the ring-current density distributions can be calculated. For open-shell molecules, also the spin current can be obtained. The ring-current pathways and ring-current strengths can be used to understand the magnetic resonance properties of the molecules, to indirectly identify the effect of non-bonded interactions on NMR chemical shifts, to design new molecules with tailored properties and to discuss molecular aromaticity. In the thesis, the magnetic criterion for aromaticity has been adopted. According to this, a molecule which has a net diatropic ring current might be aromatic. Similarly, a molecule which has a net paratropic current might be antiaromatic. If the net current is zero, the molecule is nonaromatic. The electronic structure of the investigated molecules has been resolved by quantum chemical methods. The magnetically induced currents have been calculated with the GIMIC method at the density-functional theory (DFT) level, as well as at the self-consistent field Hartree-Fock (SCF-HF), at the Møller-Plesset perturbation theory of the second order (MP2) and at the coupled-cluster singles and doubles (CCSD) levels of theory. For closed-shell molecules, accurate ring-current strengths can be obtained with a reasonable computational cost at the DFT level and with rather small basis sets. For open-shell molecules, it is shown that correlated methods such as MP2 and CCSD might be needed to obtain reliable charge and spin currents. The basis set convergence has to be checked for open-shell molecules by performing calculations with large enough basis sets. The results discussed in the thesis have been published in eight papers. In addition, some previously unpublished results on the ring currents in the endohedral fullerene Sc3C2@C80 and in coronene are presented. It is shown that dynamical effects should be taken into account when modelling magnetic resonance parameters of endohedral metallofullerenes such as Sc3C2@C80. The ring-current strengths in a series of nano-sized hydrocarbon rings are related to static polarizabilities and to H-1 nuclear magnetic resonance (NMR) shieldings. In a case study on the possible aromaticity of a Möbius-shaped [16]annulene we found that, according to the magnetic criterion, the molecule is nonaromatic. The applicability of the GIMIC method to assign the aromatic character of molecules was confirmed in a study on the ring currents in simple monocylic aromatic, homoaromatic, antiaromatic, and nonaromatic hydrocarbons. Case studies on nanorings, hexaphyrins and [n]cycloparaphenylenes show that explicit calculations are needed to unravel the ring-current delocalization pathways in complex multiring molecules. The open-shell implementation of GIMIC was applied in studies on the charge currents and the spin currents in single-ring and bi-ring molecules with open shells. The aromaticity predictions that are made based on the GIMIC results are compared to other aromaticity criteria such as H-1 NMR shieldings and shifts, electric polarizabilities, bond-length alternation, as well as to predictions provided by the traditional Hückel (4n+2) rule and its more recent extensions that account for Möbius twisted molecules and for molecules with open shells.

Conformational Studies of Some Alpha-substituted Acetones CH3COCH2XCH3 Where X=CH2, O Or S And Phenacyl Sulfide

Both semiempirical and ab initio calculations are reported for conformational studies of a series of alpha-substituted acetones CH3COCH2XCH3 where X = CH2, O or S and of phenacyl sulfide PhCOCH2SCH3. For conformational studies in the lowest triplet state of these molecules, the MINDO/3 method was employed in the unrestricted Hartree-Fock frame. Results reveal that rotation around a bond alpha to the carbonyl group is more favourable than that around the beta bond. The preferred conformations in the lowest triplet state are nearly the same as in the ground state.

Photochemical intramolecular gamma-hydrogen abstraction in a ketone and an alkene and the role of tunneling of hydrogen

Intramolecular gamma-hydrogen abstraction reactions were examined in pentane-2-one and 2-methyl-1-pentene in their lowest triplet states using the AM1 semi-empirical molecular orbital method with the complete geometry optimization in the unrestricted Hartree-Fock frame. The results reveal that the oxygen atom of the carbonyl group and the end carbon atom of the olefinic bond acquire high free valence and spin density indices in their respective lowest triplet states, leading to abstraction of hydrogen from the gamma-position relative to the carbonyl and olefinic bonds. The theoretical energy profiles fit with a polynomial and the probability of tunneling of hydrogen was estimated by the WKB (Wentzel, Kramer and Brillouin) method. The results, after thermal averaging of the rate constants, reveal that tunneling of hydrogen is significant at room temperature.

Further studies on Norrish type II reactions including a reaction in the first excited singlet state and cyclization of 1:4 biradicals

The Norrish type II processes of methyl-2,2-dimethyl- cyclopropyl ketone, alpha-alkoxy acetones and alkyl pyruvates have been examined using the AM1 semi-empirical molecular orbital method with complete geometry optimization at the partial configuration interaction level in the restricted Hartree-Fock (RHF) frame. The results reveal that the methyl-substituted cyclopropyl ketone has a constrained geometry favourable for hydrogen abstraction from the gamma-position relative to the carbonyl group in the excited singlet state. The presence of the ether oxygen atom in the beta-position relative to the carbonyl group in alkoxy acetones and alkyl pyruvates leads to increased reactivity relative to alkyl monoketones and diketones respectively. The cyclization of 1:4 biradicals has been studied in the unrestricted Hartree-Fock (UHF) frame, and the results reveal that the 1:4 biradical derived from alkoxy acetones readily cyclizes to form oxetanols. On the other hand, in the 1:4 biradicals derived from methyl-substituted cyclopropyl ketone, the three-membered ring breaks readily to form an enol intermediate. Delocalization of an odd electron in 1:4 biradicals derived from alkyl pyruvates is thought to make cyclization difficult.

Asymptotic behavior and interpretation of virtual states: The effects of confinement and of basis sets

A real-space high order finite difference method is used to analyze the effect of spherical domain size on the Hartree-Fock (and density functional theory) virtual eigenstates. We show the domain size dependence of both positive and negative virtual eigenvalues of the Hartree-Fock equations for small molecules. We demonstrate that positive states behave like a particle in spherical well and show how they approach zero. For the negative eigenstates, we show that large domains are needed to get the correct eigenvalues. We compare our results to those of Gaussian basis sets and draw some conclusions for real-space, basis-sets, and plane-waves calculations. (C) 2016 AIP Publishing LLC.

Fermionic quantum systems. Part I: Phase transitions in quantum dots. Part II: Nuclear matter on a lattice

In the first part I perform Hartree-Fock calculations to show that quantum dots (i.e., two-dimensional systems of up to twenty interacting electrons in an external parabolic potential) undergo a gradual transition to a spin-polarized Wigner crystal with increasing magnetic field strength. The phase diagram and ground state energies have been determined. I tried to improve the ground state of the Wigner crystal by introducing a Jastrow ansatz for the wave function and performing a variational Monte Carlo calculation. The existence of so called magic numbers was also investigated. Finally, I also calculated the heat capacity associated with the rotational degree of freedom of deformed many-body states and suggest an experimental method to detect Wigner crystals.

The second part of the thesis investigates infinite nuclear matter on a cubic lattice. The exact thermal formalism describes nucleons with a Hamiltonian that accommodates on-site and next-neighbor parts of the central, spin-exchange and isospin-exchange interaction. Using auxiliary field Monte Carlo methods, I show that energy and basic saturation properties of nuclear matter can be reproduced. A first order phase transition from an uncorrelated Fermi gas to a clustered system is observed by computing mechanical and thermodynamical quantities such as compressibility, heat capacity, entropy and grand potential. The structure of the clusters is investigated with the help two-body correlations. I compare symmetry energy and first sound velocities with literature and find reasonable agreement. I also calculate the energy of pure neutron matter and search for a similar phase transition, but the survey is restricted by the infamous Monte Carlo sign problem. Also, a regularization scheme to extract potential parameters from scattering lengths and effective ranges is investigated.

Linearization Method and Linear Complexity

We focus on the relationship between the linearization method and linear complexity and show that the linearization method is another effective technique for calculating linear complexity. We analyze its effectiveness by comparing with the logic circuit method. We compare the relevant conditions and necessary computational cost with those of the Berlekamp-Massey algorithm and the Games-Chan algorithm. The significant property of a linearization method is that it needs no output sequence from a pseudo-random number generator (PRNG) because it calculates linear complexity using the algebraic expression of its algorithm. When a PRNG has n [bit] stages (registers or internal states), the necessary computational cost is smaller than O(2n). On the other hand, the Berlekamp-Massey algorithm needs O(N2) where N ( 2n) denotes period. Since existing methods calculate using the output sequence, an initial value of PRNG influences a resultant value of linear complexity. Therefore, a linear complexity is generally given as an estimate value. On the other hand, a linearization method calculates from an algorithm of PRNG, it can determine the lower bound of linear complexity.

Shell-model Hamiltonian from self-consistent mean-field model: N = Z nuclei

We propose a procedure to determine the effective nuclear shell-model Hamiltonian in a truncated space from a self-consistent mean-field model, e.g., the Skyrme model. The parameters of pairing plus quadrupole-quadrupole interaction with monopole force are obtained so that the potential energy surface of the Skyrme Hartree-Fock + BCS calculation is reproduced. We test our method for N = Z nuclei in the fpg- and sd-shell regions. It is shown that the calculated energy spectra with these parameters are in a good agreement with experimental data, in which the importance of the monopole interaction is discussed. This method may represent a practical way of defining the Hamiltonian for general shell-model calculations. (C) 2009 Elsevier B.V. All rights reserved.

Theoretical study of p-d hybridization in Co perovskite

The real-space recursion method and unrestricted Hartree-Fock approximation have been applied to calculate the density of states of various Co perovskite, CeCoO3, SrCoO3 and Sr1-xCexCoO3. We have studied the magnetically ordered states of these Co perovskites in an enlarged double cell, and find its various magnetic structures due to the occupancy of 3d band and its interaction with neighboring Co ions. In this study, we have studied the p-d hybridization of the three Co perovskites, we find t(2g) electrons are localized and the flat e(g) band is responsible for the itinerant behavior, and although the rare earth elements itself contribute little to the DOS at the Fermi energy, the DOS at Fermi energy and the magnetic moment changed consequently because of different valence of Co ions in these compounds and p-d hybridization effect is very important. (C) 2009 Elsevier B.V. All rights reserved.

Improved BCS-type pairing for the relativistic mean-field theory

A density-dependent delta interaction (DDDI) is proposed in the formalism of BCS-type pairing correlations for exotic nuclei whose Fermi surfaces are close to the threshold of the unbound state. It provides the possibility to pick up those states whose wave functions are concentrated in the nuclear region by making the pairing matrix elements state dependent. On this basis, the energy level distributions, occupations, and ground-state properties are self-consistently studied in the RMF theory with deformation. Calculations are performed for the Sr isotopic chain. A good description of the total energy per nucleon, deformations, two-neutron separation energies and isotope shift from the proton drip line to the neutron drip line is found. Especially, by comparing the single-particle structure from the DDDI pairing interaction with that from the constant pairing interaction for a very neutron-rich nucleus it is demonstrated that the DDDI pairing method improves the treatment of the pairing in the continuum.

Two-electron-one-photon M1 and E2 transitions between the states of the 2p(3) and 2s(2)2p odd configurations for B-like ions with 18 <= Z <= 92

Two-electron-one-photon (TEOP) M1 and E2 transition energies, line strengths and transition probabilities between the states of the 2p(3) and 2s(2)2p odd configurations for B-like ions with 18 <= Z <= 92 have been calculated using the GRASP2K package based on the multiconfiguration Dirac-Hartree-Fock (MCDHF) method. Employing active-space techniques to expand the configuration list, we have systematically considered the valence, core-valence and core-core electron correlation effects. Breit interaction and quantum electrodynamical (QED) effects were also included to correct atomic state wavefunctions and the corresponding energies. Influences of electron correlation, Breit interaction and QED effects on transition energies and line strengths of the TEOP M1 and E2 transitions were analysed in detail. The present results were also compared with other theoretical and experimental values.

Gas-phase ion-molecule reactions of C-60 with the plasma of methyl acrylate

The gas-phase ion-molecule reactions of C-60 with the plasma generated from methyl acrylate under self-chemical ionization conditions were studied by use of a triple-quadrupole mass spectrometer. The adduct cation [C60C3H3O](+) and protonated molecular ion [C60H](+) were observed as the major product ions. The former adduct ion is formed by electrophilic reaction of C-60 with the ion [CH2=CHCO](+), a main fragment ion resulting from the methyl acrylate molecular ion [CH2=CHCOOCH3](+) through alpha cleavage. The latter ion is generated by proton transfer from protonated methyl acrylate to C-60. Semi-empirical quantum chemical calculations have been performed for the eight possible isomers of [C60C3H3O](+) at the Hartree-Fock level by use of the AMI method. The results show three types of cycloadducts as the most stable structures among the possible isomers.

Structures and vibrational spectra of C-2 and LaC2+ clusters

C-2 and LaC2+ were studied using Hartree-Fock(HF), B3LYP (Becke 3-paremeter-Lee-Yang-Parr) density functional method, second-order Moller-Plesset perturbation (MP2) and coupled cluster singles and doubles with non-iterative triples(CCSD(T)) methods. The basis set employed was LANL1DZ. Geometries, vibrational frequencies and other quantities were reported. The results showed that for C-2, all the methods performed well for low spin state (singlet), while only HF and B3LYP remained so for high spin state (triplet). For LaC2+, four isomers were presented and fully optimized. The results suggested that linear isomers with C-infinity v and D-infinity h symmetries were predicted to be saddle points on the energy surface for all the methods, while for isomers with C-2 upsilon and C-s symmetries, they were local minima except C-2 upsilon at B3LYP level, and were isoenergetic at HF, MP2 and CCSD(T) levels, near isoenergetic at B3LYP level. From the differences between HOMO and LUMO, it is also known that the isomers with C-2 upsilon and C-s symmetries offer the largest values and therefore correspond to the most stable structure. For La-C bond lengths, B3LYP gives the shortest, the order is B3LYPmethods.

Theoretical investigation of LaCn+ (n=2-8) clusters

LaCn+ (n = 2-8) have been studied using Hartree-Fock (HF) and B3LYP density functional method. The results indicated that at both levels, isomers with C-2v, C-s symmetry for n = 2, and edge insertion isomer for n = 4, 6, 8, as well as edge binding isomer for n = 3, 5, 7 were found as ground states. This is in good agreement with experimental results. The exceptional case is for n = 6 at B3LYP level, in which edge insertion and edge binding isomers were computed to be near isoenergetic. (C) 1997 Elsevier Science B.V.

Latest development of the combinatorial model of nuclear level densities

The combinatorial model of nuclear level densities has now reached a level of accuracy comparable to that of the best global analytical expressions without suffering from the limits imposed by the statistical hypothesis on which the latter expressions rely. In particular, it provides, naturally, non-Gaussian spin distribution as well as non-equipartition of parities which are known to have an impact on cross section predictions at low energies [1, 2, 3]. Our previous global models developed in Refs. [1, 2] suffered from deficiencies, in particular in the way the collective effects - both vibrational and rotational - were treated. We have recently improved this treatment using simultaneously the single-particle levels and collective properties predicted by a newly derived Gogny interaction [4], therefore enabling a microscopic description of energy-dependent shell, pairing and deformation effects. In addition for deformed nuclei, the transition to sphericity is coherently taken into account on the basis of a temperature-dependent Hartree-Fock calculation which provides at each temperature the structure properties needed to build the level densities. This new method is described and shown to give promising results with respect to available experimental data.