Uncertainty budget for determinations of mean isomer shift from Mössbauer spectra

The magnetite/maghemite content within iron oxide nanoparticles can be determined using the mean isomer shift (\(\overline {\delta }\)). However, accurate characterisation of the composition is limited by the uncertainty associated with \(\overline {\delta }\). We have identified four independent sources of uncertainty and developed a quantitative expression for the uncertainty budget. Sources of uncertainty are categorised as follows: that from the fitting of the Mössbauer spectrum (σ _fit), that of the calibration of the α-Fe reference spectrum (σ _cal), thermal corrections to the spectrum due to second order Doppler shift (SODS) (σ _Δδ ) and other experimental errors (σ _err). Each contribution is discussed in detail using ⁵⁷Fe Mössbauer spectra obtained from an iron oxide nanoparticle system at temperatures between 16 K and 295 K on different spectrometers in two different laboratories.     

Linear-scaling symmetric square-root decomposition of the overlap matrix

We present a robust linear-scaling algorithm to compute the symmetric square-root or Lowdin decomposition of the atomic-orbital overlap matrix. The method is based on Newton-Schulz iterations with a new approach to starting matrices. Calculations on 12 chemically and structurally diverse molecules demonstrate the efficiency and reliability of the method. Furthermore, the calculations show that linear scaling is achieved.     

Generalising localisation schemes of orthogonal orbitals to the localisation of non-orthogonal orbitals

ABSTRACT

When a set of orthogonal orbitals is localised using orthogonal transformations, the orbital metric is conserved as the unit matrix during the localisation. In this paper we describe how non-orthogonal orbitals may be localised by requiring that the orbital metric is conserved. In particular, we demonstrate how exponential mappings may be used to parametrise orbital transformations such that the orbital metric is conserved. Using this parametrisation, the localisation of non-orthogonal orbitals becomes a generalisation of the localisation for orthogonal orbitals, where the conservation of a unit metric is replaced by the conservation of a non-unit metric. As a result, standard orbital localisation functions and optimisation algorithms that have been used for localising orthogonal orbitals may also be used to localise non-orthogonal orbitals. Numerical illustrations show that the thickness of the orbital tails may be reduced when a set of non-orthogonal orbitals is localised compared to the orthogonal counterpart, and further that the physical interpretation of occupied orbitals for representing chemical bonds is preserved.     

Shear and dielectric responses of propylene carbonate, tripropylene glycol, and a mixture of two secondary amides

Propylene carbonate and a mixture of two secondary amides, N-methylformamide and N-ethylacetamide, are investigated by means of broadband dielectric and mechanical shear spectroscopy. The similarities between the rheological and the dielectric responses of these liquids and of the previously investigated tripropylene glycol are discussed within a simple approach that employs an electrical circuit for describing the frequency-dependent behavior of viscous materials. The circuit is equivalent to the Gemant-DiMarzio-Bishop model, but allows for a negative capacitive element. The circuit can be used to calculate the dielectric from the mechanical response and vice versa. Using a single parameter for a given system, good agreement between model calculations and experimental data is achieved for the entire relaxation spectra, including secondary relaxations and the Debye-like dielectric peak in the secondary amides. In addition, the predictions of the shoving model are confirmed for the investigated liquids.     

Communication: Identical temperature dependence of the time scales of several linear-response functions of two glass-forming liquids

The frequency-dependent dielectric constant, shear and adiabatic bulk moduli, longitudinal thermal expansion coefficient, and longitudinal specific heat have been measured for two van der Waals glass-forming liquids, tetramethyl-tetraphenyl-trisiloxane (DC704) and 5-polyphenyl-4-ether. Within the experimental uncertainties the loss-peak frequencies of the measured response functions have identical temperature dependence over a range of temperatures, for which the Maxwell relaxation time varies more than nine orders of magnitude. The time scales are ordered from fastest to slowest as follows: Shear modulus, adiabatic bulk modulus, dielectric constant, longitudinal thermal expansion coefficient, and longitudinal specific heat. The ordering is discussed in light of the recent conjecture that van der Waals liquids are strongly correlating, i.e., approximate single-parameter liquids.     

Isomorph Invariance of Higher-Order Structural Measures in Four Lennard–Jones Systems

In the condensed liquid phase, both single- and multicomponent Lennard–Jones (LJ) systems obey the “hidden-scale-invariance” symmetry to a good approximation. Defining an isomorph as a line of constant excess entropy in the thermodynamic phase diagram, the consequent approximate isomorph invariance of structure and dynamics in appropriate units is well documented. However, although all measures of the structure are predicted to be isomorph invariant, with few exceptions only the radial distribution function (RDF) has been investigated. This paper studies the variation along isomorphs of the nearest-neighbor geometry quantified by the occurrence of Voronoi structures, Frank–Kasper bonds, icosahedral local order, and bond-orientational order. Data are presented for the standard LJ system and for three binary LJ mixtures (Kob–Andersen, Wahnström, NiY2). We find that, while the nearest-neighbor geometry generally varies significantly throughout the phase diagram, good invariance is observed along the isomorphs. We conclude that higher-order structural correlations are no less isomorph invariant than is the RDF.     

Shape resonances and partial photoemission cross sections of solid SF6 and CCl4

Photoelectron energy distribution curves from solid films of SF₆ and CCl₄ have been measured in the photon energy range 10 ? hr ? 40 eV using synchrotron radiation. The binding energies, peak widths and relative partial cross sections have been determined. In the photoelectron spectra a 1:1 correspondence to the gas phase is observed for the occupied molecular orbitals, and a straightforward assignment of the occupied valence bands emerges. Furthermore, the cross sections of the individual orbitals show for both samples great similarities to the gas phase. For SF₆ detailed structures are visible in the cross sections which are only partly interpreted as shape resonances. A new assignment for the 6t_1u shape resonance is proposed and the resonance energies are related to X-ray absorption and electron scattering data. Furthermore a comparison of the total photoemission cross section to the optical reflection spectrum of solid SF₆ is presented. For CCl₄ less structures are ob served in the partial cross sections. They are all interpreted as shape resonances. An energetic scheme of the virtual orbitals is proposed for CCl₄.     

A determinantal approach to spin-orbit configuration interaction

A general configuration interaction (CI) algorithm incorporating one- and two-electron spin-orbit operators is presented. The algorithm is determinant based and enables the use of highly vectorized non-relativistic algorithms in the most operation-intensive part. Excitations between α and β spin orbitals are avoided in the time consuming parts by performing separate S ₊ and S ₋ operations. The relativistic CI expansions are often very large, so the algorithms require only the presence of segments of vectors in memory. Double-group symmetry is fully accounted for and time-reversal symmetry is exploited for both even and odd numbers of electrons. Received: 12 October 1996 / Accepted: 24 January 1997     

A relativistic 4-component general-order multi-reference coupled cluster method: initial implementation and application to HBr

We present the initial implementation of a determinant-based general-order coupled cluster method which fully accounts for relativistic effects within the four-component framework. The method opens the way for the treatment of multi-reference problems through a state-selective expansion of the model space. The evaluation of the coupled cluster vector function is carried out via relativistic configuration interaction expansions. The implementation is based on a large-scale configuration interaction technique, which may efficiently treat long determinant expansions of more than 10⁸ terms. We demonstrate the capabilities of the new method in calculations of complete potential energy curves of the HBr molecule. The inclusion of spin–orbit interaction and higher excitations than coupled cluster double excitations, either by multi-reference model spaces or the inclusion of full iterative triple excitations, lead to highly accurate results for spectral constants of HBr. An erratum to this article can be found at     

Linear-scaling implementation of molecular electronic self-consistent field theory

A linear-scaling implementation of Hartree-Fock and Kohn-Sham self-consistent field (SCF) theories is presented and illustrated with applications to molecules consisting of more than 1000 atoms. The diagonalization bottleneck of traditional SCF methods is avoided by carrying out a minimization of the Roothaan-Hall (RH) energy function and solving the Newton equations using the preconditioned conjugate-gradient (PCG) method. For rapid PCG convergence, the Lowdin orthogonal atomic orbital basis is used. The resulting linear-scaling trust-region Roothaan-Hall (LS-TRRH) method works by the introduction of a level-shift parameter in the RH Newton equations. A great advantage of the LS-TRRH method is that the optimal level shift can be determined at no extra cost, ensuring fast and robust convergence of both the SCF iterations and the level-shifted Newton equations. For density averaging, the authors use the trust-region density-subspace minimization (TRDSM) method, which, unlike the traditional direct inversion in the iterative subspace (DIIS) scheme, is firmly based on the principle of energy minimization. When combined with a linear-scaling evaluation of the Fock/Kohn-Sham matrix (including a boxed fitting of the electron density), LS-TRRH and TRDSM methods constitute the linear-scaling trust-region SCF (LS-TRSCF) method. The LS-TRSCF method compares favorably with the traditional SCF/DIIS scheme, converging smoothly and reliably in cases where the latter method fails. In one case where the LS-TRSCF method converges smoothly to a minimum, the SCF/DIIS method converges to a saddle point.