Multiscale Modeling and Simulation of Stepped Crystal Surfaces

A primary goal of this dissertation is to understand the links between mathematical models that describe crystal surfaces at three fundamental length scales: The scale of individual atoms, the scale of collections of atoms forming crystal defects, and macroscopic scale. Characterizing connections between different classes of models is a critical task for gaining insight into the physics they describe, a long-standing objective in applied analysis, and also highly relevant in engineering applications. The key concept I use in each problem addressed in this thesis is coarse graining, which is a strategy for connecting fine representations or models with coarser representations. Often this idea is invoked to reduce a large discrete system to an appropriate continuum description, e.g. individual particles are represented by a continuous density. While there is no general theory of coarse graining, one closely related mathematical approach is asymptotic analysis, i.e. the description of limiting behavior as some parameter becomes very large or very small. In the case of crystalline solids, it is natural to consider cases where the number of particles is large or where the lattice spacing is small. Limits such as these often make explicit the nature of links between models capturing different scales, and, once established, provide a means of improving our understanding, or the models themselves. Finding appropriate variables whose limits illustrate the important connections between models is no easy task, however. This is one area where computer simulation is extremely helpful, as it allows us to see the results of complex dynamics and gather clues regarding the roles of different physical quantities. On the other hand, connections between models enable the development of novel multiscale computational schemes, so understanding can assist computation and vice versa. Some of these ideas are demonstrated in this thesis. The important outcomes of this thesis include: (1) a systematic derivation of the step-flow model of Burton, Cabrera, and Frank, with corrections, from an atomistic solid-on-solid-type models in 1+1 dimensions; (2) the inclusion of an atomistically motivated transport mechanism in an island dynamics model allowing for a more detailed account of mound evolution; and (3) the development of a hybrid discrete-continuum scheme for simulating the relaxation of a faceted crystal mound. Central to all of these modeling and simulation efforts is the presence of steps composed of individual layers of atoms on vicinal crystal surfaces. Consequently, a recurring theme in this research is the observation that mesoscale defects play a crucial role in crystal morphological evolution.

MINING AND VERIFICATION OF TEMPORAL EVENTS WITH APPLICATIONS IN COMPUTER MICRO-ARCHITECTURE RESEARCH

Computer simulation programs are essential tools for scientists and engineers to understand a particular system of interest. As expected, the complexity of the software increases with the depth of the model used. In addition to the exigent demands of software engineering, verification of simulation programs is especially challenging because the models represented are complex and ridden with unknowns that will be discovered by developers in an iterative process. To manage such complexity, advanced verification techniques for continually matching the intended model to the implemented model are necessary. Therefore, the main goal of this research work is to design a useful verification and validation framework that is able to identify model representation errors and is applicable to generic simulators. The framework that was developed and implemented consists of two parts. The first part is First-Order Logic Constraint Specification Language (FOLCSL) that enables users to specify the invariants of a model under consideration. From the first-order logic specification, the FOLCSL translator automatically synthesizes a verification program that reads the event trace generated by a simulator and signals whether all invariants are respected. The second part consists of mining the temporal flow of events using a newly developed representation called State Flow Temporal Analysis Graph (SFTAG). While the first part seeks an assurance of implementation correctness by checking that the model invariants hold, the second part derives an extended model of the implementation and hence enables a deeper understanding of what was implemented. The main application studied in this work is the validation of the timing behavior of micro-architecture simulators. The study includes SFTAGs generated for a wide set of benchmark programs and their analysis using several artificial intelligence algorithms. This work improves the computer architecture research and verification processes as shown by the case studies and experiments that have been conducted.

Simulation of photoinduced processes in organic chromophores: mapping photochemically relevant states for accurate dynamics

The simulation of ultrafast photoinduced processes is a fundamental step towards the understanding of the underlying molecular mechanism and interpretation/prediction of experimental data. Performing a computer simulation of a complex photoinduced process is only possible introducing some approximations but, in order to obtain reliable results, the need to reduce the complexity must balance with the accuracy of the model, which should include all the relevant degrees of freedom and a quantitatively correct description of the electronic states involved in the process. This work presents new computational protocols and strategies for the parameterisation of accurate models for photochemical/photophysical processes based on state-of-the-art multiconfigurational wavefunction-based methods. The required ingredients for a dynamics simulation include potential energy surfaces (PESs) as well as electronic state couplings, which must be mapped across the wide range of geometries visited during the wavepacket/trajectory propagation. The developed procedures allow to obtain solid and extended databases reducing as much as possible the computational cost, thanks to, e.g., specific tuning of the level of theory for different PES regions and/or direct calculation of only the needed components of vectorial quantities (like gradients or nonadiabatic couplings). The presented approaches were applied to three case studies (azobenzene, pyrene, visual rhodopsin), all requiring an accurate parameterisation but for different reasons. The resulting models and simulations allowed to elucidate the mechanism and time scale of the internal conversion, reproducing or even predicting new transient experiments. The general applicability of the developed protocols to systems with different peculiarities and the possibility to parameterise different types of dynamics on an equal footing (classical vs purely quantum) prove that the developed procedures are flexible enough to be tailored for each specific system, and pave the way for exact quantum dynamics with multiple degrees of freedom.

Photoelastic And Finite Element Analyses Of Occlusal Loads In Mandibular Body.

This study proposed to evaluate the mandibular biomechanics in the posterior dentition based on experimental and computational analyses. The analyses were performed on a model of human mandible, which was modeled by epoxy resin for photoelastic analysis and by computer-aided design for finite element analysis. To standardize the evaluation, specific areas were determined at the lateral surface of mandibular body. The photoelastic analysis was configured through a vertical load on the first upper molar and fixed support at the ramus of mandible. The same configuration was used in the computer simulation. Force magnitudes of 50, 100, 150, and 200 N were applied to evaluate the bone stress. The stress results presented similar distribution in both analyses, with the more intense stress being at retromolar area and oblique line and alveolar process at molar level. This study presented the similarity of results in the experimental and computational analyses and, thus, showed the high importance of morphology biomechanical characterization at posterior dentition.

Resultados preliminares de um algoritmo para ablação de lente de contato personalizada

OBJETIVO: Desenvolver simulação computadorizada de ablação para produzir lentes de contato personalizadas a fim de corrigir aberrações de alta ordem. MÉTODOS: Usando dados reais de um paciente com ceratocone, mensurados em um aberrômetro ("wavefront") com sensor Hartmann-Shack, foram determinados as espessuras de lentes de contato que compensam essas aberrações assim como os números de pulsos necessários para fazer ablação as lentes especificamente para este paciente. RESULTADOS: Os mapas de correção são apresentados e os números dos pulsos foram calculados, usando feixes com a largura de 0,5 mm e profundidade de ablação de 0,3 µm. CONCLUSÕES: Os resultados simulados foram promissores, mas ainda precisam ser aprimorados para que o sistema de ablação "real" possa alcançar a precisão desejada.

AcMus: an open, integrated platform for room acoustics research

This article describes the design, implementation, and experiences with AcMus, an open and integrated software platform for room acoustics research, which comprises tools for measurement, analysis, and simulation of rooms for music listening and production. Through use of affordable hardware, such as laptops, consumer audio interfaces and microphones, the software allows evaluation of relevant acoustical parameters with stable and consistent results, thus providing valuable information in the diagnosis of acoustical problems, as well as the possibility of simulating modifications in the room through analytical models. The system is open-source and based on a flexible and extensible Java plug-in framework, allowing for cross-platform portability, accessibility and experimentation, thus fostering collaboration of users, developers and researchers in the field of room acoustics.

Uma perspectiva computacional sobre catálise enzimática

Enzymes are extremely efficient catalysts. Here, part of the mechanisms proposed to explain this catalytic power will be compared to quantitative experimental results and computer simulations. Influence of the enzymatic environment over species along the reaction coordinate will be analysed. Concepts of transition state stabilisation and reactant destabilisation will be confronted. Divided site model and near-attack conformation hypotheses will also be discussed. Molecular interactions such as covalent catalysis, general acid-base catalysis, electrostatics, entropic effects, steric hindrance, quantum and dynamical effects will also be analysed as sources of catalysis. Reaction mechanisms, in particular that catalysed by protein tyrosine phosphatases, illustrate the concepts.

Surface plasmon resonance of gold nanoparticles formed by cathodic arc plasma ion implantation into polymer

Shallow subsurface layers of gold nanoclusters were formed in polymethylmethacrylate (PMMA) polymer by very low energy (49 eV) gold ion implantation. The ion implantation process was modeled by computer simulation and accurately predicted the layer depth and width. Transmission electron microscopy (TEM) was used to image the buried layer and individual nanoclusters; the layer width was similar to 6-8 nm and the cluster diameter was similar to 5-6 nm. Surface plasmon resonance (SPR) absorption effects were observed by UV-visible spectroscopy. The TEM and SPR results were related to prior measurements of electrical conductivity of Au-doped PMMA, and excellent consistency was found with a model of electrical conductivity in which either at low implantation dose the individual nanoclusters are separated and do not physically touch each other, or at higher implantation dose the nanoclusters touch each other to form a random resistor network (percolation model). (C) 2009 American Vacuum Society. [DOI: 10.1116/1.3231449]

Gold-implanted shallow conducting layers in polymethylmethacrylate

PMMA (polymethylmethacrylate) was ion implanted with gold at very low energy and over a range of different doses using a filtered cathodic arc metal plasma system. A nanometer scale conducting layer was formed, fully buried below the polymer surface at low implantation dose, and evolving to include a gold surface layer as the dose was increased. Depth profiles of the implanted material were calculated using the Dynamic TRIM computer simulation program. The electrical conductivity of the gold-implanted PMMA was measured in situ as a function of dose. Samples formed at a number of different doses were subsequently characterized by Rutherford backscattering spectrometry, and test patterns were formed on the polymer by electron beam lithography. Lithographic patterns were imaged by atomic force microscopy and demonstrated that the contrast properties of the lithography were well maintained in the surface-modified PMMA.

A microscopic view of substitution reactions solvated by ionic liquids

The solvation effect of the ionic liquid 1-N-butyl-3-methylimidazolium hexafluorophosphate on nucleophilic substitution reactions of halides toward the aliphatic carbon of methyl p-nitrobenzenesulfonate (pNBS) was investigated by computer simulations. The calculations were performed by using a hybrid quantum-mechanical/molecular-mechanical (QM/MM) methodology. A semiempirical Hamiltonian was first parametrized on the basis of comparison with ab initio calculations for Cl(-) and Br(-) reaction with pNBS at gas phase. In condensed phase, free energy profiles were obtained for both reactions. The calculated reaction barriers are in agreement with experiment. The structure of species solvated by the ionic liquid was followed along the reaction progress from the reagents, through the transition state, to the final products. The simulations indicate that this substitution reaction in the ionic liquid is slower than in nonpolar molecular solvents proper to significant stabilization of the halide anion by the ionic liquid in comparison with the transition state with delocalized charge. Solute-solvent interactions in the first solvation shell contain several hydrogen bonds that are formed or broken in response to charge density variation along the reaction coordinate. The detailed structural analysis can be used to rationalize the design of new ionic liquids with tailored solvation properties. (c) 2008 American Institute of Physics.

Low-frequency Raman spectra and fragility of imidazolium ionic liquids

Raman spectra within the 5-200 cm(-1) range have been recorded as a function of temperature for different ionic liquids based on imidazolium cations. A correlation has been found between fragility and the temperature dependence of the strength of fast relaxational motions. Understanding quasielastic scattering as the relaxational contribution to ionic mean-squared displacement elucidates some effects on ionic liquids' fragility resulting from modifications in the chemical structure. (C) 2010 American Institute of Physics. [doi:10.1063/1.3462962]

Polarization effects in molecular dynamics simulations of glass-formers Ca(NO(3))(2)center dot nH(2)O, n=4, 6, and 8

Thermodynamics, equilibrium structure, and dynamics of glass-forming liquids Ca(NO(3))(2)center dot nH(2)O, n=4, 6, and 8, have been investigated by molecular dynamics (MD) simulations. A polarizable model was considered for H(2)O and NO(3)- on the basis of previous fluctuating charge models for pure water and the molten salt 2Ca(NO(3))(2)center dot 3KNO(3). Similar thermodynamic properties have been obtained with nonpolarizable and polarizable models. The glass transition temperature, T(g), estimated from MD simulations was dependent on polarization, in particular the dependence of T(g) with electrolyte concentration. Significant polarization effects on equilibrium structure were observed in cation-cation, cation-anion, and water-water structures. Polarization increases the diffusion coefficient of H(2)O, but does not change significantly the diffusion coefficients of ions. Viscosity decreases upon inclusion of polarization, but the conductivity calculated with the polarizable model is smaller than the nonpolarizable model because polarization enhances anion-cation interactions.

Intermolecular vibrations and fast relaxations in supercooled ionic liquids

Short-time dynamics of ionic liquids has been investigated by low-frequency Raman spectroscopy (4 < omega < 100 cm(-1)) within the supercooled liquid range. Raman spectra are reported for ionic liquids with the same anion, bis(trifluoromethylsulfonyl)imide, and different cations: 1-butyl-3-methylimidazolium, 1-hexyl-3-methylimidazolium, 1-butyl-1-methylpiperidinium, trimethylbutylammonium, and tributylmethylammonium. It is shown that low-frequency Raman spectroscopy provides similar results as optical Kerr effect (OKE) spectroscopy, which has been used to study intermolecular vibrations in ionic liquids. The comparison of ionic liquids containing aromatic and non-aromatic cations identifies the characteristic feature in Raman spectra usually assigned to librational motion of the imidazolium ring. The strength of the fast relaxations (quasi-elastic scattering, QES) and the intermolecular vibrational contribution (boson peak) of ionic liquids with non-aromatic cations are significantly lower than imidazolium ionic liquids. A correlation length assigned to the boson peak vibrations was estimated from the frequency of the maximum of the boson peak and experimental data of sound velocity. The correlation length related to the boson peak (similar to 19 angstrom) does not change with the length of the alkyl chain in imidazolium cations, in contrast to the position of the first-sharp diffraction peak observed in neutron and X-ray scattering measurements of ionic liquids. The rate of change of the QES intensity in the supercooled liquid range is compared with data of excess entropy, free volume, and mean-squared displacement recently reported for ionic liquids. The temperature dependence of the QES intensity in ionic liquids illustrates relationships between short-time dynamics and long-time structural relaxation that have been proposed for glass-forming liquids. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3604533]

Building a Open Source Framework for Virtual Medical Training

This paper presents a framework to build medical training applications by using virtual reality and a tool that helps the class instantiation of this framework. The main purpose is to make easier the building of virtual reality applications in the medical training area, considering systems to simulate biopsy exams and make available deformation, collision detection, and stereoscopy functionalities. The instantiation of the classes allows quick implementation of the tools for such a purpose, thus reducing errors and offering low cost due to the use of open source tools. Using the instantiation tool, the process of building applications is fast and easy. Therefore, computer programmers can obtain an initial application and adapt it to their needs. This tool allows the user to include, delete, and edit parameters in the functionalities chosen as well as storing these parameters for future use. In order to verify the efficiency of the framework, some case studies are presented.

Thermally activated peroxydisulfate in the presence of additives: A clean method for the degradation of pollutants

The kinetics and mechanism of the thermal activation of peroxydisulfate, in the temperature range from 60 to 80 degrees C, was investigated in the presence and absence of sodium formate as an additive to turn the oxidizing capacity of the reaction mixture into a reductive one. Trichloroacetic acid, TCA, whose degradation by a reductive mechanism is well reported in the literature, was used as a probe. The chemistry of thermally activated peroxydisulfate is described by a reaction scheme involving free radical generation. The proposed mechanism is evaluated by a computer simulation of the concentration profiles obtained under different experimental conditions. In the presence of formate, SO(4)(center dot-) radicals yield CO(2)(center dot-), which are the main species available for degrading TCA. Under the latter conditions, TCA is more efficiently depleted than in the absence of formate, but otherwise identical conditions of temperature and [S(2)O(8)(2-)]. We therefore conclude that activated peroxydisulfate in the presence of formate as an additive is a convenient method for the mineralization of substrates that are refractory to oxidation. such as perchlorinated hydrocarbons and TCA. This method has the advantage that leaves no toxic residues. (C) 2009 Elsevier Ltd. All rights reserved.