Dissertations / Theses on the topic 'Computational physics|Computational chemistry'

Block copolymers are composed of chemically distinct polymer chains that can be covalently linked in a variety of sequences and architectures. They are ubiquitous as ingredients of consumer products and also have applications in advanced plastics, drug delivery, advanced membranes, and next generation nano-lithographic patterning. The wide spectrum of possible block copolymer applications is a consequence of block copolymer self-assembly into periodic, meso-scale morphologies as a function of varying block composition and architecture in both melt and solution states, and the broad spectrum of physical properties that such mesophases afford.

Materials exploration and discovery has traditionally been pursued through an iterative process between experimental and theoretical/computational collaborations. This process is often implemented in a trial-and-error fashion, and from the computational perspective of generating phase diagrams, usually requires some existing knowledge about the competitive phases for a given system. Self-Consistent Field Theory (SCFT) simulations have proven to be both qualitatively and quantitatively accurate in the determination, or forward mapping, of block copolymer phases of a given system. However, it is possible to miss candidates. This is because SCFT simulations are highly dependent on their initial configurations, and the ability to map phase diagrams requires a priori knowledge of what the competing candidate morphologies are. The unguided search for the stable phase of a block copolymer of a given composition and architecture is a problem of global optimization. SCFT by itself is a local optimization method, so we can combine it with population-based heuristic algorithms geared at global optimization to facilitate forward mapping. In this dissertation, we discuss the development of two such methods: Genetic Algorithm + SCFT (GA-SCFT) and Particle Swarm Optimization + SCFT (PSO-SCFT). Both methods allow a population of configurations to explore the space associated with the numerous states accessible to a block copolymer of a given composition and architecture.

GA-SCFT is a real-space method in which a population of SCFT field configurations “evolves” over time. This is achieved by initializing the population randomly, allowing the configurations to relax to local basins of attraction using SCFT simulations, then selecting fit members (lower free energy structures) to recombine their fields and undergo mutations to generate a new “generation” of structures that iterate through this process. We present results from benchmark testing of this GA-SCFT technique on the canonical AB diblock copolymer melt, for which the theoretical phase diagram has long been established. The GA-SCFT algorithm successfully predicts many of the conventional mesophases from random initial conditions in large, 3-dimensional simulation cells, including hexagonally-packed cylinders, BCC-packed spheres, and lamellae, over a broad composition range and weak to moderate segregation strength. However, the GA-SCFT method is currently not effective at discovery of network phases, such as the Double-Gyroid (GYR) structure.

PSO-SCFT is a reciprocal space approach in which Fourier components of SCFT fields near the principal shell are manipulated. Effectively, PSO-SCFT facilitates the search through a space of reciprocal-space SCFT seeds which yield a variety of morphologies. Using intensive free energy as a fitness metric by which to compare these morphologies, the PSO-SCFT methodology allows us to agnostically identify low-lying competitive and stable morphologies. We present results for applying PSO-SCFT to conformationally symmetric diblock copolymers and a miktoarm star polymer, AB₄, which offers a rich variety of competing sphere structures. Unlike the GA-SCFT method we previously presented, PSO-SCFT successfully predicts the double gyroid morphology in the AB-diblock. Furthermore, PSO-SCFT successfully recovers the A ₁₅ morphology at a composition where it is expected to be stable in the miktoarm system, as well as several competitive metastable candidates, and a new sphere morphology belonging to the hexagonal space group 191, which has not been seen before in polymer systems. Thus, we believe the PSO-SCFT method provides a promising platform for screening for competitive structures in a given block copolymer system.

This dissertation contains investigations based on Molecular Dynamics (MD) of a variety of systems, from small atomic clusters to polymers in solution and in their condensed phases. The overall research is divided in three parts. First, I tested a new thermostat in the literature on the thermal equilibration of a small cluster of Lennard-Jones (LJ) atoms. The proposed thermostat is a Hamiltonian thermostat based on a logarithmic oscillator with the outstanding property that the mean value of its kinetic energy is constant independent of the mass and energy. I inspected several weak-coupling interaction models between the LJ cluster and the logarithmic oscillator in 3D. In all cases I show that this coupling gives rise to a kinetic motion of the cluster center of mass without transferring kinetic energy to the interatomic vibrations. This is a failure of the published thermostat because the temperature of the cluster is mainly due to vibrations in small atomic clusters This logarithmic oscillator cannot be used to thermostat any atomic or molecular system, small or large.

The second part of the dissertation is the investigation of the inherent structure of the polymer polyethylene glycol (PEG) solvated in three different solvents: water, water with 4% ethanol, and ethyl acetate. PEG with molecular weight of 2000 Da (PEG₂₀₀₀) is a polymer with many applications from industrial manufacturing to medicine that in bulk is a paste. However, its structure in very dilute solutions deserved a thorough study, important for the onset of aggregation with other polymer chains. I introduced a modification to the GROMOS 54A7 force field parameters for modeling PEG₂₀₀₀ and ethyl acetate. Both force fields are new and have now been incorporated into the database of known residues in the molecular dynamics package Gromacs. This research required numerous high performance computing MD simulations in the ARGO cluster of GMU for systems with about 100,000 solvent molecules. My findings show that PEG₂₀₀₀ in water acquires a ball-like structure without encapsulating solvent molecules. In addition, no hydrogen bonds were formed. In water with 4% ethanol, PEG₂₀₀₀ acquires also a ball-like structure but the polymer ends fluctuate folding outward and onward, although the general shape is still a compact ball-like structure.

In contrast, PEG₂₀₀₀ in ethyl acetate is quite elongated, as a very flexible spaghetti that forms kinks that unfold to give rise to folds and kinks in other positions along the polymer length. The behavior resembles an ideal polymer in a &thetas; solvent. A Principal Component Analysis (PCA) of the minima composing the inherent structure evidences the presence of two distinct groups of ball-like structures of PEG₂₀₀₀ in water and water with 4% ethanol. These groups give a definite signature to the solvated structure of PEG₂₀₀₀ in these two solvents. In contrast, PCA reveals several groups of avoided states for PEG₂₀₀₀ in ethyl acetate that disqualify the possibility of being an ideal polymer in a &thetas; solvent.

The third part of the dissertation is a work in progress, where I investigate the condensed phase of PEG₂₀₀₀ and study the interface between the condensed phase and the three different solvents under study. With a strategy of combining NPT MD simulations at different temperatures and pressures, PEG ₂₀₀₀ condensed phase displays the experimental density within a 1% discrepancy at 300 K and 1 atm. This is a very encouraging result on this ongoing project.

We use non-adiabatic ab initio molecular dynamics to study the influence of substituent side groups on the photoactive unit (Z)-hexa-1,3,5-triene (HT). The Time-Dependent Density Functional Theory Surface Hopping method (TDDFT-SH) is used to investigate the influence of substituted isopropyl and methyl groups on the excited state dynamics. The 1,4 and 2,5-substituted molecules are simulated: 2,5-dimethylhexa-1,3,5-triene (DMHT), 2-isopropyl-5-methyl-1,3,5-hexatriene (2,5-IMHT), 3,7-dimethylocta-1,3,5-triene (1,4-IMHT), and 2,5-diisopropyl-1,3,5-hexatriene (DIHT). We find that HT and 1,4-IMHT have the lowest ring-closing branching ratios of 5.3% and 1.0%, respectively. For the 2,5-substituted derivatives, the branching ratio increases with increasing size of the substituents, exhibiting yields of 9.78%, 19%, and 24% for DMHT, 2,5-IMHT, and DIHT, respectively. The reaction channels are shown to prefer certain conformation configurations at excitation, where the ring-closing reaction tends to originate from the gauche-Z-gauche (gZg) rotamer almost exclusively. In addition, there is a conformational dependency on absorption, gZg conformers have on average lower S₁ ← S₀ excitation energies that the other rotamers. Furthermore, we develop a method to calculate a predicted quantum yield that is in agreement with the wavelength-dependence observed in experiment for DMHT. In addition, the quantum yield method also predicts DIHT to have the highest CHD yield of 0.176 at 254 nm and 0.390 at 290 nm.

Additionally, we study the vitamin D derivative Tachysterol (Tachy) which exhibits similar photochemical properties as HT and its derivatives. We find the reaction channels of Tachy also have a conformation dependency, where the reactive products toxisterol-D1 (2.3%), previtamin D (1.4%) and cyclobutene toxisterol (0.7%) prefer cEc, cEt, and tEc configurations at excitation, leaving the tEt completely non-reactive. The rotamers similarly have a dependence on absorption as well, where the cEc configuration has the lowest energy S ₁ ← S₀ excitation of the rotamers. The wavelength dependence of the rotamers should lead to selective properties of these molecules at excitation. An excitation to the red-shifted side of the maximum absorption peak will on average lead to excitations of the gZg rotamers more exclusively.

This thesis presents briefly the application of quantum mechanics on systems ofchemical interest, i.e., the field of quantum chemistry and computational chemistry.The molecules of the ten-electron series, hydrogen fluoride, water, ammonia,methane and neon, are taken as computational examples. Some applications ofquantum chemistry are then shown on these systems, with emphasis on the natureof the molecular bonds. Conceptual methods of chemistry and theoreticalchemistry for these systems are shown to be valid with some restrictions, as theseinterpretations does not represent physically measurable entities.The orbitals and orbital energies of neon is studied, the binding van der Waalsinteractionresulting in a Ne2 molecule is studied with a theoretical bond lengthof 3.23 °A and dissociation energy of 81.75 μEh. The equilibrium geometries ofFH, H2O, NH3 and CH4 are studied and the strength and character of the bondsinvolved evaluated using bond order, dipole moment, Mulliken population analysisand L¨owdin population analysis. The concept of electronegativity is studied in thecontext of electron transfer. Lastly, the barrier of inversion for NH3 is studied, withan obtained barrier height of 8.46 mEh and relatively constant electron transfer.