MULTISCALE EXAMINATION AND MODELING OF ELECTRON TRANSPORT IN NANOSCALE MATERIALS AND DEVICES

For half a century the integrated circuits (ICs) that make up the heart of electronic devices have been steadily improving by shrinking at an exponential rate. However, as the current crop of ICs get smaller and the insulating layers involved become thinner, electrons leak through due to quantum mechanical tunneling. This is one of several issues which will bring an end to this incredible streak of exponential improvement of this type of transistor device, after which future improvements will have to come from employing fundamentally different transistor architecture rather than fine tuning and miniaturizing the metal-oxide-semiconductor field effect transistors (MOSFETs) in use today. Several new transistor designs, some designed and built here at Michigan Tech, involve electrons tunneling their way through arrays of nanoparticles. We use a multi-scale approach to model these devices and study their behavior. For investigating the tunneling characteristics of the individual junctions, we use a first-principles approach to model conduction between sub-nanometer gold particles. To estimate the change in energy due to the movement of individual electrons, we use the finite element method to calculate electrostatic capacitances. The kinetic Monte Carlo method allows us to use our knowledge of these details to simulate the dynamics of an entire device— sometimes consisting of hundreds of individual particles—and watch as a device ‘turns on’ and starts conducting an electric current. Scanning tunneling microscopy (STM) and the closely related scanning tunneling spectroscopy (STS) are a family of powerful experimental techniques that allow for the probing and imaging of surfaces and molecules at atomic resolution. However, interpretation of the results often requires comparison with theoretical and computational models. We have developed a new method for calculating STM topographs and STS spectra. This method combines an established method for approximating the geometric variation of the electronic density of states, with a modern method for calculating spin-dependent tunneling currents, offering a unique balance between accuracy and accessibility.

Crystal chemistry and hydrogen bonding of rustumite Ca10(Si2O7)2(SiO4)(OH)2Cl2 with variable OH, Cl, F

Three samples of the skarn mineral rustumite Ca10(Si2O7)2(SiO4)(OH)2Cl2, space group C2/c, a ≈7.6, b ≈ 18.5, c ≈ 15.5 Å, β ≈ 104°, with variable OH, Cl, F content were investigated by electron microprobe, single-crystal X-ray structure refinements, and Raman spectroscopy. “Rust1LCl” is a low chlorine rustumite Ca10(Si2O7)2(SiO4)(OH1.88F0.12)(Cl1.28,OH0.72) from skarns associated with the Rize batholith near Ikizedere, Turkey. “Rust2F” is a F-bearing rustumite Ca10(Si2O7)2(SiO4)(OH1.13F0.87) (Cl1 96OH0.04) from xenoliths in ignimbrites of the Upper Chegem Caldera, Northern Caucasus, Russia. “Rust3LClF” represents a low-Cl, F-bearing rustumite Ca10(Si2O7)2(SiO4)0.87(H4O4)0.13(OH1.01F0.99) (Cl1.00 OH1.00) from altered merwinite skarns of the Birkhin massif, Baikal Lake area, Eastern Siberia, Russia. Rustumite from Birkhin massif is characterized by a significant hydrogarnet-like or fluorine substitution at the apices of the orthosilicate group, leading to specific atomic displacements. The crystal structures including hydrogen positions have been refined from single-crystal X-ray data to R1 = 0.0205 (Rust1_LCl), R1 = 0.0295 (Rust2_F), and R1 = 0.0243 (Rust3_LCl_F), respectively. Depletion in Cl and replacement by OH is associated with smaller unit-cell dimensions. The substitution of OH by F leads to shorter hydrogen bonds O-H⋯F instead of O-H⋯OH. Raman spectra for all samples have been measured and confirm slight strengthening of the hydrogen bonds with uptake of F.This study discusses the complex crystal chemistry of the skarn mineral rustumite and may provide a wider understanding of the chemical reactions related to contact metamorphism of limestones.

Calculation of crystal optical properties from molecular electron density

We have recently developed a method to obtain distributed atomic polarizabilities adopting a partitioning of the molecular electron density (for example, the Quantum Theory of Atoms in Molecules, [1]), calculated with or without an applied electric field. The procedure [2] allows to obtained atomic polarizability tensors, which are perfectly exportable, because quite representative of an atom in a given functional group. Among the many applications of this idea, the calculation of crystal susceptibility is easily available, either from a rough estimation (the polarizability of the isolated molecule is used) or from a more precise estimation (the polarizability of a molecule embedded in a cluster representing the first coordination sphere is used). Lorentz factor is applied to include the long range effect of packing, which is enhancing the molecular polarizability. Simple properties like linear refractive index or the gyration tensor can be calculated at relatively low costs and with good precision. This approach is particularly useful within the field of crystal engineering of organic/organometallic materials, because it would allow a relatively easy prediction of a property as a function of the packing, thus allowing "reverse crystal engineering". Examples of some amino acid crystals and salts of amino acids [3] will be illustrated, together with other crystallographic or non-crystallographic applications. For example, the induction and dispersion energies of intermolecular interactions could be calculated with superior precision (allowing anisotropic van der Waals interactions). This could allow revision of some commonly misunderstood intermolecular interactions, like the halogen bonding (see for example the recent remarks by Stone or Gilli [4]). Moreover, the chemical reactivity of coordination complexes could be reinvestigated, by coupling the conventional analysis of the electrostatic potential (useful only in the circumstances of hard nucleophilic/electrophilic interaction) with the distributed atomic polarizability. The enhanced reactivity of coordinated organic ligands would be better appreciated. [1] R. F. W. Bader, Atoms in Molecules: A Quantum Theory. Oxford Univ. Press, 1990. [2] A. Krawczuk-Pantula, D. Pérez, K. Stadnicka, P. Macchi, Trans. Amer. Cryst. Ass. 2011, 1-25 [3] A. S. Chimpri1, M. Gryl, L. H.R. Dos Santos1, A. Krawczuk, P. Macchi Crystal Growth & Design, in the press. [4] a) A. J. Stone, J. Am. Chem. Soc. 2013, 135, 7005−7009; b) V. Bertolasi, P. Gilli, G. Gilli Crystal Growth & Design, 2013, 12, 4758-4770.

Anti-PSMA antibody-coupled gold nanorods detection by optical and electron microscopies

While cancer is one of the greatest challenges to public health care, prostate cancer was chosen as cancer model to develop a more accurate imaging assessment than those currently available. Indeed, an efficient imaging technique which considerably improves the sensitivity and specificity of the diagnostic and predicting the cancer behavior would be extremely valuable. The concept of optoacoustic imaging using home-made functionalized gold nanoparticles coupled to an antibody targeting PSMA (prostate specific membrane antigen) was evaluated on different cancer cell lines to demonstrate the specificity of the designed platform. Two commonly used microscopy techniques (indirect fluorescence and scanning electron microscopy) showed their straightforwardness and versatility for the nanoparticle binding investigations regardless the composition of the investigated nanoobjects. Moreover most of the research laboratories and centers are equipped with fluorescence microscopes, so indirect fluorescence using Quantum dots can be used for any active targeting nanocarriers (polymers, ceramics, metals, etc.). The second technique based on backscattered electron is not only limited to gold nanoparticles but also suits for any study of metallic nanoparticles as the electronic density difference between the nanoparticles and binding surface stays high enough. Optoacoustic imaging was finally performed on a 3D cellular model to assess and prove the concept of the developed platform.

Search for Quantum Black Hole Production in High-Invariant-Mass Lepton plus Jet Final States Using pp Collisions at root s=8 TeV and the ATLAS Detector

This Letter presents a search for quantum black-hole production using 20.3 fb(-1) of data collected with the ATLAS detector in pp collisions at the LHC at root s = 8 TeV. The quantum black holes are assumed to decay into a final state characterized by a lepton (electron or muon) and a jet. In either channel, no event with a lepton-jet invariant mass of 3.5 TeV or more is observed, consistent with the expected background. Limits are set on the product of cross sections and branching fractions for the lepton + jet final states of quantum black holes produced in a search region for invariant masses above 1 TeV. The combined 95% confidence level upper limit on this product for quantum black holes with threshold mass above 3.5 TeV is 0.18 fb. This limit constrains the threshold quantum black-hole mass to be above 5.3 TeV in the model considered.

A Cruciform Electron Donor-Acceptor Semiconductor with Solid-State Red Emission: 1D/2D Optical Waveguides and Highly Sensitive/Selective Detection of H2S Gas

In this paper, a new cruciform donor–acceptor molecule 2,2'-((5,5'-(3,7-dicyano-2,6-bis(dihexylamino)benzo[1,2-b:4,5-b']difuran-4,8-diyl)bis(thiophene-5,2-diyl))bis (methanylylidene))dimalononitrile (BDFTM) is reported. The compound exhibits both remarkable solid-state red emission and p-type semiconducting behavior. The dual functions of BDFTM are ascribed to its unique crystal structure, in which there are no intermolecular face-to-face π–π interactions, but the molecules are associated by intermolecular CN…π and H-bonding interactions. Firstly, BDFTM exhibits aggregation-induced emission; that is, in solution, it is almost non-emissive but becomes significantly fluorescent after aggregation. The emission quantum yield and average lifetime are measured to be 0.16 and 2.02 ns, respectively. Crystalline microrods and microplates of BDFTM show typical optical waveguiding behaviors with a rather low optical loss coefficient. Moreover, microplates of BDFTM can function as planar optical microcavities which can confine the emitted photons by the reflection at the crystal edges. Thin films show an air-stable p-type semiconducting property with a hole mobility up to 0.0015 cm2V−1s−1. Notably, an OFET with a thin film of BDFTM is successfully utilized for highly sensitive and selective detection of H2S gas (down to ppb levels).

Fate of accidental symmetries of the relativistic hydrogen atom in a spherical cavity

The non-relativistic hydrogen atom enjoys an accidental SO(4) symmetry, that enlarges the rotational SO(3) symmetry, by extending the angular momentum algebra with the Runge–Lenz vector. In the relativistic hydrogen atom the accidental symmetry is partially lifted. Due to the Johnson–Lippmann operator, which commutes with the Dirac Hamiltonian, some degeneracy remains. When the non-relativistic hydrogen atom is put in a spherical cavity of radius R with perfectly reflecting Robin boundary conditions, characterized by a self-adjoint extension parameter γ, in general the accidental SO(4) symmetry is lifted. However, for R=(l+1)(l+2)a (where a is the Bohr radius and l is the orbital angular momentum) some degeneracy remains when γ=∞ or γ = 2/R. In the relativistic case, we consider the most general spherically and parity invariant boundary condition, which is characterized by a self-adjoint extension parameter. In this case, the remnant accidental symmetry is always lifted in a finite volume. We also investigate the accidental symmetry in the context of the Pauli equation, which sheds light on the proper non-relativistic treatment including spin. In that case, again some degeneracy remains for specific values of R and γ.

Distributed Atomic Polarizabilities of Amino Acids and their Hydrogen-Bonded Aggregates

With the purpose of rational design of optical materials, distributed atomic polarizabilities of amino acid molecules and their hydrogen-bonded aggregates are calculated in order to identify the most efficient functional groups, able to buildup larger electric susceptibilities in crystals. Moreover, we carefully analyze how the atomic polarizabilities depend on the one-electron basis set or the many-electron Hamiltonian, including both wave function and density functional theory methods. This is useful for selecting the level of theory that best combines high accuracy and low computational costs, very important in particular when using the cluster method to estimate susceptibilities of molecular-based materials.

Hydrogen atom in space with a compactified extra dimension and potential defined by Gauss' law

We investigate the consequences of one extra spatial dimension for the stability and energy spectrum of the non-relativistic hydrogen atom with a potential defined by Gauss' law, i.e. proportional to 1 /| x | 2 . The additional spatial dimension is considered to be either infinite or curled-up in a circle of radius R. In both cases, the energy spectrum is bounded from below for charges smaller than the same critical value and unbounded from below otherwise. As a consequence of compactification, negative energy eigenstates appear: if R is smaller than a quarter of the Bohr radius, the corresponding Hamiltonian possesses an infinite number of bound states with minimal energy extending at least to the ground state of the hydrogen atom.

Materials properties from electron density distributions: Molecular magnetism and linear optical properties

The general goal of this thesis is correlating observable properties of organic and metal-organic materials with their ground-state electron density distribution. In a long-term view, we expect to develop empirical or semi-empirical approaches to predict materials properties from the electron density of their building blocks, thus allowing to rationally engineering molecular materials from their constituent subunits, such as their functional groups. In particular, we have focused on linear optical properties of naturally occurring amino acids and their organic and metal-organic derivatives, and on magnetic properties of metal-organic frameworks. For analysing the optical properties and the magnetic behaviour of the molecular or sub-molecular building blocks in materials, we mostly used the more traditional QTAIM partitioning scheme of the molecular or crystalline electron densities, however, we have also investigated a new approach, namely, X-ray Constrained Extremely Localized Molecular Orbitals (XC-ELMO), that can be used in future to extracted the electron densities of crystal subunits. With the purpose of rationally engineering linear optical materials, we have calculated atomic and functional group polarizabilities of amino acid molecules, their hydrogen-bonded aggregates and their metal-organic frameworks. This has enabled the identification of the most efficient functional groups, able to build-up larger electric susceptibilities in crystals, as well as the quantification of the role played by intermolecular interactions and coordinative bonds on modifying the polarizability of the isolated building blocks. Furthermore, we analysed the dependence of the polarizabilities on the one-electron basis set and the many-electron Hamiltonian. This is useful for selecting the most efficient level of theory to estimate susceptibilities of molecular-based materials. With the purpose of rationally design molecular magnetic materials, we have investigated the electron density distributions and the magnetism of two copper(II) pyrazine nitrate metal-organic polymers. High-resolution X-ray diffraction and DFT calculations were used to characterize the magnetic exchange pathways and to establish relationships between the electron densities and the exchange-coupling constants. Moreover, molecular orbital and spin-density analyses were employed to understand the role of different magnetic exchange mechanisms in determining the bulk magnetic behaviour of these materials. As anticipated, we have finally investigated a modified version of the X-ray constrained wavefunction technique, XC-ELMOs, that is not only a useful tool for determination and analysis of experimental electron densities, but also enables one to derive transferable molecular orbitals strictly localized on atoms, bonds or functional groups. In future, we expect to use XC-ELMOs to predict materials properties of large systems, currently challenging to calculate from first-principles, such as macromolecules or polymers. Here, we point out advantages, needs and pitfalls of the technique. This work fulfils, at least partially, the prerequisites to understand materials properties of organic and metal-organic materials from the perspective of the electron density distribution of their building blocks. Empirical or semi-empirical evaluation of optical or magnetic properties from a preconceived assembling of building blocks could be extremely important for rationally design new materials, a field where accurate but expensive first-principles calculations are generally not used. This research could impact the community in the fields of crystal engineering, supramolecular chemistry and, of course, electron density analysis.

Advanced Physical Chemistry Problems (II), Hydrogen Atom and Hyrdrogen Molecular Ion

Slightly advanced problems in Physical Chemistry, herein concerning the H-atom and the Hydrogen Molecular Cation, are presented and discussed.

Specific conductivity and hydrogen ions measured on two ice core (NEEM and NGRIP)

A stratigraphy-based chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core has been derived by transferring the annual layer counted Greenland Ice Core Chronology 2005 (GICC05) and its model extension (GICC05modelext) from the NGRIP core to the NEEM core using 787 match points of mainly volcanic origin identified in the electrical conductivity measurement (ECM) and dielectrical profiling (DEP) records. Tephra horizons found in both the NEEM and NGRIP ice cores are used to test the matching based on ECM and DEP and provide five additional horizons used for the timescale transfer. A thinning function reflecting the accumulated strain along the core has been determined using a Dansgaard-Johnsen flow model and an isotope-dependent accumulation rate parameterization. Flow parameters are determined from Monte Carlo analysis constrained by the observed depth-age horizons. In order to construct a chronology for the gas phase, the ice age-gas age difference (Delta age) has been reconstructed using a coupled firn densification-heat diffusion model. Temperature and accumulation inputs to the Delta age model, initially derived from the water isotope proxies, have been adjusted to optimize the fit to timing constraints from d15N of nitrogen and high-resolution methane data during the abrupt onset of Greenland interstadials. The ice and gas chronologies and the corresponding thinning function represent the first chronology for the NEEM core, named GICC05modelext-NEEM-1. Based on both the flow and firn modelling results, the accumulation history for the NEEM site has been reconstructed. Together, the timescale and accumulation reconstruction provide the necessary basis for further analysis of the records from NEEM.

Seawater carbonate chemistry, calcification rate, oxygen production, maximum quantum yield, symbiont density, chlorophyll concentration and crystal width of Halimeda macroloba, Halimeda cylindracea and Marginopora vertebralis during experiments, 2011

The effects of elevated CO2 and temperature on photosynthesis and calcification in the calcifying algae Halimeda macroloba and Halimeda cylindracea and the symbiont-bearing benthic foraminifera Marginopora vertebralis were investigated through exposure to a combination of four temperatures (28°C, 30°C, 32°C, and 34°C) and four CO2 levels (39, 61, 101, and 203 Pa; pH 8.1, 7.9, 7.7, and 7.4, respectively). Elevated CO2 caused a profound decline in photosynthetic efficiency (FV : FM), calcification, and growth in all species. After five weeks at 34°C under all CO2 levels, all species died. Chlorophyll (Chl) a and b concentration in Halimeda spp. significantly decreased in 203 Pa, 32°C and 34°C treatments, but Chl a and Chl c2 concentration in M. vertebralis was not affected by temperature alone, with significant declines in the 61, 101, and 203 Pa treatments at 28°C. Significant decreases in FV : FM in all species were found after 5 weeks of exposure to elevated CO2 (203 Pa in all temperature treatments) and temperature (32°C and 34°C in all pH treatments). The rate of oxygen production declined at 61, 101, and 203 Pa in all temperature treatments for all species. The elevated CO2 and temperature treatments greatly reduced calcification (growth and crystal size) in M. vertebralis and, to a lesser extent, in Halimeda spp. These findings indicate that 32°C and 101 Pa CO2, are the upper limits for survival of these species on Heron Island reef, and we conclude that these species will be highly vulnerable to the predicted future climate change scenarios of elevated temperature and ocean acidification.

Intermediate band materials for more efficient solar energy use: quantum modelling and experimental realizations

The intermediate band (IB) solar cell (Fig. 1) has been proposed [1] to increase photovoltaic efficiency by a factor above 1.5, based on the absorption of two sub-bandgap photons to promote an electron across the bandgap. To realize this principle, that can be applied also to obtain efficient photocatalysis with sunlight, we proposed in recent years several materials where a metal or heavy element, substituting for an electropositive atom in a known semiconductor that has an appropriate band gap width (around 2 eV), forms inside the gap the partially filled levels needed for this aim