Vibration and dissociation coupling behind strong shock waves. Paper presented at the Symposium on Dynamics of Manned Lifting Planetary Entry,

Cover title.

Base-induced decomposition of alkyl hydroperoxides in the gas phase. Part 3. Kinetics and dynamics in HO-+CH3OOH, C2H5OOH, and tert-C4H9OOH reactions

The E-CO(2) elimination reactions of alkyl hydroperoxides proceed via abstraction of an (x-hydrogen by a base: X- + (RRHCOOH)-R-1-H-2 -> HX + (RRC)-R-1-C-2=O + HO-. Efficiencies and product distributions for the reactions of the hydroxide anion with methyl, ethyl, and tert-butyl hydroperoxides are studied in the gas phase. On the basis of experiments using three isotopic analogues, HO- + CH3OOH, HO- + CD3OOH, and H18O- + CH3OOH. the overall intrinsic reaction efficiency is determined to be 80% or greater. The E(CO)2 decomposition is facile for these methylperoxide reactions, and predominates over competing proton transfer at the hydroperoxide moiety. The CH3CH2OOH reaction displays a similar E(CO)2 reactivity, whereas proton transfer and the formation of HOO- are the exclusive pathways observed for (CH3)(3)COOH, which has no (x-hydrogen. All results are consistent with the E-CO(2) mechanism, transition state structure, and reaction energy diagrams calculated using the hybrid density functional B3LYP approach. Isotope labeling for HO- + CH3OOH also reveals some interaction between H2O and HO- within the E(CO)2 product complex [H2O center dot center dot center dot CH2=O center dot center dot center dot HO-]. There is little evidence, however. for the formation of the most exothermic products H2O + CH2(OH)O-, which would arise from nuclephilic condensation of CH2=O and HO-. The results suggest that the product dynamics are not totally statistical but are rather direct after the E-CO(2) transition state. The larger HO- + CH3CH2OOH system displays more statistical behavior during complex dissociation.

Analysis of dynamics and mechanism of ligand binding to Artocarpus integrifolia agglutinin. A 13C and 19F NMR study

Binding of 13C-labeled N-acetylgalactosamine (13C-GalNAc) and N-trifluoroacetylgalactosamine (19F-GalNAc) to Artocarpus integrifolia agglutinin has been studied using 13C and 19F nuclear magnetic resonance spectroscopy, respectively. Binding of these saccharides resulted in broadening of the resonances, and no change in chemical shift was observed, suggesting that the alpha- and beta-anomers of 13C-GalNAc and 19F-GalNAc experience a magnetically equivalent environment in the lectin combining site. The alpha- and beta-anomers of 13C-GalNAc and 19F-GalNAc were found to be in slow exchange between free and protein bound states. Binding of 13C-GalNAc was studied as a function of temperature. From the temperature dependence of the line broadening, the thermodynamic and kinetic parameters were evaluated. The association rate constants obtained for the alpha-anomers of 13C-GalNAc and 19F-GalNAc (k+1 = 1.01 x 10(5) M-1.s-1 and 0.698 x 10(5) M-1.s-1, respectively) are in close agreement with those obtained for the corresponding beta-anomers (k+1 = 0.95 x 10(5) M-1.s-1 and 0.65 x 10(5) M-1.s-1, respectively), suggesting that the two anomers bind to the lectin by a similar mechanism. In addition these values are several orders of magnitude slower than those obtained for diffusion controlled processes. The dissociation rate constants obtained are 49.9, 56.9, 42, and 43 s-1, respectively, for the alpha- and beta-anomers of 13C-GalNAc and 19F-GalNAc. A two-step mechanism has been proposed for the interaction of 13C-GalNAc and 19F-GalNAc with A. integrifolia lectin in view of the slow association rates and high activation entropies. The thermodynamic parameters obtained for the association and dissociation reactions suggest that the binding process is entropically favored and that there is a small enthalpic contribution.

Photodissociation of isomeric dichloroethylenes in the ultraviolet: Effect of the second chlorine atom substitution on the dynamics

We report the Cl-*(P-2(1/2)) production dynamics in the near-UV dissociation of three isomers (cis-, gem-, and trans-) of dichloroethylene using the conventional resonance enhanced multiphoton ionization technique. Substantial amounts of Cl-* are produced in the wavelength range 222-304 nm. The Cl-* quantum yield (phi(*)) i maximum at 304 nm for all the isomers and phi(*)(cis) is markedly higher than phi(*)(gem) and phi(*)(trans) except at 222 nm. Existence of both direct and indirect dissociation pathways at these wavelengths complicates the Cl* production dynamics. The higher value of phi(*)(cis) originates from a large contribution from direct dissociation via the (n, sigma(*)) state.

Dynamics of I-asterisk(P-2(1/2)) production inthe ultraviolet photodissociation of alkyn iodides

The relative quantum yields, phi*, for the production of I*(P-2(1/2)) at 266, 280, and similar to 305 nm are reported for a series of primary alkyl iodides using the technique of two-photon laser-induced fluorescence for the detection of I(P-2(3/2)) and I*(P-2(1/2)) atoms. Results are analyzed by invoking the impulsive energy disposal model, which summarizes the dynamics of dissociation as a single parameter. Comparison of our data with those calculated by a more sophisticated time-dependent quantum mechanical model is also made. Near the red edge of the alkyl iodide A band, absorption contribution from the (3)Q(1) state is important and the dynamics near the (3)Q(0)-(1)Q(1) curve-crossing region seem to be influenced by the kinematics of the dissociation process

Dynamics of crystal size distributions with size-dependent rates

The growth and dissolution dynamics of nonequilibrium crystal size distributions (CSDs) can be determined by solving the governing population balance equations (PBEs) representing reversible addition or dissociation. New PBEs are considered that intrinsically incorporate growth dispersion and yield complete CSDs. We present two approaches to solving the PBEs, a moment method and a numerical scheme. The results of the numerical scheme agree with the moment technique, which can be solved exactly when powers on mass-dependent growth and dissolution rate coefficients are either zero or one. The numerical scheme is more general and can be applied when the powers of the rate coefficients are non-integers or greater than unity. The influence of the size dependent rates on the time variation of the CSDs indicates that as equilibrium is approached, the CSDs become narrow when the exponent on the growth rate is less than the exponent on the dissolution rate. If the exponent on the growth rate is greater than the exponent on the dissolution rate, then the polydispersity continues to broaden. The computation method applies for crystals large enough that interfacial stability issues, such as ripening, can be neglected. (C) 2002 Elsevier Science B.V. All rights reserved.

Photodissociation dynamics of CH2ICl at 222, 236, 266, 280, and-304 nm

Dynamics of I*(P-2(1/2)) formation from CH2ICl dissociation has-been investigated at five different ultraviolet excitation wavelengths, e.g., 222, 236, 266, 280, and similar to304 nm. The quantum yield of I*((2)p(1/2)) production, phi*, has been measured by monitoring nascent I(P-2(3/2)) and I* concentrations using a resonance enhanced multiphoton ionization detection scheme. The measured quantum yield as a function of excitation energy follows the same trend as that of methyl iodide except at 236 run. The photodissociation dynamics of CH2ICl also involves three upper states similar to methyl iodide, and a qualitative correlation diagram has been constructed to account for the observed quantum yield. From the difference in behavior at 236 nm, it appears that the crossing region between the two excited states ((3)Q(0) and (1)Q(1)) is located near the exit valley away from the Franck Condon excitation region. The B- and C-band transitions do not participate in the dynamics, and the perturbation of the methyl iodide states due to Cl-I interaction is relatively weak at the photolysis wavelengths employed in this investigation.

Quaternary association in beta-prism I fold plant lectins: Insights from X-ray crystallography, modelling and molecular dynamics

Dimeric banana lectin and calsepa, tetrameric artocarpin and octameric heltuba are mannose-specific beta-prism I fold lectins of nearly the same tertiary structure. MD simulations on individual subunits and the oligomers provide insights into the changes in the structure brought about in the protomers on oligomerization, including swapping of the N-terminal stretch in one instance. The regions that undergo changes also tend to exhibit dynamic flexibility during MD simulations. The internal symmetries of individual oligomers are substantially retained during the calculations. Energy minimization and simulations were also carried out on models using all possible oligomers by employing the four different protomers. The unique dimerization pattern observed in calsepa could be traced to unique substitutions in a peptide stretch involved in dimerization. The impossibility of a specific mode of oligomerization involving a particular protomer is often expressed in terms of unacceptable steric contacts or dissociation of the oligomer during simulations. The calculations also led to a rationale for the observation of a heltuba tetramer in solution although the lectin exists as an octamer in the crystal, in addition to providing insights into relations among evolution, oligomerization and ligand binding.

An Antitubulin Agent BCFMT Inhibits Proliferation of Cancer Cells and Induces Cell Death by Inhibiting Microtubule Dynamics

Using cell based screening assay, we identified a novel anti-tubulin agent (Z)-5-((5-(4-bromo-3-chlorophenyl)furan-2-yl)methylene)-2-thioxothiazoli din-4-one (BCFMT) that inhibited proliferation of human cervical carcinoma (HeLa) (IC50, 7.2 +/- 1.8 mu M), human breast adenocarcinoma (MCF-7) (IC50, 10.0 +/- 0.5 mu M), highly metastatic breast adenocarcinoma (MDA-MB-231) (IC50, 6.0 +/- 1 mu M), cisplatin-resistant human ovarian carcinoma (A2780-cis) (IC50, 5.8 +/- 0.3 mu M) and multi-drug resistant mouse mammary tumor (EMT6/AR1) (IC50, 6.5 +/- 1 mu M) cells. Using several complimentary strategies, BCFMT was found to inhibit cancer cell proliferation at G2/M phase of the cell cycle apparently by targeting microtubules. In addition, BCFMT strongly suppressed the dynamics of individual microtubules in live MCF-7 cells. At its half maximal proliferation inhibitory concentration (10 mu M), BCFMT reduced the rates of growing and shortening phases of microtubules in MCF-7 cells by 37 and 40%, respectively. Further, it increased the time microtubules spent in the pause (neither growing nor shortening detectably) state by 135% and reduced the dynamicity (dimer exchange per unit time) of microtubules by 70%. In vitro, BCFMT bound to tubulin with a dissociation constant of 8.3 +/- 1.8 mu M, inhibited tubulin assembly and suppressed GTPase activity of microtubules. BCFMT competitively inhibited the binding of BODIPY FL-vinblastine to tubulin with an inhibitory concentration (K-i) of 5.2 +/- 1.5 mu M suggesting that it binds to tubulin at the vinblastine site. In cultured cells, BCFMT-treatment depolymerized interphase microtubules, perturbed the spindle organization and accumulated checkpoint proteins (BubR1 and Mad2) at the kinetochores. BCFMT-treated MCF-7 cells showed enhanced nuclear accumulation of p53 and its downstream p21, which consequently activated apoptosis in these cells. The results suggested that BCFMT inhibits proliferation of several types of cancer cells including drug resistance cells by suppressing microtubule dynamics and indicated that the compound may have chemotherapeutic potential.

Ab Initio Molecular Dynamics Simulations of Amino Acids in Aqueous Solutions: Estimating pK(a) Values from Metadynamics Sampling

Changes in the protonation and deprotonation of amino acid residues in proteins play a key role in many biological processes and pathways. Here, we report calculations of the free-energy profile for the protonation deprotonation reaction of the 20 canonical alpha amino acids in aqueous solutions using ab initio Car-Parrinello molecular dynamics simulations coupled with metad-ynamics sampling. We show here that the calculated change in free energy of the dissociation reaction provides estimates of the multiple pK(a) values of the amino acids that are in good agreement with experiment. We use the bond-length-dependent number of the protons coordinated to the hydroxyl oxygen of the carboxylic and the amine groups as the collective variables to explore the free-energy profiles of the Bronsted acid-base chemistry of amino acids in aqueous solutions. We ensure that the amino acid undergoing dissociation is solvated by at least three hydrations shells with all water molecules included in the simulations. The method works equally well for amino acids with neutral, acidic and basic side chains and provides estimates of the multiple pK(a) values with a mean relative error, with respect to experimental results, of 0.2 pK(a) units.

Forced extension of P-selectin construct using steered molecular dynamics

P-selectin, a 70-nm-long cellular adhesive molecule, possesses elastic and extensible properties when neutrophils roll over the activated endotheliam of blood vessel in inflammatory reaction. Transient formation and dissociation of P-selectin/ligand bond on applied force of blood flow induces the extension of P-selectin and relevant ligands. Steered molecular dynamics simulations were performed to stretch a single P-selectin construct consisting of a lectin (Lec) domain and an epithelial growth factor (EGF)-like domain, where P-selectin construct was forced to extend in water with pulling velocities of 0.005-0.05 nm/ps and with constant forces of 1000-2500 pN respectively. Resulting force-extension profiles exhibited a dual-peak pattern on various velocities, while both plateaus and shoulders appeared in the extension-time profiles on various forces. The force or extension profiles along stretching pathways were correlated to the conformational changes, suggesting that the structural collapses of P-selectin Lec/EGF domains were mainly attributed to the burst of hydrogen bonds within the major beta sheet of EGF domain and the disruptions of two hydrophobic cores of Lee domain. This work furthers the understanding of forced dissociation of P-selectin/ligand bond.

Molecular dynamics simulation of interaction of a dislocation array from a crack tip with grain boundaries

The interaction of a dislocation array emitted from a crack tip under mode II loading with asymmetric tilt grain boundaries (GBs) is analysed by the molecular dynamics method. The GBs can generally be described by planar and linear matching zones and unmatching zones. All GBs are observed to emit dislocations. The GBs migrated easily due to their planar and linear matching structure and asymmetrical type. The diffusion induced by stress concentration is found to promote the GB migration. The transmissions of dislocations are either along the matched plane or along another plane depending on tilt angle theta. Alternate processes of stress concentration and stress relaxation take place ahead of the pileup. The stress concentration can be released either by transmission of dislocations, by atom diffusion along GBs, or by migration of GBs by formation of twinning bands. The simulated results also unequivocally demonstrate two processes, i.e. asymmetrical GBs evolving into symmetrical ones and unmatching zones evolving into matching ones during the loading process.

Photodissociation and reaction dynamics of chlorine-containing species important in stratospheric ozone chemistry

Chlorine oxide species have received considerable attention in recent years due to their central role in the balance of stratospheric ozone. Many questions pertaining to the behavior of such species still remain unanswered and plague the ability of researchers to develop accurate chemical models of the stratosphere. Presented in this thesis are three experiments that study various properties of some specific chlorine oxide species.

In the first chapter, the reaction between ClONO_2 and protonated water clusters is investigated to elucidate a possible reaction mechanism for the heterogeneous reaction of chlorine nitrate on ice. The ionic products were various forms of protonated nitric acid, NO_2 +(H_20)_m, m = 0, 1, 2. These products are analogous to products previously reported in the literature for the neutral reaction occurring on ice surfaces. Our results support the hypothesis that the heterogeneous reaction is acid-catalyzed.

In the second chapter, the photochemistry of ClONO_2 was investigated at two wavelengths, 193 and 248 nm, using the technique of photofragmentation translational spectroscopy. At both wavelengths, the predominant dissociation pathways were Cl + NO_3 and ClO + NO_2. Channel assignments were confirmed by momentum matching the counterfragments from each channel. A one-dimensional stratospheric model using the new 248 nm branching ratio determined how our results would affect the predicted Cl_x and NO_x partitioning in the stratosphere.

Chapter three explores the photodissociation dynamics of Cl_2O at 193, 248 and 308 nm. At 193 nm, we found evidence for the concerted reaction channel, Cl_2 + O. The ClO + Cl channel was also accessed, however, the majority of the ClO fragments were formed with sufficient internal energies for spontaneous secondary dissociation to occur. At 248 and 308 nm, we only observed only the ClO + Cl channel. . Some of the ClO formed at 248 nm was formed internally hot and spontaneously dissociated. Bimodal translational energy distributions of the ClO and Cl products indicate two pathways leading to the same product exist.

Appendix A, B and C discuss the details of data analysis techniques used in Chapters 1 and 2. The development of a molecular beam source of ClO dimer is presented in Appendix D.

First principles based multiparadigm modeling of electronic structures and dynamics

Electronic structures and dynamics are the key to linking the material composition and structure to functionality and performance.

An essential issue in developing semiconductor devices for photovoltaics is to design materials with optimal band gaps and relative positioning of band levels. Approximate DFT methods have been justified to predict band gaps from KS/GKS eigenvalues, but the accuracy is decisively dependent on the choice of XC functionals. We show here for CuInSe₂ and CuGaSe₂, the parent compounds of the promising CIGS solar cells, conventional LDA and GGA obtain gaps of 0.0-0.01 and 0.02-0.24 eV (versus experimental values of 1.04 and 1.67 eV), while the historically first global hybrid functional, B3PW91, is surprisingly the best, with band gaps of 1.07 and 1.58 eV. Furthermore, we show that for 27 related binary and ternary semiconductors, B3PW91 predicts gaps with a MAD of only 0.09 eV, which is substantially better than all modern hybrid functionals, including B3LYP (MAD of 0.19 eV) and screened hybrid functional HSE06 (MAD of 0.18 eV).

The laboratory performance of CIGS solar cells (> 20% efficiency) makes them promising candidate photovoltaic devices. However, there remains little understanding of how defects at the CIGS/CdS interface affect the band offsets and interfacial energies, and hence the performance of manufactured devices. To determine these relationships, we use the B3PW91 hybrid functional of DFT with the AEP method that we validate to provide very accurate descriptions of both band gaps and band offsets. This confirms the weak dependence of band offsets on surface orientation observed experimentally. We predict that the CBO of perfect CuInSe₂/CdS interface is large, 0.79 eV, which would dramatically degrade performance. Moreover we show that band gap widening induced by Ga adjusts only the VBO, and we find that Cd impurities do not significantly affect the CBO. Thus we show that Cu vacancies at the interface play the key role in enabling the tunability of CBO. We predict that Na further improves the CBO through electrostatically elevating the valence levels to decrease the CBO, explaining the observed essential role of Na for high performance. Moreover we find that K leads to a dramatic decrease in the CBO to 0.05 eV, much better than Na. We suggest that the efficiency of CIGS devices might be improved substantially by tuning the ratio of Na to K, with the improved phase stability of Na balancing phase instability from K. All these defects reduce interfacial stability slightly, but not significantly.

A number of exotic structures have been formed through high pressure chemistry, but applications have been hindered by difficulties in recovering the high pressure phase to ambient conditions (i.e., one atmosphere and room temperature). Here we use dispersion-corrected DFT (PBE-ulg flavor) to predict that above 60 GPa the most stable form of N₂O (the laughing gas in its molecular form) is a 1D polymer with an all-nitrogen backbone analogous to cis-polyacetylene in which alternate N are bonded (ionic covalent) to O. The analogous trans-polymer is only 0.03-0.10 eV/molecular unit less stable. Upon relaxation to ambient conditions both polymers relax below 14 GPa to the same stable non-planar trans-polymer, accompanied by possible electronic structure transitions. The predicted phonon spectrum and dissociation kinetics validate the stability of this trans-poly-NNO at ambient conditions, which has potential applications as a new type of conducting polymer with all-nitrogen chains and as a high-energy oxidizer for rocket propulsion. This work illustrates in silico materials discovery particularly in the realm of extreme conditions.

Modeling non-adiabatic electron dynamics has been a long-standing challenge for computational chemistry and materials science, and the eFF method presents a cost-efficient alternative. However, due to the deficiency of FSG representation, eFF is limited to low-Z elements with electrons of predominant s-character. To overcome this, we introduce a formal set of ECP extensions that enable accurate description of p-block elements. The extensions consist of a model representing the core electrons with the nucleus as a single pseudo particle represented by FSG, interacting with valence electrons through ECPs. We demonstrate and validate the ECP extensions for complex bonding structures, geometries, and energetics of systems with p-block character (C, O, Al, Si) and apply them to study materials under extreme mechanical loading conditions.

Despite its success, the eFF framework has some limitations, originated from both the design of Pauli potentials and the FSG representation. To overcome these, we develop a new framework of two-level hierarchy that is a more rigorous and accurate successor to the eFF method. The fundamental level, GHA-QM, is based on a new set of Pauli potentials that renders exact QM level of accuracy for any FSG represented electron systems. To achieve this, we start with using exactly derived energy expressions for the same spin electron pair, and fitting a simple functional form, inspired by DFT, against open singlet electron pair curves (H₂ systems). Symmetric and asymmetric scaling factors are then introduced at this level to recover the QM total energies of multiple electron pair systems from the sum of local interactions. To complement the imperfect FSG representation, the AMPERE extension is implemented, and aims at embedding the interactions associated with both the cusp condition and explicit nodal structures. The whole GHA-QM+AMPERE framework is tested on H element, and the preliminary results are promising.

Molecular Dynamics Simulation on Characteristics of Decomposition of Methane Hydrate in the Presence of Hydrogen Peroxide Solution

In this work, the characteristics of the decomposition of methane hydrate Structure I (SI) in the presence of hydrogen peroxide solution is investigated using the molecular dynamics simulation. The mechanism of the transformation process from the solid hydrate to the liquid is analyzed with the effect of hydrogen peroxide (HP) solution. In addition, the effect of ethylene glycol (EG) with the same molar concentration with HP on the methane hydrate dissociation is also studied. The results illustrate that both HP and EG promote well the hydrate dissociation. The work provides the important reference value for the experimental investigation into the promotion effect of HP on the hydrate dissociation.