Object Kinetic Monte Carlo calculations of irradiated Fe-Cr dilute alloys: The effect of the interaction radius between substitutional Cr and self-interstitial Fe

ObjectKineticMonteCarlo models allow for the study of the evolution of the damage created by irradiation to time scales that are comparable to those achieved experimentally. Therefore, the essential ObjectKineticMonteCarlo parameters can be validated through comparison with experiments. However, this validation is not trivial since a large number of parameters is necessary, including migration energies of point defects and their clusters, binding energies of point defects in clusters, as well as the interactionradii. This is particularly cumbersome when describing an alloy, such as the Fe–Cr system, which is of interest for fusion energy applications. In this work we describe an ObjectKineticMonteCarlo model for Fe–Cr alloys in the dilute limit. The parameters used in the model come either from density functional theory calculations or from empirical interatomic potentials. This model is used to reproduce isochronal resistivity recovery experiments of electron irradiateddiluteFe–Cr alloys performed by Abe and Kuramoto. The comparison between the calculated results and the experiments reveal that an important parameter is the capture radius between substitutionalCr and self-interstitialFe atoms. A parametric study is presented on the effect of the capture radius on the simulated recovery curves.

Work Fluctuations in Bosonic Josephson Junctions

We calculate the first two moments and full probability distribution of the work performed on a system of bosonic particles in a two-mode Bose-Hubbard Hamiltonian when the self-interaction term is varied instantaneously or with a finite-time ramp. In the instantaneous case, we show how the irreversible work scales differently depending on whether the system is driven to the Josephson or Fock regime of the bosonic Josephson junction. In the finite-time case, we use optimal control techniques to substantially decrease the irreversible work to negligible values. Our analysis can be implemented in present-day experiments with ultracold atoms and we show how to relate the work statistics to that of the population imbalance of the two modes.

Theoretical investigation of thermal tweezers for parallel manipulation of atoms and nanoparticles on surfaces

A major focus of research in nanotechnology is the development of novel, high throughput techniques for fabrication of arbitrarily shaped surface nanostructures of sub 100 nm to atomic scale. A related pursuit is the development of simple and efficient means for parallel manipulation and redistribution of adsorbed atoms, molecules and nanoparticles on surfaces – adparticle manipulation. These techniques will be used for the manufacture of nanoscale surface supported functional devices in nanotechnologies such as quantum computing, molecular electronics and lab-on-achip, as well as for modifying surfaces to obtain novel optical, electronic, chemical, or mechanical properties. A favourable approach to formation of surface nanostructures is self-assembly. In self-assembly, nanostructures are grown by aggregation of individual adparticles that diffuse by thermally activated processes on the surface. The passive nature of this process means it is generally not suited to formation of arbitrarily shaped structures. The self-assembly of nanostructures at arbitrary positions has been demonstrated, though these have typically required a pre-patterning treatment of the surface using sophisticated techniques such as electron beam lithography. On the other hand, a parallel adparticle manipulation technique would be suited for directing the selfassembly process to occur at arbitrary positions, without the need for pre-patterning the surface. There is at present a lack of techniques for parallel manipulation and redistribution of adparticles to arbitrary positions on the surface. This is an issue that needs to be addressed since these techniques can play an important role in nanotechnology. In this thesis, we propose such a technique – thermal tweezers. In thermal tweezers, adparticles are redistributed by localised heating of the surface. This locally enhances surface diffusion of adparticles so that they rapidly diffuse away from the heated regions. Using this technique, the redistribution of adparticles to form a desired pattern is achieved by heating the surface at specific regions. In this project, we have focussed on the holographic implementation of this approach, where the surface is heated by holographic patterns of interfering pulsed laser beams. This implementation is suitable for the formation of arbitrarily shaped structures; the only condition is that the shape can be produced by holographic means. In the simplest case, the laser pulses are linearly polarised and intersect to form an interference pattern that is a modulation of intensity along a single direction. Strong optical absorption at the intensity maxima of the interference pattern results in approximately a sinusoidal variation of the surface temperature along one direction. The main aim of this research project is to investigate the feasibility of the holographic implementation of thermal tweezers as an adparticle manipulation technique. Firstly, we investigate theoretically the surface diffusion of adparticles in the presence of sinusoidal modulation of the surface temperature. Very strong redistribution of adparticles is predicted when there is strong interaction between the adparticle and the surface, and the amplitude of the temperature modulation is ~100 K. We have proposed a thin metallic film deposited on a glass substrate heated by interfering laser beams (optical wavelengths) as a means of generating very large amplitude of surface temperature modulation. Indeed, we predict theoretically by numerical solution of the thermal conduction equation that amplitude of the temperature modulation on the metallic film can be much greater than 100 K when heated by nanosecond pulses with an energy ~1 mJ. The formation of surface nanostructures of less than 100 nm in width is predicted at optical wavelengths in this implementation of thermal tweezers. Furthermore, we propose a simple extension to this technique where spatial phase shift of the temperature modulation effectively doubles or triples the resolution. At the same time, increased resolution is predicted by reducing the wavelength of the laser pulses. In addition, we present two distinctly different, computationally efficient numerical approaches for theoretical investigation of surface diffusion of interacting adparticles – the Monte Carlo Interaction Method (MCIM) and the random potential well method (RPWM). Using each of these approaches we have investigated thermal tweezers for redistribution of both strongly and weakly interacting adparticles. We have predicted that strong interactions between adparticles can increase the effectiveness of thermal tweezers, by demonstrating practically complete adparticle redistribution into the low temperature regions of the surface. This is promising from the point of view of thermal tweezers applied to directed self-assembly of nanostructures. Finally, we present a new and more efficient numerical approach to theoretical investigation of thermal tweezers of non-interacting adparticles. In this approach, the local diffusion coefficient is determined from solution of the Fokker-Planck equation. The diffusion equation is then solved numerically using the finite volume method (FVM) to directly obtain the probability density of adparticle position. We compare predictions of this approach to those of the Ermak algorithm solution of the Langevin equation, and relatively good agreement is shown at intermediate and high friction. In the low friction regime, we predict and investigate the phenomenon of ‘optimal’ friction and describe its occurrence due to very long jumps of adparticles as they diffuse from the hot regions of the surface. Future research directions, both theoretical and experimental are also discussed.

Porphyrins with metal, metalloid or phosphorus atoms directly bonded to the carbon periphery

Organometallic porphyrins with a metal, metalloid or phosphorus fragment directly attached to their carbon framework emerged for the first time in 1976, and these macrocycles have been intensively investigated in the past decade. The present review summarises for the first time all reported examples as well as applications of these systems.

Kinetic energy of He atoms in liquid He-4-He-3 mixtures

Deep inelastic neutron scattering measurements on liquid 3He-4He mixtures in the normal phase have been performed on the VESUVIO spectrometer at the ISIS pulsed neutron source at exchanged wave vectors of about q≃120.0Å-1. The neutron Compton profiles J(y) of the mixtures were measured along the T=1.96K isotherm for 3He concentrations, x, ranging from 0.1 to 1.0 at saturated vapor pressures. Values of kinetic energies 〈T〉 of 3He and 4He atoms as a function of x, 〈T〉(x), were extracted from the second moment of J(y). The present determinations of 〈T〉(x) confirm previous experimental findings for both isotopes and, in the case of 3He, a substantial disagreement with theory is found. In particular 〈T〉(x) for the 3He atoms is found to be independent of concentration yielding a value 〈T〉3(x=0.1)≃12K, much lower than the value suggested by the most recent theoretical estimates of approximately 19 K.

Atomic hydrogen diffusion in novel magnesium nanostructures : the impact of incorporated subsurface carbon atoms

Ab initio Density Functional Theory (DFT) calculations are performed to study the diffusion of atomic hydrogen on a Mg(0001) surface and their migration into the subsurface layers. A carbon atom located initially on a Mg(0001) surface can migrate into the sub-surface layer and occupy a fcc site, with charge transfer to the C atom from neighboring Mg atoms. The cluster of postively charged Mg atoms surrounding a sub-surface C is then shown to facilitate the dissociative chemisorption of molecular hydrogen on the Mg(0001) surface, and the surface migration and subsequent diffusion into the subsurface of atomic hydrogen. This helps rationalize the experimentally-observed improvement in absorption kinetics of H2 when graphite or single walled carbon nanotubes (SWCNT) are introduced into the Mg powder during ball milling.

Nuclear instruments and methods in physics research section B : beam interactions with materials and atoms

The impact-induced deposition of Al13 clusters with icosahedral structure on Ni(0 0 1) surface was studied by molecular dynamics (MD) simulation using Finnis–Sinclair potentials. The incident kinetic energy (Ein) ranged from 0.01 to 30 eV per atom. The structural and dynamical properties of Al clusters on Ni surfaces were found to be strongly dependent on the impact energy. At much lower energy, the Al cluster deposited on the surface as a bulk molecule. However, the original icosahedral structure was transformed to the fcc-like one due to the interaction and the structure mismatch between the Al cluster and Ni surface. With increasing the impinging energy, the cluster was deformed severely when it contacted the substrate, and then broken up due to dense collision cascade. The cluster atoms spread on the surface at last. When the impact energy was higher than 11 eV, the defects, such as Al substitutions and Ni ejections, were observed. The simulation indicated that there exists an optimum energy range, which is suitable for Al epitaxial growth in layer by layer. In addition, at higher impinging energy, the atomic exchange between Al and Ni atoms will be favourable to surface alloying.

Reduction of Au3+ ions by activated surface atoms of platinum

In this work it is demonstrated that Pt electrodes can be activated by cathodic polarisation in the hydrogen evolution region which makes it prone to oxidation at potentials below that of bulk oxide formation. When an activated Pt electrode is placed in an aqueous HAuCl4 solution the electroless deposition of Au onto the surface of the electrode is observed and confirmed by cyclic voltammetry and XPS measurements. It is demonstrated that the oxidation of active Pt surface atoms provides the driving force for the spontaneous reduction of Au3+ ions into metallic Au to generate a Pt/Au surface which is highly active for the electro-oxidation of ethanol.

From single crystal surfaces to single atoms : investigating active sites in electrocatalysis

Electrocatalytic processes will undoubtedly be at the heart of energising future transportation and technology with the added importance of being able to create the necessary fuels required to do so in an environmentally friendly and cost effective manner. For this to be successful two almost mutually exclusive surface properties need to be reconciled, namely producing highly active/reactive surface sites that exhibit long term stability. This article reviews the various approaches which have been undertaken to study the elusive nature of these active sites on metal surfaces which are considered as adatoms or clusters of adatoms with low coordination number. This includes the pioneering studies at extended well defined stepped single crystal surfaces using cyclic voltammetry up to the highly sophisticated in situ electrochemical imaging techniques used to study chemically synthesised nanomaterials. By combining the information attained from single crystal surfaces, individual nanoparticles of defined size and shape, density functional theory calculations and new concepts such as mesoporous multimetallic thin films and single atom electrocatalysts new insights into the design and fabrication of materials with highly active but stable active sites can be achieved. The area of electrocatalysis is therefore not only a fascinating and exciting field in terms of realistic technological and economical benefits but also from the fundamental understanding that can be acquired by studying such an array of interesting materials.

Location of hydrogen atoms in hydronium jarosite

Various models for the crystal structure of hydronium jarosite were determined from Rietveld refinements against neutron powder diffraction patterns collected at ambient temperature and also single-crystal X-ray diffraction data. The possibility of a lower symmetry space group for hydronium jarosite that has been suggested by the literature was investigated. It was found the space group is best described as R3¯m, the same for other jarosite minerals. The hydronium oxygen atom was found to occupy the 3¯m site (3a Wyckoff site). Inadequately refined hydronium bond angles and bond distances without the use of restraints are due to thermal motion and disorder of the hydronium hydrogen atoms across numerous orientations. However, the acquired data do not permit a precise determination of these orientations; the main feature up/down disorder of hydronium is clear. Thus, the highest symmetry model with the least disorder necessary to explain all data was chosen: The hydronium hydrogen atoms were modeled to occupy an m (18 h Wyckoff site) with 50 % fractional occupancy, leading to disorder across two orientations. A rigid body description of the hydronium ion rotated by 60° with H–O–H bond angles of 112° and O–H distances of 0.96 Å was optimal. This rigid body refinement suggests that hydrogen bonds between hydronium hydrogen atoms and basal sulfate oxygen atoms are not predominant. Instead, hydrogen bonds are formed between hydronium hydrogen atoms and hydroxyl oxygen atoms. The structure of hydronium alunite is expected to be similar given that alunite supergroup minerals are isostructural.

Photophysical characterization of light-emitting poly(indenofluorene)s

Time-resolved photoluminescence spectroscopy experiments of three poly(2,8-indenofluorene) derivatives bearing different pendant groups are presented. A comparison of the photophysical properties of dilute solutions and thin films provides information on the chemical purity of the materials. The photophysical properties of poly(2,8-indenofluorene)s are correlated with the morphological characteristics of their corresponding films. Wide-angle X-ray scattering experiments reveal the order in these materials at the molecular level. The spectroscopic results confirm the positive impact of a new synthetic approach on the spectral purity of the poly(indenofluorene)s. It is concluded that complete side-chain substitution of the bridgehead carbon atoms C-6 and C-12 in the indenofluorene unit, prior to indenofluorene ring formation, reduces the probability of keto formation. Due to the intrinsic chemical purity of the arylated derivative, identification of a long-delayed spectral feature, other than the known keto band, is possible in the case of thin films. Controlled doping experiments on the arylated derivative with trace amounts of an indenofluorene-monoketone provide quantitative information on the rates of two major photophysical processes, namely, singlet photoluminescence emission and singlet photoluminescence quenching. These results allow the determination of the minimum keto concentration that can affect the intrinsic photophysical properties of this polymer. The data suggest that photoluminescence quenching operates in the doped films according to the Stern-Volmer formalism.

A novel population balance model for the dilute acid hydrolysis of hemicellulose

Acid hydrolysis is a popular pretreatment for removing hemicellulose from lignocelluloses in order to produce a digestible substrate for enzymatic saccharification. In this work, a novel model for the dilute acid hydrolysis of hemicellulose within sugarcane bagasse is presented and calibrated against experimental oligomer profiles. The efficacy of mathematical models as hydrolysis yield predictors and as vehicles for investigating the mechanisms of acid hydrolysis is also examined. Experimental xylose, oligomer (degree of polymerisation 2 to 6) and furfural yield profiles were obtained for bagasse under dilute acid hydrolysis conditions at temperatures ranging from 110C to 170C. Population balance kinetics, diffusion and porosity evolution were incorporated into a mathematical model of the acid hydrolysis of sugarcane bagasse. This model was able to produce a good fit to experimental xylose yield data with only three unknown kinetic parameters ka, kb and kd. However, fitting this same model to an expanded data set of oligomeric and furfural yield profiles did not successfully reproduce the experimental results. It was found that a ``hard-to-hydrolyse'' parameter, $\alpha$, was required in the model to ensure reproducibility of the experimental oligomer profiles at 110C, 125C and 140C. The parameters obtained through the fitting exercises at lower temperatures were able to be used to predict the oligomer profiles at 155C and 170C with promising results. The interpretation of kinetic parameters obtained by fitting a model to only a single set of data may be ambiguous. Although these parameters may correctly reproduce the data, they may not be indicative of the actual rate parameters, unless some care has been taken to ensure that the model describes the true mechanisms of acid hydrolysis. It is possible to challenge the robustness of the model by expanding the experimental data set and hence limiting the parameter space for the fitting parameters. The novel combination of ``hard-to-hydrolyse'' and population balance dynamics in the model presented here appears to stand up to such rigorous fitting constraints.

Limitations of a laboratory scale model in predicting optimal pilot scale conditions for dilute acid pretreatment of sugarcane bagasse

Pilot and industrial scale dilute acid pretreatment data can be difficult to obtain due to the significant infrastructure investment required. Consequently, models of dilute acid pretreatment by necessity use laboratory scale data to determine kinetic parameters and make predictions about optimal pretreatment conditions at larger scales. In order for these recommendations to be meaningful, the ability of laboratory scale models to predict pilot and industrial scale yields must be investigated. A mathematical model of the dilute acid pretreatment of sugarcane bagasse has previously been developed by the authors. This model was able to successfully reproduce the experimental yields of xylose and short chain xylooligomers obtained at the laboratory scale. In this paper, the ability of the model to reproduce pilot scale yield and composition data is examined. It was found that in general the model over predicted the pilot scale reactor yields by a significant margin. Models that appear very promising at the laboratory scale may have limitations when predicting yields on a pilot or industrial scale. It is difficult to comment whether there are any consistent trends in optimal operating conditions between reactor scale and laboratory scale hydrolysis due to the limited reactor datasets available. Further investigation is needed to determine whether the model has some efficacy when the kinetic parameters are re-evaluated by parameter fitting to reactor scale data, however, this requires the compilation of larger datasets. Alternatively, laboratory scale mathematical models may have enhanced utility for predicting larger scale reactor performance if bulk mass transport and fluid flow considerations are incorporated into the fibre scale equations. This work reinforces the need for appropriate attention to be paid to pilot scale experimental development when moving from laboratory to pilot and industrial scales for new technologies.

A comprehensive mechanistic kinetic model for dilute acid hydrolysis of switchgrass cellulose to glucose, 5-HMF and levulinic acid

Switchgrass was treated by 1% (w/w) H₂SO₄in batch tube reactors at temperatures ranging from 140–220°C for up to 60 minutes. In this study, release patterns of glucose, 5-hydroxymethylfurfural (5-HMF), and levulinic acid from switchgrass cellulose were investigated through a mechanistic kinetic model. The predictions were consistent with the measured products of interest when new parameters reflecting the effects of reaction limitations, such as cellulose crystallinity, acid soluble lignin–glucose complex (ASL–glucose) and humins that cannot be quantitatively analyzed, were included. The new mechanistic kinetic model incorporating these parameters simulated the experimental data with R² above 0.97. Results showed that glucose yield was most sensitive to variations in the parameter regarding the cellulose crystallinity at low temperatures (140–180°C), while the impact of crystallinity on the glucose yield became imperceptible at elevated temperatures (200–220 °C). Parameters related to the undesired products (e.g. ASL–glucose and humins) were the most sensitive factors compared with rate constants and other additional parameters in impacting the levulinic acid yield at elevated temperatures (200–220°C), while their impacts were negligible at 140–180°C. These new findings provide a more rational explanation for the kinetic changes in dilute acid pretreatment performance and suggest that the influences of cellulose crystallinity and undesired products including ASL–glucose and humins play key roles in determining the generation of glucose, 5-HMF and levulinic acid from biomass-derived cellulose.