Application of Reptation Quantum Monte Carlo and Related Methods to LiH

This work investigates mathematical details and computational aspects of Metropolis-Hastings reptation quantum Monte Carlo and its variants, in addition to the Bounce method and its variants. The issues that concern us include the sensitivity of these algorithms' target densities to the position of the trial electron density along the reptile, time-reversal symmetry of the propagators, and the length of the reptile. We calculate the ground-state energy and one-electron properties of LiH at its equilibrium geometry for all these algorithms. The importance sampling is performed with a single-determinant large Slater-type orbitals (STO) basis set. The computer codes were written to exploit the efficiencies engineered into modern, high-performance computing software. Using the Bounce method in the calculation of non-energy-related properties, those represented by operators that do not commute with the Hamiltonian, is a novel work. We found that the unmodified Bounce gives good ground state energy and very good one-electron properties. We attribute this to its favourable time-reversal symmetry in its target density's Green's functions. Breaking this symmetry gives poorer results. Use of a short reptile in the Bounce method does not alter the quality of the results. This suggests that in future applications one can use a shorter reptile to cut down the computational time dramatically.

Optical properties of organic superconductor K-(BETS)2FeBr4 /

K-(BETS)2FeBr4 is a quasi-2D charge transfer organic metal with interesting electronic and magnetic properties. It undergoes a transition to an antiferromagnetic (AF) state at ambient pressure at the Neel temperature (T^^) = 2.5 K, as well as to a superconducting (SC) state at 1.1 K [1]. The temperature dependence of the electrical resistivity shows a small decrease at T;v indicating the resistivity drops as a result of the onset of the ordering of Fe'*''" spins. A sharp drop in the resistivity at 1.1 K is due to its superconducting transition. The temperature dependence of the susceptibility indicates an antiferromagnetic spin structure with the easy axis parallel to the a-axis. The specific heat at zero-field shows a large peak at about 2.4 K, which corresponds to the antiferromagnetic transition temperature (Tat) and no anomaly is observed around the superconducting transition temperature (1.1 K) demonstrating that the magnetically ordered state is not destroyed by the appearance of another phase transition (the superconducting transition) in the 7r-electron layers [1], [2]. This work presents an investigation of how the low frequency electromagnetic response is affected by the antiferromagnetic and superconducting states, as well as the onset of strong correlation. The location of the easy axis of three samples was determined and polarized thermal reflectance measurements of these «-(BETS)2FeBr4 samples oriented with their vertical axis along the a- and c axes were then carried out using a *He refrigerator cryostat and a Martin-Puplett type polarizing interferometer at various temperatures (T = 0.5 K, 1.4 K. 1.9 K, 2.8 K) above and below the superconducting state and/or antiferromagnetic state. Comparison of the SC state to the normal state along the o- and c-axes indicates a rising thermal reflectance at low frequencies (below 10 cm"' ) which may be a manifestation of the superconducting energy gap. A dip-Hke feature is detected at low frequencies (below 15 cm"') in the thermal reflectance plots which probe the antiferromagnetic state along the two axes, and may be due to the opening of a gap in the excitation spectrum as a result of the antiferromagnetism. In another set of experiments, thermal reflectance measurements carried out along the a- and c-axes at higher temperatures (10 K-80 K) show that the reflectivity decreases with increasing temperature to 60 K (the coherence temperature) above which it increases again. Comparison of the thermal reflectance plots along the a- and c-axes at higher temperatures reveals an anisotropy between these two axes. The Hagen-Rubens thermal reflectance plots corresponding to an average over the ac-plane were calculated using experimental hterature resistivity values. Comparison of the Hagen-Rubens plots with the experimental thermal reflectance along the a- and c-axes indicates that both exhibit the general trend of a decrease in thermal reflectance with increasing frequency, however the calculated Hagen-Rubens thermal reflectance at different temperatures is much lower than the experimental curves.

Thermodynamic properties and Debye-Waller factor of fee materials

We have calculated the equation of state and the various thermodynamic properties of monatomic fcc crystals by minimizing the Helmholtz free energy derived in the high temperature limit for the quasiharmonic theory, QH, and the lowest-order (cubic and quartic), 'A2, anharmonic terms of the perturbation theory, PT. The total energy in each case is obtained by adding the static energy. The calculation of the thermal properties was carried out for a nearest-neighbour central-force model of the fcc lattice by means of the appropriate thermodynamic relations. We have calculated the lattice constant, the thermal expansion, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the isothermal and adiabatic bulk moduli, and the Griineisen parameter, for the rare-gas solids Kr and Xe, and gold. Morse potential and modified Morse potential were each used to represent the atomic interaction for the three fcc materials. For most of the calculated thermodynamic properties from the QH theory, the results for Kr and Xe with the modified Morse potential show an improvement over the results for the Morse potential when compared with the experimental data. However, the results of the 'A 2 equation of state with the modified Morse potential are in good agreement with experiment only in the case of the specific heat at constant volume and at constant pressure. For Au we have calculated the lattice contribution from the QH and 'A 2 PT and the electronic contribution to the thermal properties. The electronic contribution was taken into account by using the free electron model. The results of the thermodynamic properties calculated with the modified Morse potential were similar to those obtained with the Morse potential. U sing the minimized equation of state we also calculated the Mossbauer recoilless fraction for Kr and Xe and the Debye-Waller factor (DWF) for Pb, AI, eu, Ag, and Au. The Mossbauer recoilless fraction was obtained for the above two potentials and Lennard-Jones potential. The L-J potential gives the best agreement with experiment for Kr. No experimental data exists for Xe. At low temperature the calculated DWF results for Pb, AI, and eu show a good agreement with experimental values, but at high temperature the experimental DWF results increase very rapidly. For Ag the computed values were below the expected results at all temperatures. The DWF results of the modified Morse potential for Pb, AI, eu and Ag were slightly better than those of the Morse potential. In the case of Au the calculated values were in poor agreement with experimental results. We have calculated the quasiharmonic phonon dispersion curves for Kr, Xe, eu, Ag, and Au. The calculated and experimental results of the frequencies agree quite well for all the materials except for Au where the longitudinal modes show serious discrepancies with the experimental results. In addition, the two lowest-order anharmonic contributions to the phonon frequency were derived using the Green's function method. The A 2 phonon dispersion curves have been calculated only for eu, and the results were similar to those of the QH dispersion curves. Finally, an expression for the Griineisen parameter "( has been derived from the anharmonic frequencies, and calculated for these materials. The "( results are comparable with those obtained from the thermodynamic definition.

Triarylamminium radical-cation complexes: design and magnetic properties and Base-catalyzed hydrosilylation: mechanism and substrate scope

1. Triarylamminium radical-cation complexes. The detailed study of manganese, copper and nickel metal-radical complexes with triarylamminium ligands was conducted. Stable, neutral and pseudo-octahedral coordination monometallic complexes with simple monodentate 2,2`-bipyridine ligand containing a redox-active N,N`-(4,4`-dimethoxydiphenyl-amino) substituent were synthesized and fully characterized. The one-electron oxidation process and formation of persistent radical-cation complexes was observed by cyclic voltammetry and spectroelectrochemical measurements. Evans method measurements were performed with radical-cation complexes generated by chemical one-electron oxidation with NOPF6 in acetonitrile. The experimental results indicate ferromagnetic coupling between metal and triarylamminium cation in manganese (II) complex and antiferromagnetic coupling in nickel (II) complex. This data is supported by DFT calculations which also lend weight to the  spin polarization mechanism as an operative model for magnetic exchange coupling. Neutral bimetallic complexes with a new ditopic ligand were synthesized and fully characterized, including magnetic and electrochemical studies. Chemical oxidation of these precursor complexes did not generate radical-cations, but dicationic complexes, which was confirmed by UV-vis and EPR-experiments, as well as varied temperature magnetic measurements. DFT calculations for radical-cation complexes are included. A synthetic pathway for polytopic ligand with multiple redox-active triarylamine sites was developed. The structure of the ligand is presumably suitable for -spin polarization exchange model and allows for production of polymetallic complexes having high spin ground states. 2. Base-catalyzed hydrosilylation. A simple reductive base-catalyzed hydrosilation of aldehydes and ketones was adapted to the use of the cheap, safe, and non-toxic polymethylhydrosiloxane (PMHS) instead of the common PhSiH3 and (EtO)3SiH, which present significant cost and safety concerns, respectively. The conversion of silane into pentacoordinate silicate species upon addition of a base was studied in details for the cases of phenyl silane and PMHS and is believed to be essential for the hydrosilylation process. We discovered that nucleophiles (a base or fluoride-anion) induced the rearrangement of PMHS and TMDS into light silanes: MeSiH3 and Me2SiH2, respectively. The reductive properties of PMHS under basic conditions can be attributed to the formation of methyl silane and its conversion into a silicate species. A procedure for the generation of methyl silane and its use in further efficient reductions of aldehydes and ketones has been developed. The protocol was extended to the selective reduction of esters and tertiary amides into alcohols and aldimines into amines with good isolated yields and reduction of heterocyclic compounds was attempted.

Phonon spectra and temperature variation of thermodynamic properties of fcc metals via Finnis-Sinclair type many body potentials: Sutton-Chen and improved Sutton-Chen models

Volume(density)-independent pair-potentials cannot describe metallic cohesion adequately as the presence of the free electron gas renders the total energy strongly dependent on the electron density. The embedded atom method (EAM) addresses this issue by replacing part of the total energy with an explicitly density-dependent term called the embedding function. Finnis and Sinclair proposed a model where the embedding function is taken to be proportional to the square root of the electron density. Models of this type are known as Finnis-Sinclair many body potentials. In this work we study a particular parametrization of the Finnis-Sinclair type potential, called the "Sutton-Chen" model, and a later version, called the "Quantum Sutton-Chen" model, to study the phonon spectra and the temperature variation thermodynamic properties of fcc metals. Both models give poor results for thermal expansion, which can be traced to rapid softening of transverse phonon frequencies with increasing lattice parameter. We identify the power law decay of the electron density with distance assumed by the model as the main cause of this behaviour and show that an exponentially decaying form of charge density improves the results significantly. Results for Sutton-Chen and our improved version of Sutton-Chen models are compared for four fcc metals: Cu, Ag, Au and Pt. The calculated properties are the phonon spectra, thermal expansion coefficient, isobaric heat capacity, adiabatic and isothermal bulk moduli, atomic root-mean-square displacement and Gr\"{u}neisen parameter. For the sake of comparison we have also considered two other models where the distance-dependence of the charge density is an exponential multiplied by polynomials. None of these models exhibits the instability against thermal expansion (premature melting) as shown by the Sutton-Chen model. We also present results obtained via pure pair potential models, in order to identify advantages and disadvantages of methods used to obtain the parameters of these potentials.