Simulation of Smart Materials

The aim of this thesis is the elucidation of structure-properties relationship of molecular semiconductors for electronic devices. This involves the use of a comprehensive set of simulation techniques, ranging from quantum-mechanical to numerical stochastic methods, and also the development of ad-hoc computational tools. In more detail, the research activity regarded two main topics: the study of electronic properties and structural behaviour of liquid crystalline (LC) materials based on functionalised oligo(p-phenyleneethynylene) (OPE), and the investigation on the electric field effect associated to OFET operation on pentacene thin film stability. In this dissertation, a novel family of substituted OPE liquid crystals with applications in stimuli-responsive materials is presented. In more detail, simulations can not only provide evidence for the characterization of the liquid crystalline phases of different OPEs, but elucidate the role of charge transfer states in donor-acceptor LCs containing an endohedral metallofullerene moiety. Such systems can be regarded as promising candidates for organic photovoltaics. Furthermore, exciton dynamics simulations are performed as a way to obtain additional information about the degree of order in OPE columnar phases. Finally, ab initio and molecular mechanics simulations are used to investigate the influence of an applied electric field on pentacene reactivity and stability. The reaction path of pentacene thermal dimerization in the presence of an external electric field is investigated; the results can be related to the fatigue effect observed in OFETs, that show significant performance degradation even in the absence of external agents. In addition to this, the effect of the gate voltage on a pentacene monolayer are simulated, and the results are then compared to X-ray diffraction measurements performed for the first time on operating OFETs.

Charge induced enhancement of adsorption for hydrogen storage materials

The rising concerns about environmental pollution and global warming have facilitated research interest in hydrogen energy as an alternative energy source. To apply hydrogen for transportations, several issues have to be solved, within which hydrogen storage is the most critical problem. Lots of materials and devices have been developed; however, none is able to meet the DOE storage target. The primary issue for hydrogen physisorption is a weak interaction between hydrogen and the surface of solid materials, resulting negligible adsorption at room temperature. To solve this issue, there is a need to increase the interaction between the hydrogen molecules and adsorbent surface. In this study, intrinsic electric dipole is investigated to enhance the adsorption energy. The results from the computer simulation of single ionic compounds with hydrogen molecules to form hydrogen clusters showed that electrical charge of substances plays an important role in generation of attractive interaction with hydrogen molecules. In order to further examine the effects of static interaction on hydrogen adsorption, activated carbon with a large surface area was impregnated with various ionic salts including LiCl, NaCl, KCl, KBr, and NiCl and their performance for hydrogen storage was evaluated by using a volumetric method. Corresponding computer simulations have been carried out by using DFT (Density Functional Theory) method combined with point charge arrays. Both experimental and computational results prove that the adsorption capacity of hydrogen and its interaction with the solid materials increased with electrical dipole moment. Besides the intrinsic dipole, an externally applied electric field could be another means to enhance hydrogen adsorption. Hydrogen adsorption under an applied electric field was examined by using porous nickel foil as electrodes. Electrical signals showed that adsorption capacity increased with the increasing of gas pressure and external electric voltage. Direct measurement of the amount of hydrogen adsorption was also carried out with porous nickel oxides and magnesium oxides using the piezoelectric material PMN-PT as the charge supplier due to the pressure. The adsorption enhancement from the PMN-PT generated charges is obvious at hydrogen pressure between 0 and 60 bars, where the hydrogen uptake is increased at about 35% for nickel oxide and 25% for magnesium oxide. Computer simulation reveals that under the external electric field, the electron cloud of hydrogen molecules is pulled over to the adsorbent site and can overlap with the adsorbent electrons, which in turn enhances the adsorption energy Experiments were also carried out to examine the effects of hydrogen spillover with charge induced enhancement. The results show that the overall storage capacity in nickel oxide increased remarkably by a factor of 4.

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.

Response of liquid xenon to Compton electrons down to 1.5 keV

The response of liquid xenon to low-energy electronic recoils is relevant in the search for dark-matter candidates which interact predominantly with atomic electrons in the medium, such as axions or axionlike particles, as opposed to weakly interacting massive particles which are predicted to scatter with atomic nuclei. Recently, liquid-xenon scintillation light has been observed from electronic recoils down to 2.1 keV, but without applied electric fields that are used in most xenon dark-matter searches. Applied electric fields can reduce the scintillation yield by hindering the electron-ion recombination process that produces most of the scintillation photons. We present new results of liquid xenon's scintillation emission in response to electronic recoils as low as 1.5 keV, with and without an applied electric field. At zero field, a reduced scintillation output per unit deposited energy is observed below 10 keV, dropping to nearly 40% of its value at higher energies. With an applied electric field of 450 V/cm, we observe a reduction of the scintillation output to about 75% relative to the value at zero field. We see no significant energy dependence of this value between 1.5 and 7.8 keV. With these results, we estimate the electronic-recoil energy thresholds of ZEPLIN-III, XENON10, XENON100, and XMASS to be 2.8, 2.5, 2.3, and 1.1 keV, respectively, validating their excellent sensitivity to low-energy electronic recoils.

Computer simulations of sample preconcentration in carrier-free systems and isoelectric focusing in microchannels using simple ampholytes.

In this work, electrophoretic preconcentration of protein and peptide samples in microchannels was studied theoretically using the 1D dynamic simulator GENTRANS, and experimentally combined with MS. In all configurations studied, the sample was uniformly distributed throughout the channel before power application, and driving electrodes were used as microchannel ends. In the first part, previously obtained experimental results from carrier-free systems are compared to simulation results, and the effects of atmospheric carbon dioxide and impurities in the sample solution are examined. Simulation provided insight into the dynamics of the transport of all components under the applied electric field and revealed the formation of a pure water zone in the channel center. In the second part, the use of an IEF procedure with simple well defined amphoteric carrier components, i.e. amino acids, for concentration and fractionation of peptides was investigated. By performing simulations a qualitative description of the analyte behavior in this system was obtained. Neurotensin and [Glu1]-Fibrinopeptide B were separated by IEF in microchannels featuring a liquid lid for simple sample handling and placement of the driving electrodes. Component distributions in the channel were detected using MALDI- and nano-ESI-MS and data were in agreement with those obtained by simulation. Dynamic simulations are demonstrated to represent an effective tool to investigate the electrophoretic behavior of all components in the microchannel.

Broken symmetry and critical transport properties of random metals

Recent experimental data on the conductivity σ+(T), T → 0, on the metallic side of the metal–insulator transition in ideally random (neutron transmutation-doped) 70Ge:Ga have shown that σ+(0) ∝ (N − Nc)μ with μ = ½, confirming earlier ultra-low-temperature results for Si:P. This value is inconsistent with theoretical predictions based on diffusive classical scaling models, but it can be understood by a quantum-directed percolative filamentary amplitude model in which electronic basis states exist which have a well-defined momentum parallel but not normal to the applied electric field. The model, which is based on a new kind of broken symmetry, also explains the anomalous sign reversal of the derivative of the temperature dependence in the critical regime.

The role of London forces in defining noncentrosymmetric order of high dipole moment–high hyperpolarizability chromophores in electrically poled polymeric thin films

Graphs of second harmonic generation coefficients and electro-optic coefficients (measured by ellipsometry, attenuated total reflection, and two-slit interference modulation) as a function of chromophore number density (chromophore loading) are experimentally observed to exhibit maxima for polymers containing chromophores characterized by large dipole moments and polarizabilities. Modified London theory is used to demonstrated that this behavior can be attributed to the competition of chromophore-applied electric field and chromophore–chromophore electrostatic interactions. The comparison of theoretical and experimental data explains why the promise of exceptional macroscopic second-order optical nonlinearity predicted for organic materials has not been realized and suggests routes for circumventing current limitations to large optical nonlinearity. The results also suggest extensions of measurement and theoretical methods to achieve an improved understanding of intermolecular interactions in condensed phase materials including materials prepared by sequential synthesis and block copolymer methods.

Intermolecular and Intramolecular Charge Transfer in Functional Molecular Materials

This thesis is devoted to the investigation of inter and intramolecular charge transfer (CT) in molecular functional materials and specifically organic dyes and CT crystals. An integrated approach encompassing quantum-chemical calculations, semiempirical tools, theoretical models and spectroscopic measurements is applied to understand structure-property relationships governing the low-energy physics of these materials. Four main topics were addressed: 1) Spectral properties of organic dyes. Charge-transfer dyes are constituted by electron donor (D) and electron acceptor (A) units linked through bridge(s) to form molecules with different symmetry and dimensionality. Their low-energy physics is governed by the charge resonance between D and A groups and is effectively described by a family of parametric Hamiltonians known as essential-state models. These models account for few electronic states, corresponding to the main resonance structures of the relevant dye, leading to a simple picture that is completed introducing the coupling of the electronic system to molecular vibrations, treated in a non-adiabatic way, and an effective classical coordinate, describing polar solvation. In this work a specific essential-state model was proposed and parametrized for the dye Brilliant Green. The central issue in this work has been the definition of the diabatic states, a not trivial task for a multi-branched chromophore. In a second effort, we have used essential-state models for the description of the early-stage dynamics of excited states after ultrafast excitation. Crucial to this work is the fully non-adiabatic treatment of the coupled electronic and vibrational motion, allowing for a reliable description of the dynamics of systems showing a multistable, broken-symmetry excited state. 2) Mixed-stack CT salts. Mixed-stack (MS) CT crystals are an interesting class of multifunctional molecular materials, where D and A molecules arrange themselves to form stacks, leading to delocalized electrons in one dimension. The interplay between the intermolecular CT, electrostatic interactions, lattice phonons and molecular vibrations leads to intriguing physical properties that include (photoinduced) phase transitions, multistability, antiferromagnetism, ferroelectricity and potential multiferroicity. The standard microscopic model to describe this family of materials is the Modified Hubbard model accounting for electron-phonon coupling (Peierls coupling), electron-molecular vibrations coupling (Holstein coupling) and electrostatic interactions. We adopt and validate a method, based on DFT calculations on dimeric DA structures, to extract relevant model parameters. The approach offers a powerful tool to shed light on the complex physics of MS-CT salts. 3) Charge transfer in organic radical dipolar dyes. In collaboration with the group of Prof. Jaume Veciana (ICMAB- Barcellona), we have studied spectral properties of a special class of CT dyes with D-bridge-A structure where the acceptor group is a stable radical (of the perchlorotriphenylmethyl, PTM, family), leading to an open-shell CT dyes. These materials are of interest since they associate the electronic and optical properties of CT dyes with magnetic properties from the unpaired electron. The first effort was devoted to the parametrization of the relevant essential-state model. Two strategies were adopted, one based on the calculation of the low-energy spectral properties, the other based on the variation of ground state properties with an applied electric field. 4) The spectral properties of organic nanoparticles based on radical species are investigated in collaboration with Dr. I. Ratera (ICMAB- Barcellona). Intriguing spectroscopic behavior was observed pointing to the presence of excimer states. In an attempt to rationalize these findings, extensive calculations (TD-DFT and ZINDO) were performed. The results for the isolated dyes are validated against experimental spectra in solution. To address intermolecular interactions we studied dimeric structures in the gas phase, but the preliminary results obtained do not support excimer formation.

Terahertz Response of Excitons in Nanorod Heterostructures

Quantum-confined systems are one of the most promising ways to enable us to control a material's interactions with light. Nanorods in particular offer the right dimensions for exploring and manipulating the terahertz region of the spectrum. In this thesis, we model excitons confined inside a nanorod using the envelope function approximation. A region-matching transfer matrix method allows us to simulate excitonic states inside arbitrary heterostructures grown along the length of the rod. We apply the method to colloidal CdSe rods 70 nm in length and under 10 nm in diameter, capped with ligands of DDPA and pyridine. We extend past studies on these types of rods by taking into account their dielectric permittivity mismatch. Compared to previous calculations and experimentally measured terahertz absorption, we predict a higher energy main 1S$z$ to 2P$z$ transition peak. This indicates that the rods are likely larger in diameter than previously thought. We also investigate a nanorod with GaAs/Al$_{0.3}$Ga$_{0.7}$As coupled double dots. The excitonic transitions were found to be manipulable by varying the strength of an applied electric field. We employ quasi-static state population distributions to simulate the effects of exciton relaxation from optically active states to dim ground states. A critical value of the applied field, corresponding to the exciton binding energy of ~18 meV, was found to dramatically alter the terahertz absorption due to state mixing. Above this critical field, more nuanced shifts in transition energies were observed, and gain from radiative relaxation to the ground state is predicted.

Numerical method for determination of the NMR frequency of the single-qubit operation in a silicon-based solid-state quantum computer

A numerical method is introduced to determine the nuclear magnetic resonance frequency of a donor (P-31) doped inside a silicon substrate under the influence of an applied electric field. This phosphorus donor has been suggested for operation as a qubit for the realization of a solid-state scalable quantum computer. The operation of the qubit is achieved by a combination of the rotation of the phosphorus nuclear spin through a globally applied magnetic field and the selection of the phosphorus nucleus through a locally applied electric field. To realize the selection function, it is required to know the relationship between the applied electric field and the change of the nuclear magnetic resonance frequency of phosphorus. In this study, based on the wave functions obtained by the effective-mass theory, we introduce an empirical correction factor to the wave functions at the donor nucleus. Using the corrected wave functions, we formulate a first-order perturbation theory for the perturbed system under the influence of an electric field. In order to calculate the potential distributions inside the silicon and the silicon dioxide layers due to the applied electric field, we use the multilayered Green's functions and solve an integral equation by the moment method. This enables us to consider more realistic, arbitrary shape, and three-dimensional qubit structures. With the calculation of the potential distributions, we have investigated the effects of the thicknesses of silicon and silicon dioxide layers, the relative position of the donor, and the applied electric field on the nuclear magnetic resonance frequency of the donor.

Investigation of photorefractive subharmonics in the absence of wavemixing

Using a new optical configuration free from the influence of photorefractive optical nonlinearity, we investigate the main characteristics of the spatial subharmonic K/2 excited in a Bi12SiO20 crystal by a light-intensity pattern with wave vector K and frequency O. It is shown that in a large region of intensity and applied electric field the optimum value O of the frequency corresponds to the conditions of parametric excitation of the weakly damped eigenmodes of the medium: the space-charge waves. The threshold and above-threshold characteristics of the subharmonic regime are in good agreement with the theory.

Investigation of photorefractive subharmonics in the absence of wavemixing

Using a new optical configuration free from the influence of photorefractive optical nonlinearity, we investigate the main characteristics of the spatial subharmonic K/2 excited in a Bi12SiO20 crystal by a light-intensity pattern with wave vector K and frequency O. It is shown that in a large region of intensity and applied electric field the optimum value O of the frequency corresponds to the conditions of parametric excitation of the weakly damped eigenmodes of the medium: the space-charge waves. The threshold and above-threshold characteristics of the subharmonic regime are in good agreement with the theory.

Two-wave mixing in BTO in the presence of detuning

The gain of a weak signal beam is measured in the usual two-wave mixing configuration in a photorefractive BTO crystal by varying a set of parameters, namely the applied electric field, the detuning, the beam ratio and the input intensity of the pump beam. All the experimental results are shown to be in good agreement with a theory based on that of Kukhtarev et al.

Two-wave mixing in BTO in the presence of detuning

The gain of a weak signal beam is measured in the usual two-wave mixing configuration in a photorefractive BTO crystal by varying a set of parameters, namely the applied electric field, the detuning, the beam ratio and the input intensity of the pump beam. All the experimental results are shown to be in good agreement with a theory based on that of Kukhtarev et al.

A Stochastic Approach for Investigation Ultrafast Phenomena in Semiconductors

2002 Mathematics Subject Classification: 65C05