Atomic harmonic generation in time-dependent R-matrix theory

We have developed the capability to determine accurate harmonic spectra for multielectron atoms within time-dependent R-matrix (TDRM) theory. Harmonic spectra can be calculated using the expectation value of the dipole length, velocity, or acceleration operator. We assess the calculation of the harmonic spectrum from He irradiated by 390-nm laser light with intensities up to 4 x 10(14) W cm(-2) using each form, including the influence of the multielectron basis used in the TDRM code. The spectra are consistent between the different forms, although the dipole acceleration calculation breaks down at lower harmonics. The results obtained from TDRM theory are compared with results from the HELIUM code, finding good quantitative agreement between the methods. We find that bases which include pseudostates give the best comparison with the HELIUM code, but models comprising only physical orbitals also produce accurate results.

Structural effects on the pH-dependent fluorescence of naphthalenic derivatives and consequences for sensing/switching

Naphthalenic compounds are a rich resource for designers of fluorescent sensing/switching/logic systems. The degree of internal charge transfer (ICT) character in the fluorophore excited states can vary from negligible to substantial. Naphthalene-1,8;4,5-diimides (11–13), 1,8-naphthalimides (16) and 4-chloro-1,8-naphthalimides (15) are of the former type. The latter type is represented by the 4-alkylamino-1,8-naphthalimides (1). Whether ICT-based or not, these serve as the fluorophore in ‘fluorophore-spacer-receptor’ switching systems where PET holds sway until the receptor is bound to H+. On the other hand, 4-dialkylamino-1,8-naphthalimides (3–4) show modest H+-induced fluorescence switching unless the 4-dialkylamino group is a part of a small ring (5). Electrostatic destabilization of a non-emissive twisted internal charge transfer (ICT) excited state is the origin of this behaviour. An evolution to the non-emissive twisted ICT excited state is responsible for the weak emission of the model compound 6 (and related structures 7 and 8) across the pH range. Twisted ICT excited states are also implicated in the switch 9 and its model compound 10, which are based on the 6-dialkylamino-3H-benzimidazo[2,1-a]benz[d,e]isoquinolin-3-one fluorophore.

Time-dependent three-dimensional spectrum synthesis for Type Ia supernovae

A Monte Carlo code (artis) for modelling time-dependent three-dimensional spectral synthesis in chemically inhomogeneous models of Type Ia supernova ejecta is presented. Following the propagation of ?-ray photons, emitted by the radioactive decay of the nucleosynthesis products, energy is deposited in the supernova ejecta and the radiative transfer problem is solved self-consistently, enforcing the constraint of energy conservation in the comoving frame. Assuming a photoionization-dominated plasma, the equations of ionization equilibrium are solved together with the thermal balance equation adopting an approximate treatment of excitation. Since we implement a fully general treatment of line formation, there are no free parameters to adjust. Thus, a direct comparison between synthetic spectra and light curves, calculated from hydrodynamic explosion models, and observations is feasible. The code is applied to the well-known W7 explosion model and the results tested against other studies. Finally, the effect of asymmetric ejecta on broad-band light curves and spectra is illustrated using an elliptical toy model. © 2009 RAS.

Rheological and heat transfer behaviour of the ionic liquid, [C(4)mim][NTf2]

Systematic experiments have been carried out on the thermal and rheological behaviour of the ionic liquid, 1-butyl-3-methylimidazolium bis{(trifluoromethyl)sulfonyl} imide, [C(4)mim][NTf2], and, for the first time, on the forced convective heat transfer of an ionic liquid under the laminar flow conditions. The results show that the thermal conductivity of the ionic liquid is similar to 0.13 W m(-1) K-1, which is almost independent of temperature between 25 and 40 degrees C. Rheological measurements show that the [C(4)mim][NTf2] liquid is a Newtonian fluid with its shear viscosity decreasing with increasing temperature according to the exponential law over a temperature range of 20-90 degrees C. The convective heat transfer experiments demonstrate that the thermal entrance length of the ionic liquid is very large due to its high viscosity and low thermal conductivity. The convective heat transfer coefficient is observed to be much lower than that of distilled water under the same conditions. The convective heat transfer data are also found to fit well to the convectional Shah's equation under the conditions of this work. (C) 2007 Elsevier Inc. All rights reserved.

A multiconfigurational time-dependent Hartree-Fock method for excited electronic states. I. General formalism and application to open-shell states

The solution of the time-dependent Schrodinger equation for systems of interacting electrons is generally a prohibitive task, for which approximate methods are necessary. Popular approaches, such as the time-dependent Hartree-Fock (TDHF) approximation and time-dependent density functional theory (TDDFT), are essentially single-configurational schemes. TDHF is by construction incapable of fully accounting for the excited character of the electronic states involved in many physical processes of interest; TDDFT, although exact in principle, is limited by the currently available exchange-correlation functionals. On the other hand, multiconfigurational methods, such as the multiconfigurational time-dependent Hartree-Fock (MCTDHF) approach, provide an accurate description of the excited states and can be systematically improved. However, the computational cost becomes prohibitive as the number of degrees of freedom increases, and thus, at present, the MCTDHF method is only practical for few-electron systems. In this work, we propose an alternative approach which effectively establishes a compromise between efficiency and accuracy, by retaining the smallest possible number of configurations that catches the essential features of the electronic wavefunction. Based on a time-dependent variational principle, we derive the MCTDHF working equation for a multiconfigurational expansion with fixed coefficients and specialise to the case of general open-shell states, which are relevant for many physical processes of interest. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3600397]

Influence of microstructure on the cutting behaviour of silicon

We use molecular dynamics simulation to study the mechanisms of plasticity during cutting of monocrystalline and polycrystalline silicon. Three scenarios are considered: (i) cutting a single crystal silicon workpiece with a single crystal diamond tool, (ii) cutting a polysilicon workpiece with a single crystal diamond tool, and (iii) cutting a single crystal silicon workpiece with a polycrystalline diamond tool. A long-range analytical bond order potential is used in the simulations, providing a more accurate picture of the atomic-scale mechanisms of brittle fracture, ductile plasticity, and structural changes in silicon. The MD simulation results show a unique phenomenon of brittle cracking typically inclined at an angle of 45° to 55° to the cut surface, leading to the formation of periodic arrays of nanogrooves in monocrystalline silicon, which is a new insight into previously published results. Furthermore, during cutting, silicon is found to undergo solid-state directional amorphisation without prior Si-I to Si-II (beta tin) transformation, which is in direct contrast to many previously published MD studies on this topic. Our simulations also predict that the propensity for amorphisation is significantly higher in single crystal silicon than in polysilicon, signifying that grain boundaries eases the material removal process.

Calculated and observed long term performance of Leneghans embankment

This paper presents the finite element (FE) analysis of the consolidation of the foundation of an embankment constructed over soft clay deposit which shows significant time dependent behaviour and was improved with prefabricated vertical drains. To assess the capability of a simple elastic viscoplastic (EVP) model to predict the long term performance of a geotechnical structure constructed on soft soils, a well documented (Leneghans) embankment was analyzed to predict its long term behaviour characteristics. Two fully coupled two dimensional (2D) plane strain FE analyses have been carried out. In one of these, the foundation of the embankment was modelled with a relatively simpler time dependent EVP model and in the other one, for comparison purposes, the foundation soil was modelled with elasto-plastic Modified Cam-clay (MCC) model. Details of the analyses and the results are discussed in comparison with the field performance. Predictions from the creep (EVP) model were found to be better than those from Elasto-plastic (MCC) analysis. However, the creep analysis requires an additional parameter and additional computational time and resources. © 2011 Taylor & Francis.

Analysis of the compressive behaviour of the three-dimensional printed porous titanium for dental implants using a modified cellular solid model

A set of cylindrical porous titanium test samples were produced using the three-dimensional printing and sintering method with samples sintered at 900 °C, 1000 °C, 1100 °C, 1200 °C or 1300 °C. Following compression testing, it was apparent that the stress-strain curves were similar in shape to the curves that represent cellular solids. This is despite a relative density twice as high as what is considered the threshold for defining a cellular solid. As final sintering temperature increased, the compressive behaviour developed from being elastic-brittle to elastic-plastic and while Young's modulus remained fairly constant in the region of 1.5 GPa, there was a corresponding increase in 0.2% proof stress of approximately 40-80 MPa. The cellular solid model consists of two equations that predict Young's modulus and yield or proof stress. By fitting to experimental data and consideration of porous morphology, appropriate changes to the geometry constants allow modification of the current models to predict with better accuracy the behaviour of porous materials with higher relative densities (lower porosity).

A multiconfigurational time-dependent Hartree-Fock method for excited electronic states. II. Coulomb interaction effects in single conjugated polymer chains

Conjugated polymers have attracted considerable attention in the last few decades due to their potential for optoelectronic applications. A key step that needs optimisation is charge carrier separation following photoexcitation. To understand better the dynamics of the exciton prior to charge separation, we have performed simulations of the formation and dynamics of localised excitations in single conjugated polymer strands. We use a nonadiabatic molecular dynamics method which allows for the coupled evolution of the nuclear degrees of freedom and of multiconfigurational electronic wavefunctions. We show the relaxation of electron-hole pairs to form excitons and oppositely charged polaron pairs and discuss the modifications to the relaxation process predicted by the inclusion of the Coulomb interaction between the carriers. The issue of charge photogeneration in conjugated polymers in dilute solution is also addressed. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3600404]

Thermal properties and eutectic behaviour of dapivirine in combination with steroid hormones and other antiretrovirals

Background: Combination drug products can display thermal behaviour that is more complex than for the corresponding single drug products. For example, the contraceptive vaginal ring (VR) Nuvaring contains a eutectic (lowest melting) composition of etonogestrel (ETN) and ethinyl estradiol. Here we report the predisposition of dapivirine (DPV) to form reduced melting/eutectic mixtures when combined with other contraceptive hormones and antiretrovirals, and discuss the implications for development of combination microbicide and multipurpose prevention technology (MPT) products.
Methods: Binary mixtures of DPV with darunavir (DRV), levonorgestrel (LNG), ETN or maraviroc (MVC) were prepared either by physical mixing or by solvent evaporation. Selected binary mixtures were also incorporated into silicone elastomer (SE) VR devices. Thermal behavior of the mixtures was analyzed using differential scanning calorimetry (DSC) operating in standard heating ramp mode (10 °C/min). DSC data were used to construct two component phase diagrams for each binary system.
Results: Drug mixtures typically showed reduced melting transitions for both drug components, with clear evidence for a eutectic mixture at a well-defined intermediate composition. Eutectic temperatures and compositions for the various mixtures were: 40% DPV / 60% ETN - 170°C; 25% DPV / 75% MVC - 172°C; 65% DPV / 35% LNG - 192°C. In each case, the eutectic composition was also detected when the drug mixtures were incorporated into SE VRs. For the DPV/DRV system, the thermal behaviour is complicated by desolvation from the darunavir ethanolate polymorph.
Conclusions: When DPV is combined with small molecular weight hydrophobic drugs, the melting temperature for both drugs is typically reduced to a degree dependent on the composition of the mixture. At specified compositions, a low melting eutectic system results. The formation of eutectic behavior in binary drug systems needs to be carefully characterised in order to define product performance and drug release.

Harmonic generation by noble-gas atoms in the near-IR regime using ab initio time-dependent R-matrix theory

We demonstrate the capability of ab initio time-dependent R-matrix theory to obtain accurate harmonic generation spectra of noble-gas atoms at near-IR wavelengths between 1200 and 1800 nm and peak intensities up to 1.8 × 10^(14) W/cm^(2). To accommodate the excursion length of the ejected electron, we use an angular-momentum expansion up to Lmax=279. The harmonic spectra show evidence of atomic structure through the presence of a Cooper minimum in harmonic generation for Kr, and of multielectron interaction through the giant resonance for Xe. The theoretical spectra agree well with those obtained experimentally.

REDUCED ORDER MODELLING OF THE THERMAL BEHAVIOUR OF AN OFFICE SPACE

Reduced Order Models (ROMs) have proven to be a valid and efficient approach to model the thermal behaviour of building zones. The main issues associated with the use of zonal/lumped models are how to (1) divide the domain (lumps) and (2) evaluate the pa- rameters which characterise the lump-to-lump exchange of energy and momentum. The object of this research is to develop a methodology for the generation of ROMs from CFD models. The lumps of the ROM and their average property values are automatically ex- tracted from the CFD models through user defined constraints. This methodology has been applied to validated CFD models of a zone of the Environmental Research Insti- tute (ERI) Building in University College Cork (UCC). The ROM predicts temperature distribution in the domain with an average error lower than 2%. It is computationally efficient with an execution time of 3.45 seconds. Future steps in this research will be the development of the procedure to automatically extract the parameters which define lump-to-lump energy and momentum exchange. At the moment these parameters are evaluated through the minimisation of a cost function. The ROMs will also be utilised to predict the transient thermal behaviour of the building zone.

Resistance to chemotherapy in ovarian cancer is miRrored by the microRNA miR-433-dependent dysregulation of the cell cycle.

An understanding of the mechanisms underlying the development of resistance to chemotherapy treatment is a gateway to the introduction of novel therapies and improved outcomes for women presenting with ovarian cancer (OC). The desired apoptotic death post-chemotherapy depends on an intact and fully functioning cell cycle machinery.

In this study we demonstrate that stable expression of miR-433 renders OC cells more resistant to paclitaxel treatment. Interestingly, only cells with the highest miR-433 survived paclitaxel suggesting the possible role of miR-433 in cancer recurrence. Importantly, for the first time we demonstrate that miR 433 induces cellular senescence, exemplified by a flattened morphology, the downregulation of phosphorylated Retinoblastoma (p Rb) and increased β galactosidase activity. Surprisingly, miR 433 induced senescence was independent of two well recognised senescent drivers: p21 and p16. Further in silico analysis followed by in vitro experiments identified CKD6 as a novel miR-433 target gene possibly explaining the observed p21 and p16-independent induction of cellular senescence. Another in silico identified miR-433 target gene was CDC27, a protein involved in the regulation of the cell cycle during mitosis. We demonstrate that the overexpression of pre-miR-433 leads to the downregulation of CDC27 in vitro revealing a novel interaction between miR-433 and CDC27, an integral cell cycle regulating protein.

Interestingly, miR-433 expressing cells also demonstrated an ability to impact their tumour microenvironment. We show that miR-433 is present in exosomes released from miR-433 overexpressing and high miR-433 naïve cells. Moreover, growth condition media (GCM) harvested from cells with high miR-433 have higher levels of IL-6 and IL-8, two key cytokines involved in the senescence associated secretory phenotype (SASP). Importantly, GCM from miR-433-enriched cells repressed the growth of co-cultured cells with initial studies showing a GCM-dependent induction of chemoresistance.

In conclusion, data in this study highlights how the aberrant expression miR-433 contributes to chemoresistance in OC cells. We postulate that standard chemotherapy, particularly paclitaxel, used to treat women with OC may have an attenuated ability to kill cells harbouring increased levels of miR-433, allowing for a subsequent chemoresistant phenotype post-therapy.

Electrical conduction and rheological behaviour of composites of poly(epsilon-caprolactone) and MWCNTs

Two mechanisms of conduction were identified from temperature dependent (120 K-340 K) DC electrical resistivity measurements of composites of poly(c-caprolactone) (PCL) and multi-walled carbon nanotubes (MWCNTs). Activation of variable range hopping (VRH) occurred at lower temperatures than that for temperature fluctuation induced tunneling (TFIT). Experimental data was in good agreement with the VRH model in contrast to the TFIT model, where broadening of tunnel junctions and increasing electrical resistivity at T > T-g is a consequence of a large difference in the coefficients of thermal expansion of PCL and MWCNTs. A numerical model was developed to explain this behavior accounting for a thermal expansion effect by supposing the large increase in electrical resistivity corresponds to the larger relative deformation due to thermal expansion associated with disintegration of the conductive MWCNT network. MWCNTs had a significant nucleating effect on PCL resulting in increased PCL crystallinity and an electrically insulating layer between MWCNTs. The onset of rheological percolation at similar to 0.18 vol% MWCNTs was clearly evident as storage modulus, G' and complex viscosity, vertical bar eta*vertical bar increased by several orders of magnitude. From Cole-Cole and Van Gurp-Palmen plots, and extraction of crossover points (G(c)) from overlaying plots of G' and G '' as a function of frequency, the onset of rheological percolation at 0.18 vol% MWCNTs was confirmed, a similar MWCNT loading to that determined for electrical percolation. 

Predicting low velocity impact damage and Compression-After-Impact (CAI) behaviour of composite laminates

Low-velocity impact damage can drastically reduce the residual strength of a composite structure even when the damage is barely visible. The ability to computationally predict the extent of damage and compression-after-impact (CAI) strength of a composite structure can potentially lead to the exploration of a larger design space without incurring significant time and cost penalties. A high-fidelity three-dimensional composite damage model, to predict both low-velocity impact damage and CAI strength of composite laminates, has been developed and implemented as a user material subroutine in the commercial finite element package, ABAQUS/Explicit. The intralaminar damage model component accounts for physically-based tensile and compressive failure mechanisms, of the fibres and matrix, when subjected to a three-dimensional stress state. Cohesive behaviour was employed to model the interlaminar failure between plies with a bi-linear traction–separation law for capturing damage onset and subsequent damage evolution. The virtual tests, set up in ABAQUS/Explicit, were executed in three steps, one to capture the impact damage, the second to stabilize the specimen by imposing new boundary conditions required for compression testing, and the third to predict the CAI strength. The observed intralaminar damage features, delamination damage area as well as residual strength are discussed. It is shown that the predicted results for impact damage and CAI strength correlated well with experimental testing without the need of model calibration which is often required with other damage models.