Molecular modelling of organic functional materials

The investigation of the mechanisms lying behind the (photo-)chemical processes is fundamental to address and improve the design of new organic functional materials. In many cases, dynamics simulations represent the only tool to capture the system properties emerging from complex interactions between many molecules. Despite the outstanding progresses in calculation power, the only way to carry out such computational studies is to introduce several approximations with respect to a fully quantum mechanical (QM) description. This thesis presents an approach that combines QM calculations with a classical Molecular Dynamics (MD) approach by means of accurate QM-derived force fields. It is based on a careful selection of the most relevant molecular degrees of freedom, whose potential energy surface is calculated at QM level and reproduced by the analytic functions of the force field, as well as by an accurate tuning of the approximations introduced in the model of the process to be simulated. This is made possible by some tools developed purposely, that allow to obtain and test the FF parameters through comparison with the QM frequencies and normal modes. These tools were applied in the modelling of three processes: the npi* photoisomerisation of azobenzene, where the FF description was extended to the excited state too and the non-adiabatic events were treated stochastically with Tully fewest switching algorithm; the charge separation in donors-acceptors bulk heterojunction organic solar cells, where a tight-binding Hamiltonian was carefully parametrised and solved by means of a code, also written specifically; the effect of the protonation state on the photoisomerisation quantum yield of the aryl-azoimidazolium unit of the axle molecule of a rotaxane molecular shuttle. In each case, the QM-based MD models that were specifically developed gave noteworthy information about the investigated phenomena, proving to be a fundamental key for a deeper comprehension of several experimental evidences.

Brownian motion of multidimensional systems in nonpotential velocity-dependent fields of force

We study the Brownian motion in velocity-dependent fields of force. Our main result is a Smoluchowski equation valid for moderate to high damping constants. We derive that equation by perturbative solution of the Langevin equation and using functional derivative techniques.

Investigation of Normal Force and Moment Coefficients for an AUV at Nonlinear Angle of Attack and Sideslip Range

This paper presents a comparative study of computational fluid dynamics (CFD) and analytical and semiempirical (ASE) methods applied to the prediction of the normal force and moment coefficients of an autonomous underwater vehicle (AUV). Both methods are applied to the. bare hull of the vehicle and to the body-hydroplane combination. The results are validated through experiments in a towing tank. It is shown that the CFD approach allows for a good prediction of the coefficients over the range of angles of attack considered. In contrast with the traditional ASE formulations used in naval and aircraft fields, an improved methodology is introduced that takes advantage of the qualitative information obtained from CFD flow visualizations.

The stochastic design of force-minimized compact magnets for high-field magnetic resonance imaging applications

New designs for force-minimized compact high-field clinical MRI magnets are described. The design method is a modified simulated annealing (SA) procedure which includes Maxwell forces in the error function to be minimized. This permits an automated force reduction in the magnet designs while controlling the overall dimensions of the system. As SA optimization requires many iterations to achieve a final design, it is important that each iteration in the procedure is rapid. We have therefore developed a rapid force calculation algorithm. Novel designs for short 3- and 4-T clinical MRI systems are presented in which force reduction has been invoked. The final designs provide large homogeneous regions and reduced stray fields in remarkable short magnets. A shielded 4-T design that is approximately 30% shorter than current designs is presented. This novel magnet generates a full 50-cm diameter homogeneous region.

Towards "green" smart materials for force and strain sensors: The case of polyaniline

Stress/strain sensors constitute a class of devices with a global ever-growing market thanks to their use in many fields of modern life. They are typically constituted by thin metal foils deposited on flexible supports. However, the low inherent resistivity and limited flexibility of their constituents make them inadequate for several applications, such as measuring large movements in robotic systems and biological tissues. As an alternative to the traditional compounds, in the present work we will show the advantages to employ a smart material, polyaniline (PANI), prepared by an innovative environmentally friendly route, for force/strain sensor applications wherein simple processing, environmental friendliness and sensitivity are particularly required.

Nanometer-scale oxidation of Si(100) surfaces by tapping mode atomic force microscopy

The nanometer¿scale oxidation of Si(100) surfaces in air is performed with an atomic force microscope working in tapping mode. Applying a positive voltage to the sample with respect to the tip, two kinds of modifications are induced on the sample: grown silicon oxide mounds less than 5 nm high and mounds higher than 10 nm (which are assumed to be gold depositions). The threshold voltage necessary to produce the modification is studied as a function of the average tip¿to¿sample distance.

Microcontroller-driven fluid-injection system for atomic force microscopy.

We present a programmable microcontroller-driven injection system for the exchange of imaging medium during atomic force microscopy. Using this low-noise system, high-resolution imaging can be performed during this process of injection without disturbance. This latter circumstance was exemplified by the online imaging of conformational changes in DNA molecules during the injection of anticancer drug into the fluid chamber.

Dynamic electrostatic force microscopy in liquid media

We present the implementation of dynamic electrostatic force microscopy in liquid media. This implementation enables the quantitative imaging of local dielectric properties of materials in electrolyte solutions with nanoscale spatial resolution. Local imaging capabilities are obtained by probing the frequency-dependent and ionic concentration-dependent electrostatic forces at high frequency (>1 MHz), while quantification of the interaction forces is obtained with finite-element numerical calculations. The results presented open a wide range of possibilities in a number of fields where the dielectric properties of materials need to be probed at the nanoscale and in a liquid environment.

Dynamic electrostatic force microscopy in liquid media

We present the implementation of dynamic electrostatic force microscopy in liquid media. This implementation enables the quantitative imaging of local dielectric properties of materials in electrolyte solutions with nanoscale spatial resolution. Local imaging capabilities are obtained by probing the frequency-dependent and ionic concentration-dependent electrostatic forces at high frequency (>1 MHz), while quantification of the interaction forces is obtained with finite-element numerical calculations. The results presented open a wide range of possibilities in a number of fields where the dielectric properties of materials need to be probed at the nanoscale and in a liquid environment.

Sketches of the Battle Fields of Ontario from the War of 1812, July 15, 1876.

Accompanying caption from the Canadian Illustrated News, July 15, 1876: “We publish today a page of sketches consisting of the following battle fields in Ontario :--Lundy’s Lane where, without doubt, the hardest fought battle of 1812-15 took place, and in which more troops were engaged than in any other engagement of that war : the battle field of Stony Creek where the Canadians and Indians made a night attack on the Americans and achieved a victory over a greatly superior force and obliged the Americans to retreat back to the shelter of Old Fort George which was the scene of many engagements during the war. Beaver Dam battle field is just in the suburbs of the thriving village of Thorold, and the monument covers the remains of several soldiers whose bodies were unearthed during the building of the new Welland Canal at that place.”

Studies on the dynamics and rheology of dilute suspensions of Brownian particles in simple shear flow under the action of an external force

We present a novel approach to computing the orientation moments and rheological properties of a dilute suspension of spheroids in a simple shear flow at arbitrary Peclct number based on a generalised Langevin equation method. This method differs from the diffusion equation method which is commonly used to model similar systems in that the actual equations of motion for the orientations of the individual particles are used in the computations, instead of a solution of the diffusion equation of the system. It also differs from the method of 'Brownian dynamics simulations' in that the equations used for the simulations are deterministic differential equations even in the presence of noise, and not stochastic differential equations as in Brownian dynamics simulations. One advantage of the present approach over the Fokker-Planck equation formalism is that it employs a common strategy that can be applied across a wide range of shear and diffusion parameters. Also, since deterministic differential equations are easier to simulate than stochastic differential equations, the Langevin equation method presented in this work is more efficient and less computationally intensive than Brownian dynamics simulations.We derive the Langevin equations governing the orientations of the particles in the suspension and evolve a procedure for obtaining the equation of motion for any orientation moment. A computational technique is described for simulating the orientation moments dynamically from a set of time-averaged Langevin equations, which can be used to obtain the moments when the governing equations are harder to solve analytically. The results obtained using this method are in good agreement with those available in the literature.The above computational method is also used to investigate the effect of rotational Brownian motion on the rheology of the suspension under the action of an external force field. The force field is assumed to be either constant or periodic. In the case of con- I stant external fields earlier results in the literature are reproduced, while for the case of periodic forcing certain parametric regimes corresponding to weak Brownian diffusion are identified where the rheological parameters evolve chaotically and settle onto a low dimensional attractor. The response of the system to variations in the magnitude and orientation of the force field and strength of diffusion is also analyzed through numerical experiments. It is also demonstrated that the aperiodic behaviour exhibited by the system could not have been picked up by the diffusion equation approach as presently used in the literature.The main contributions of this work include the preparation of the basic framework for applying the Langevin method to standard flow problems, quantification of rotary Brownian effects by using the new method, the paired-moment scheme for computing the moments and its use in solving an otherwise intractable problem especially in the limit of small Brownian motion where the problem becomes singular, and a demonstration of how systems governed by a Fokker-Planck equation can be explored for possible chaotic behaviour.

Protein Folding Simulations: Confinement, External Fields and Sequence Design

The present Thesis looks at the problem of protein folding using Monte Carlo and Langevin simulations, three topics in protein folding have been studied: 1) the effect of confining potential barriers, 2) the effect of a static external field and 3) the design of amino acid sequences which fold in a short time and which have a stable native state (global minimum). Regarding the first topic, we studied the confinement of a small protein of 16 amino acids known as 1NJ0 (PDB code) which has a beta-sheet structure as a native state. The confinement of proteins occurs frequently in the cell environment. Some molecules called Chaperones, present in the cytoplasm, capture the unfolded proteins in their interior and avoid the formation of aggregates and misfolded proteins. This mechanism of confinement mediated by Chaperones is not yet well understood. In the present work we considered two kinds of potential barriers which try to mimic the confinement induced by a Chaperon molecule. The first kind of potential was a purely repulsive barrier whose only effect is to create a cavity where the protein folds up correctly. The second kind of potential was a barrier which includes both attractive and repulsive effects. We performed Wang-Landau simulations to calculate the thermodynamical properties of 1NJ0. From the free energy landscape plot we found that 1NJ0 has two intermediate states in the bulk (without confinement) which are clearly separated from the native and the unfolded states. For the case of the purely repulsive barrier we found that the intermediate states get closer to each other in the free energy landscape plot and eventually they collapse into a single intermediate state. The unfolded state is more compact, compared to that in the bulk, as the size of the barrier decreases. For an attractive barrier modifications of the states (native, unfolded and intermediates) are observed depending on the degree of attraction between the protein and the walls of the barrier. The strength of the attraction is measured by the parameter $\epsilon$. A purely repulsive barrier is obtained for $\epsilon=0$ and a purely attractive barrier for $\epsilon=1$. The states are changed slightly for magnitudes of the attraction up to $\epsilon=0.4$. The disappearance of the intermediate states of 1NJ0 is already observed for $\epsilon =0.6$. A very high attractive barrier ($\epsilon \sim 1.0$) produces a completely denatured state. In the second topic of this Thesis we dealt with the interaction of a protein with an external electric field. We demonstrated by means of computer simulations, specifically by using the Wang-Landau algorithm, that the folded, unfolded, and intermediate states can be modified by means of a field. We have found that an external field can induce several modifications in the thermodynamics of these states: for relatively low magnitudes of the field ($<2.06 \times 10^8$ V/m) no major changes in the states are observed. However, for higher magnitudes than ($6.19 \times 10^8$ V/m) one observes the appearance of a new native state which exhibits a helix-like structure. In contrast, the original native state is a $\beta$-sheet structure. In the new native state all the dipoles in the backbone structure are aligned parallel to the field. The design of amino acid sequences constitutes the third topic of the present work. We have tested the Rate of Convergence criterion proposed by D. Gridnev and M. Garcia ({\it work unpublished}). We applied it to the study of off-lattice models. The Rate of Convergence criterion is used to decide if a certain sequence will fold up correctly within a relatively short time. Before the present work, the common way to decide if a certain sequence was a good/bad folder was by performing the whole dynamics until the sequence got its native state (if it existed), or by studying the curvature of the potential energy surface. There are some difficulties in the last two approaches. In the first approach, performing the complete dynamics for hundreds of sequences is a rather challenging task because of the CPU time needed. In the second approach, calculating the curvature of the potential energy surface is possible only for very smooth surfaces. The Rate of Convergence criterion seems to avoid the previous difficulties. With this criterion one does not need to perform the complete dynamics to find the good and bad sequences. Also, the criterion does not depend on the kind of force field used and therefore it can be used even for very rugged energy surfaces.

Effects of electric fields on ultracold Rydberg atom interactions

The behaviour of interacting ultracold Rydberg atoms in both constant electric fields and laser fields is important for designing experiments and constructing realistic models of them. In this paper, we briefly review our prior work and present new results on how electric fields affect interacting ultracold Rydberg atoms. Specifically, we address the topics of constant background electric fields on Rydberg atom pair excitation and laser-induced Stark shifts on pair excitation.

Gravitational Lorentz force and the description of the gravitational interaction

In the context of a gauge theory for the translation group, we have obtained, for a spinless particle, a gravitational analogue of the Lorentz force. Then, we have shown that this force equation can be rewritten in terms of magnitudes related to either the teleparallel or the Riemannian structures induced in spacetime by the presence of the gravitational field. In the first case, it gives a force equation, with torsion playing the role of force. In the second, it gives the usual geodesic equation of general relativity. The main conclusion is that scalar matter is able to feel any one of the above spacetime geometries, the teleparallel and the metric ones. Furthermore, both descriptions are found to be completely equivalent in the sense that they give the same physical trajectory for a spinless particle in a gravitational field.

On the nucleation of GaP/GaAs and the effect of buried stress fields

We have recently shown that spatial ordering for epitaxially grown InP dots can be obtained using the periodic stress field of compositional modulation on the InGaP buffer layer. The aim of this present work is to study the growth of films of GaP by Chemical Beam Epitaxy (CBE), with in-situ monitoring by Reflection High Energy Electron Diffraction (RHEED), on layers of unstressed and stressed GaAs. Complementary, we have studied the role of a buried InP dot array on GaP nucleation in order to obtain three-dimensional structures. In both cases, the topographical characteristics of the samples were investigated by Atomic Force Microscopy (AFM) in non-contact mode. Thus vertically-coupled quantum dots of different materials have been obtained keeping the in-place spatial ordering originated from the composition modulation. © 2006 Materials Research Society.