Universal tight binding model for chemical reactions in solution and at surfaces. III. Stoichiometric and reduced surfaces of titania and the adsorption of water

We demonstrate a model for stoichiometric and reduced titanium dioxide intended for use in molecular dynamics and other atomistic simulations and based in the polarizable ion tight binding theory. This extends the model introduced in two previous papers from molecular and liquid applications into the solid state, thus completing the task of providing a comprehensive and unified scheme for studying chemical reactions, particularly aimed at problems in catalysis and electrochemistry. As before, experimental results are given priority over theoretical ones in selecting targets for model fitting, for which we used crystal parameters and band gaps of titania bulk polymorphs, rutile and anatase. The model is applied to six low index titania surfaces, with and without oxygen vacancies and adsorbed water molecules, both in dissociated and non-dissociated states. Finally, we present the results of molecular dynamics simulation of an anatase cluster with a number of adsorbed water molecules and discuss the role of edge and corner atoms of the cluster. (C) 2014 AIP Publishing LLC.

Cyclic nanoindentation and nano-impact fatigue characterisation of functionally graded TiN/TiNi shape memory film

This paper investigates the mechanism of nanoscale fatigue of functionally graded TiN/TiNi films using nano-impact and multiple-loading-cycle nanoindentation tests. The functionally graded films were deposited on silicon substrate, in which TiNi films maintain shape memory and pseudo elastic behavior, while a modified TiN surface layer provides tribological and anti-corrosion properties. Nanomechanical tests were performed to comprehend the localized film performance and failure modes of the functionally graded film using NanoTestTM equipped with Berkovich and conical indenter between 100 μN to 500 mN loads. The loading mechanism and load history are critical to define film failure modes (i.e. backward depth deviation) including the shape memory effect of the functionally graded layer. The results are sensitive to the applied load, loading type (e.g. semi-static, dynamic) and probe geometry. Based on indentation force-depth profiles, depth-time data and post-test surface observations of films, it is concluded that the shape of the nanoindenter is critical in inducing the localized indentation stress and film failure, including shape recovery at the lower load range. Elastic-plastic finite element (FE) simulation during nanoindentation loading indicated that the location of subsurface maximum stress near the interface influences the backward depth deviation type of film failure. A standalone, molecular dynamics simulation was performed with the help of a long range potential energy function to simulate the tensile test of TiN nanowire with two different aspect ratios to investigate the theory of its failure mechanism.

Cloud benchmarking for performance

How can applications be deployed on the cloud to achieve maximum performance? This question has become significant and challenging with the availability of a wide variety of Virtual Machines (VMs) with different performance capabilities in the cloud. The above question is addressed by proposing a six step benchmarking methodology in which a user provides a set of four weights that indicate how important each of the following groups: memory, processor, computation and storage are to the application that needs to be executed on the cloud. The weights along with cloud benchmarking data are used to generate a ranking of VMs that can maximise performance of the application. The rankings are validated through an empirical analysis using two case study applications, the first is a financial risk application and the second is a molecular dynamics simulation, which are both representative of workloads that can benefit from execution on the cloud. Both case studies validate the feasibility of the methodology and highlight that maximum performance can be achieved on the cloud by selecting the top ranked VMs produced by the methodology.

Linear Optics Simulation of Quantum Non-Markovian Dynamics

The simulation of open quantum dynamics has recently allowed the direct investigation of the features of system-environment interaction and of their consequences on the evolution of a quantum system. Such interaction threatens the quantum properties of the system, spoiling them and causing the phenomenon of decoherence. Sometimes however a coherent exchange of information takes place between system and environment, memory effects arise and the dynamics of the system becomes non-Markovian. Here we report the experimental realisation of a non-Markovian process where system and environment are coupled through a simulated transverse Ising model. By engineering the evolution in a photonic quantum simulator, we demonstrate the role played by system-environment correlations in the emergence of memory effects.

Energy conserving schemes for the simulation of musical instrument contact dynamics

Collisions are an innate part of the function of many musical instruments. Due to the nonlinear nature of contact forces, special care has to be taken in the construction of numerical schemes for simulation and sound synthesis. Finite difference schemes and other time-stepping algorithms used for musical instrument modelling purposes are normally arrived at by discretising a Newtonian description of the system. However because impact forces are non-analytic functions of the phase space variables, algorithm stability can rarely be established this way. This paper presents a systematic approach to deriving energy conserving schemes for frictionless impact modelling. The proposed numerical formulations follow from discretising Hamilton׳s equations of motion, generally leading to an implicit system of nonlinear equations that can be solved with Newton׳s method. The approach is first outlined for point mass collisions and then extended to distributed settings, such as vibrating strings and beams colliding with rigid obstacles. Stability and other relevant properties of the proposed approach are discussed and further demonstrated with simulation examples. The methodology is exemplified through a case study on tanpura string vibration, with the results confirming the main findings of previous studies on the role of the bridge in sound generation with this type of string instrument.

Is Northern Ireland abnormal? An extension of the sociological debate on religion in modern Britain

This article places Northern Ireland within the unfolding sociological debate on religion in modern Britain. It measures secularization along Casanova’s three dimensions (1994): religious differentiation, decline and privatization. It finds that Northern Ireland has, in common with Britain, high levels of religious differentiation, grey areas of religious belief and little convinced secularism. However, Northern Ireland differs in that it has higher levels of religious affiliation and practice, and religion plays more roles in civil society than it does in other parts of Britain. The article explores the role of conflict in forming these religious trends, asking if they represent a persistence of the sacred, or simply mask deeper ethnic divisions. It concludes that the social dimensions of religion are just as important as the supernatural, and that they often inform each other. Finally, it suggests that the dynamics of religious change are comparable across regions and, as such, Northern Ireland might be a useful case study for British policy makers, particularly as it becomes increasingly multicultural and religiously plural.

Kinetics of phase transformations in a model with metastable fluid-fluid separation: A molecular dynamics study

By molecular dynamics (MD) simulations we study the crystallization process in a model system whose particles interact by a spherical pair potential with a narrow and deep attractive well adjacent to a hard repulsive core. The phase diagram of the model displays a solid-fluid equilibrium, with a metastable fluid-fluid separation. Our computations are restricted to fairly small systems (from 2592 to 10368 particles) and cover long simulation times, with constant energy trajectories extending up to 76x10(6) MD steps. By progressively reducing the system temperature below the solid-fluid line, we first observe the metastable fluid-fluid separation, occurring readily and almost reversibly upon crossing the corresponding line in the phase diagram. The nucleation of the crystal phase takes place when the system is in the two-fluid metastable region. Analysis of the temperature dependence of the nucleation time allows us to estimate directly the nucleation free energy barrier. The results are compared with the predictions of classical nucleation theory. The critical nucleus is identified, and its structure is found to be predominantly fcc. Following nucleation, the solid phase grows steadily across the system, incorporating a large number of localized and extended defects. We discuss the relaxation processes taking place both during and after the crystallization stage. The relevance of our simulation for the kinetics of protein crystallization under normal experimental conditions is discussed. (C) 2002 American Institute of Physics.

Simulation of the surface structure of butylmethylimidazolium ionic liquids

Molecular dynamics simulations of the liquid/vacuum surfaces of the room temperature ionic liquids [bmim][PF6], [bmim][BF4] and [bmim][Cl] have been carried out at various temperatures. The surfaces are structured with a top monolayer containing oriented cations and anions. The butyl side chains tend to face the vacuum and the methyl side chains the liquid. However, as the butyl chains are not densely packed, both anions and rings are visible from the vacuum phase. The effects of temperature and the anion on the degree of cation orientation is small, but the potential drop from the vacuum to the interior of the liquid is greater for liquids with smaller anions. We compare the simulation results with a range of experimental observations and suggest that neutron reflection from samples with protiated butyl groups would be a sensitive probe of the structure.

Simulation of interfaces between room temperature ionic liquids and other liquids

The structure and properties of the interfaces between the room temperature ionic liquid dimethylimidazolium chloride ([dmim]Cl) and different Lennard-Jones fluids and between ionic liquid and water have been studied by molecular dynamics simulations, and compared to the ionic liquid-vapour interface. Two contrasting types of interface were investigated, thermodynamically stable interfaces between ionic liquid and vapour and between ionic liquid and Lennard-Jones fluids, and diffusing interfaces between miscible phases of different compositions involving water. The density profiles of different species through the interface are presented. The cations and water molecules near the former type of interface are aligned relative to the surface, but no orientational preference was found near or in the broad diffusing interface. The ionic liquid has a negative electrostatic potential relative to vapour or Lennard-Jones fluid, but is more positive than pure water. This contrast is explained in terms of the relative importance of orientation and concentration differences in the two types of interface.

Observing time-dependent vibrational quantum dynamics in deuterium hydride molecular ions

Few-cycle laser pulses are used to "pump and probe" image the vibrational wavepacket dynamics of a HD+ molecular ion. The quantum dephasing and revival structure of the wavepacket are mapped experimentally with time-resolved photodissociation imaging. The motion of the molecule is simulated using a quantum-mechanical model predicting the observed structure. The coherence of the wavepacket is controlled by varying the duration of the intense laser pulses. By means of a Fourier transform analysis both the periodicity and relative population of the vibrational states of the excited molecular ion have been characterized.

Dynamical simulation of inelastic quantum transport

A method for correlated quantum electron-ion dynamics is combined with a method for electronic open boundaries to simulate in real time the heating, and eventual equilibration at an elevated vibrational energy, of a quantum ion under current flow in an atomic wire, together with the response of the current to the ionic heating. The method can also be used to extract inelastic current voltage corrections under steady-state conditions. However, in its present form the open-boundary method contains an approximation that limits the resolution of current-voltage features. The results of the simulations are tested against analytical results from scattering theory. Directions for the improvement of the method are summarized at the end.

A Van der Pol-Mathieu equation for the dynamics of dust grain charge in dusty plasmas

The chaotic profile of dust grain dynamics associated with dust-acoustic oscillations in a dusty plasma is considered. The collective behaviour of the dust plasma component is described via a multi-fluid model, comprising Boltzmann distributed electrons and ions, as well as an equation of continuity possessing a source term for the dust grains, the dust momentum and Poisson's equations. A Van der Pol–Mathieu-type nonlinear ordinary differential equation for the dust grain density dynamics is derived. The dynamical system is cast into an autonomous form by employing an averaging method. Critical stability boundaries for a particular trivial solution of the governing equation with varying parameters are specified. The equation is analysed to determine the resonance region, and finally numerically solved by using a fourth-order Runge–Kutta method. The presence of chaotic limit cycles is pointed out.

Postglacial formation and dynamics of North Patagonian Rainforest in the Chonos Archipelago, Southern Chile

Pollen analysis of continuous sediment cores from two lakes in the northern Chonos Archipelago (44S) in southern Chile shows a complete postglacial record of vegetation change. The fossil records indicate that deglaciation was complete in the northern Chonos by at least 13,600 14Cyr BP. Ericaceous heath and grassland persisted for more than 600 years after deglaciation under the influence of dry/cold climates and frequent burning. Nothofagus-Pilgerodendron-Podocarpus forest, with modern analogues in the southern Chonos Archipelago, was established across the northern islands by 12,400 14Cyr BP under increasingly warm and wet climates. There is no evidence for a return to cooler climates during the Younger Dryas chronozone. The rise of Tepualia stipularis and Weinmannia trichosperma as important forest components between 10,600 and 6000 14Cyr BP may be associated withclimates that were warmer than present. The collapse of Pilgerodendron communities during this time may have been triggered by a combination of factors related to disturbance frequency including tephra deposition events, fire and climate change. After 6000 14Cyr BP Pilgerodendron recovers and Nothofagus-Pilgerodendron-Tepualia forest persists until the present. European logging and burning activity may have increased the susceptibility of North Patagonian Rainforest to invasion by introduced species and to future collapse of the long-lived Pilgerodendron communities.