Atomistic-continuum modeling of ultrafast laser-induced melting of silicon targets

In this work, we present an atomistic-continuum model for simulations of ultrafast laser-induced melting processes in semiconductors on the example of silicon. The kinetics of transient non-equilibrium phase transition mechanisms is addressed with MD method on the atomic level, whereas the laser light absorption, strong generated electron-phonon nonequilibrium, fast heat conduction, and photo-excited free carrier diffusion are accounted for with a continuum TTM-like model (called nTTM). First, we independently consider the applications of nTTM and MD for the description of silicon, and then construct the combined MD-nTTM model. Its development and thorough testing is followed by a comprehensive computational study of fast nonequilibrium processes induced in silicon by an ultrashort laser irradiation. The new model allowed to investigate the effect of laser-induced pressure and temperature of the lattice on the melting kinetics. Two competing melting mechanisms, heterogeneous and homogeneous, were identified in our big-scale simulations. Apart from the classical heterogeneous melting mechanism, the nucleation of the liquid phase homogeneously inside the material significantly contributes to the melting process. The simulations showed, that due to the open diamond structure of the crystal, the laser-generated internal compressive stresses reduce the crystal stability against the homogeneous melting. Consequently, the latter can take a massive character within several picoseconds upon the laser heating. Due to the large negative volume of melting of silicon, the material contracts upon the phase transition, relaxes the compressive stresses, and the subsequent melting proceeds heterogeneously until the excess of thermal energy is consumed. A series of simulations for a range of absorbed fluences allowed us to find the threshold fluence value at which homogeneous liquid nucleation starts contributing to the classical heterogeneous propagation of the solid-liquid interface. A series of simulations for a range of the material thicknesses showed that the sample width we chosen in our simulations (800 nm) corresponds to a thick sample. Additionally, in order to support the main conclusions, the results were verified for a different interatomic potential. Possible improvements of the model to account for nonthermal effects are discussed and certain restrictions on the suitable interatomic potentials are found. As a first step towards the inclusion of these effects into MD-nTTM, we performed nanometer-scale MD simulations with a new interatomic potential, designed to reproduce ab initio calculations at the laser-induced electronic temperature of 18946 K. The simulations demonstrated that, similarly to thermal melting, nonthermal phase transition occurs through nucleation. A series of simulations showed that higher (lower) initial pressure reinforces (hinders) the creation and the growth of nonthermal liquid nuclei. For the example of Si, the laser melting kinetics of semiconductors was found to be noticeably different from that of metals with a face-centered cubic crystal structure. The results of this study, therefore, have important implications for interpretation of experimental data on the kinetics of melting process of semiconductors.

Targeted Stimuli-Responsive Dextran Conjugates for Doxorubicin Delivery to Hepatocytes

A targeted, stimuli-responsive, polymeric drug delivery vehicle is being developed in our lab to help alleviate severe side-effects caused by narrow therapeutic window drugs. Targeting specific cell types or organs via proteins, specifically, lectin-mediated targeting holds potential due to the high specificity and affinity of receptor-ligand interactions, rapid internalization, and relative ease of processing. Dextran, a commercially available, biodegradable polymer has been conjugated to doxorubicin and galactosamine to target hepatocytes in a three-step, one-pot synthesis. The loading of doxorubicin and galactose on the conjugates was determined by absorbance at 485 nm and elemental analysis, respectively. Conjugation efficiency based on the amount loaded of each reactant varies from 20% to 50% for doxorubicin and from 2% to 20% for galactosamine. Doxorubicin has also been attached to dextran through an acid-labile hydrazide bond. Doxorubicin acts by intercalating with DNA in the nuclei of cells. The fluorescence of doxorubicin is quenched when it binds to DNA. This allows a fluorescence-based cell-free assay to evaluate the efficacy of the polymer conjugates where we measure the fluorescence of doxorubicin and the conjugates in increasing concentrations of calf thymus DNA. Fluorescence quenching indicates that our conjugates can bind to DNA. The degree of binding increases with polymer molecular weight and substitution of doxorubicin. In cell culture experiments with hepatocytes, the relative uptake of polymer conjugates was evaluated using flow cytometry, and the killing efficiency was determined using the MTT cell proliferation assay. We have found that conjugate uptake is much lower than that of free doxorubicin. Lower uptake of conjugates may increase the maximum dose of drug tolerated by the body. Also, non-galactosylated conjugate uptake is lower than that of the galactosylated conjugate. Microscopy indicates that doxorubicin localizes almost exclusively at the nucleus, whereas the conjugates are present throughout the cell. Doxorubicin linked to dextran through a hydrazide bond was used to achieve improved killing efficiency. Following uptake, the doxorubicin dissociates from the polymer in an endosomal compartment and diffuses to the nucleus. The LC₅₀ of covalently linked doxorubicin is 7.4 μg/mL, whereas that of hydrazide linked doxorubicin is 4.4 μg/mL.

Evaluación de la efectividad de la marca Bodytech® en el ámbito de la percepción de marca haciendo uso de una herramienta de simulación computacional

El presente trabajo de investigación fue realizado con el propósito de modelar el proceso de percepción de marca a partir del análisis de los componentes provenientes de la marca “Bodytech”, esto con el fin de simular el proceso de percepción de marca y evaluar la efectividad de la misma. El proceso que se modela es el de percepción – razonamiento – acción y se hace con el fin de evaluar los gastos en cada uno de los “componentes” que antes mencionaron Para realizar el análisis se hizo uso de un sistema de simulación basada en agentes, el cual recibe valores de diferentes variables analizadas por medio de tres herramientas: (1) un diagrama de grupos poblacionales, (2) un diagrama de desagregación de los núcleos temáticos de la marca y (3) las conclusiones obtenidas de una entrevista que se realizó a los responsables de gestionar la marca. Dicho proceso se lleva a cabo con el fin de poder determinar los valores relacionados al gasto en cada uno de los núcleos temáticos que llevan al sistema a evaluar la percepción de marca y la efectividad de estos gastos. Posteriormente, basados en los resultados del sistema de simulación, se obtiene un escenario que puede ser entendido y parcialmente predicho que le permitirán a Bodytech tener una herramienta de valoración de percepción de su marca.

Numerical model of solid phase transformations governed by nucleation and growth: microstructure development during isothermal crystallization

A simple numerical model which calculates the kinetics of crystallization involving randomly distributed nucleation and isotropic growth is presented. The model can be applied to different thermal histories and no restrictions are imposed on the time and the temperature dependences of the nucleation and growth rates. We also develop an algorithm which evaluates the corresponding emerging grain-size distribution. The algorithm is easy to implement and particularly flexible, making it possible to simulate several experimental conditions. Its simplicity and minimal computer requirements allow high accuracy for two- and three-dimensional growth simulations. The algorithm is applied to explore the grain morphology development during isothermal treatments for several nucleation regimes. In particular, thermal nucleation, preexisting nuclei, and the combination of both nucleation mechanisms are analyzed. For the first two cases, the universal grain-size distribution is obtained. The high accuracy of the model is stated from its comparison to analytical predictions. Finally, the validity of the Kolmogorov-Johnson-Mehl-Avrami model SSSR, is verified for all the cases studied

Funcions densitat i semblança molecular quàntica: nous desenvolupaments i aplicacions

La present tesi, tot i que emmarcada dins de la teoria de les Mesures Semblança Molecular Quántica (MQSM), es deriva en tres àmbits clarament definits: - La creació de Contorns Moleculars de IsoDensitat Electrònica (MIDCOs, de l'anglès Molecular IsoDensity COntours) a partir de densitats electròniques ajustades. - El desenvolupament d'un mètode de sobreposició molecular, alternatiu a la regla de la màxima semblança. - Relacions Quantitatives Estructura-Activitat (QSAR, de l'anglès Quantitative Structure-Activity Relationships). L'objectiu en el camp dels MIDCOs és l'aplicació de funcions densitat ajustades, ideades inicialment per a abaratir els càlculs de MQSM, per a l'obtenció de MIDCOs. Així, es realitza un estudi gràfic comparatiu entre diferents funcions densitat ajustades a diferents bases amb densitats obtingudes de càlculs duts a terme a nivells ab initio. D'aquesta manera, l'analogia visual entre les funcions ajustades i les ab initio obtinguda en el ventall de representacions de densitat obtingudes, i juntament amb els valors de les mesures de semblança obtinguts prèviament, totalment comparables, fonamenta l'ús d'aquestes funcions ajustades. Més enllà del propòsit inicial, es van realitzar dos estudis complementaris a la simple representació de densitats, i són l'anàlisi de curvatura i l'extensió a macromolècules. La primera observació correspon a comprovar no només la semblança dels MIDCOs, sinó la coherència del seu comportament a nivell de curvatura, podent-se així observar punts d'inflexió en la representació de densitats i veure gràficament aquelles zones on la densitat és còncava o convexa. Aquest primer estudi revela que tant les densitats ajustades com les calculades a nivell ab initio es comporten de manera totalment anàloga. En la segona part d'aquest treball es va poder estendre el mètode a molècules més grans, de fins uns 2500 àtoms. Finalment, s'aplica part de la filosofia del MEDLA. Sabent que la densitat electrònica decau ràpidament al allunyar-se dels nuclis, el càlcul d'aquesta pot ser obviat a distàncies grans d'aquests. D'aquesta manera es va proposar particionar l'espai, i calcular tan sols les funcions ajustades de cada àtom tan sols en una regió petita, envoltant l'àtom en qüestió. Duent a terme aquest procés, es disminueix el temps de càlcul i el procés esdevé lineal amb nombre d'àtoms presents en la molècula tractada. En el tema dedicat a la sobreposició molecular es tracta la creació d'un algorisme, així com la seva implementació en forma de programa, batejat Topo-Geometrical Superposition Algorithm (TGSA), d'un mètode que proporcionés aquells alineaments que coincideixen amb la intuïció química. El resultat és un programa informàtic, codificat en Fortran 90, el qual alinea les molècules per parelles considerant tan sols nombres i distàncies atòmiques. La total absència de paràmetres teòrics permet desenvolupar un mètode de sobreposició molecular general, que proporcioni una sobreposició intuïtiva, i també de forma rellevant, de manera ràpida i amb poca intervenció de l'usuari. L'ús màxim del TGSA s'ha dedicat a calcular semblances per al seu ús posterior en QSAR, les quals majoritàriament no corresponen al valor que s'obtindria d'emprar la regla de la màxima semblança, sobretot si hi ha àtoms pesats en joc. Finalment, en l'últim tema, dedicat a la Semblança Quàntica en el marc del QSAR, es tracten tres aspectes diferents: - Ús de matrius de semblança. Aquí intervé l'anomenada matriu de semblança, calculada a partir de les semblances per parelles d'entre un conjunt de molècules. Aquesta matriu és emprada posteriorment, degudament tractada, com a font de descriptors moleculars per a estudis QSAR. Dins d'aquest àmbit s'han fet diversos estudis de correlació d'interès farmacològic, toxicològic, així com de diverses propietats físiques. - Aplicació de l'energia d'interacció electró-electró, assimilat com a una forma d'autosemblança. Aquesta modesta contribució consisteix breument en prendre el valor d'aquesta magnitud, i per analogia amb la notació de l'autosemblança molecular quàntica, assimilar-la com a cas particular de d'aquesta mesura. Aquesta energia d'interacció s'obté fàcilment a partir de programari mecanoquàntic, i esdevé ideal per a fer un primer estudi preliminar de correlació, on s'utilitza aquesta magnitud com a únic descriptor. - Càlcul d'autosemblances, on la densitat ha estat modificada per a augmentar el paper d'un substituent. Treballs previs amb densitats de fragments, tot i donar molt bons resultats, manquen de cert rigor conceptual en aïllar un fragment, suposadament responsable de l'activitat molecular, de la totalitat de l'estructura molecular, tot i que les densitats associades a aquest fragment ja difereixen degut a pertànyer a esquelets amb diferents substitucions. Un procediment per a omplir aquest buit que deixa la simple separació del fragment, considerant així la totalitat de la molècula (calcular-ne l'autosemblança), però evitant al mateix temps valors d'autosemblança no desitjats provocats per àtoms pesats, és l'ús de densitats de Forats de fermi, els quals es troben definits al voltant del fragment d'interès. Aquest procediment modifica la densitat de manera que es troba majoritàriament concentrada a la regió d'interès, però alhora permet obtenir una funció densitat, la qual es comporta matemàticament igual que la densitat electrònica regular, podent-se així incorporar dins del marc de la semblança molecular. Les autosemblances calculades amb aquesta metodologia han portat a bones correlacions amb àcids aromàtics substituïts, podent així donar una explicació al seu comportament. Des d'un altre punt de vista, també s'han fet contribucions conceptuals. S'ha implementat una nova mesura de semblança, la d'energia cinètica, la qual consisteix en prendre la recentment desenvolupada funció densitat d'energia cinètica, la qual al comportar-se matemàticament igual a les densitats electròniques regulars, s'ha incorporat en el marc de la semblança. A partir d'aquesta mesura s'han obtingut models QSAR satisfactoris per diferents conjunts moleculars. Dins de l'aspecte del tractament de les matrius de semblança s'ha implementat l'anomenada transformació estocàstica com a alternativa a l'ús de l'índex Carbó. Aquesta transformació de la matriu de semblança permet obtenir una nova matriu no simètrica, la qual pot ser posteriorment tractada per a construir models QSAR.

Anthropogenic snowfall events in the UK: examples of urban weather modification?

Snow in the UK is generally associated with synoptic or mesoscale weather systems, thus snowfall during quiescent anticyclonic conditions is surprising and might not even be forecast. Consequently it could present a hazard. Snowfall during anticyclonic freezing fog conditions at Didcot and Hereford in December 2006 is investigated here. These two snowfalls seem to present circumstances in which anthropogenically-produced aerosols could have provided ice nuclei within the freezing fog, and therefore might provide characteristic examples of Anthropogenic Snowfall Events (ASEs).

Evaluation of a large-eddy model simulation of a mixed-phase altocumulus cloud using microwave radiometer, lidar and Doppler radar data

Using the Met Office large-eddy model (LEM) we simulate a mixed-phase altocumulus cloud that was observed from Chilbolton in southern England by a 94 GHz Doppler radar, a 905 nm lidar, a dual-wavelength microwave radiometer and also by four radiosondes. It is important to test and evaluate such simulations with observations, since there are significant differences between results from different cloud-resolving models for ice clouds. Simulating the Doppler radar and lidar data within the LEM allows us to compare observed and modelled quantities directly, and allows us to explore the relationships between observed and unobserved variables. For general-circulation models, which currently tend to give poor representations of mixed-phase clouds, the case shows the importance of using: (i) separate prognostic ice and liquid water, (ii) a vertical resolution that captures the thin layers of liquid water, and (iii) an accurate representation the subgrid vertical velocities that allow liquid water to form. It is shown that large-scale ascents and descents are significant for this case, and so the horizontally averaged LEM profiles are relaxed towards observed profiles to account for these. The LEM simulation then gives a reasonable. cloud, with an ice-water path approximately two thirds of that observed, with liquid water at the cloud top, as observed. However, the liquid-water cells that form in the updraughts at cloud top in the LEM have liquid-water paths (LWPs) up to half those observed, and there are too few cells, giving a mean LWP five to ten times smaller than observed. In reality, ice nucleation and fallout may deplete ice-nuclei concentrations at the cloud top, allowing more liquid water to form there, but this process is not represented in the model. Decreasing the heterogeneous nucleation rate in the LEM increased the LWP, which supports this hypothesis. The LEM captures the increase in the standard deviation in Doppler velocities (and so vertical winds) with height, but values are 1.5 to 4 times smaller than observed (although values are larger in an unforced model run, this only increases the modelled LWP by a factor of approximately two). The LEM data show that, for values larger than approximately 12 cm s(-1), the standard deviation in Doppler velocities provides an almost unbiased estimate of the standard deviation in vertical winds, but provides an overestimate for smaller values. Time-smoothing the observed Doppler velocities and modelled mass-squared-weighted fallspeeds shows that observed fallspeeds are approximately two-thirds of the modelled values. Decreasing the modelled fallspeeds to those observed increases the modelled IWC, giving an IWP 1.6 times that observed.

Linear nucleation of polymers

The use of high-melting fibres as linear nuclei for quiescent polymeric melts is instrumental in providing the superior mechanical properties of polymeric self-composites. It also has inherent advantages in the elucidation of fundamental aspects of polymeric crystallization and self-organization, not least in allowing systematic microscopic studies of polymeric crystallization from nucleation through to the growth interface. This has demonstrated explicitly that lamellae develop in two distinct ways, for slower and faster growth, depending on whether fold packing has or has not time to order before the next molecular layer is added with only the former leading to banded growth in linear polyethylene. Other gains in understanding concern cellulation and morphological instability, internuclear interference, isothermal lamellar thickening and banded growth being a consequence of the partial relief of initial surface stress. (c) 2006 Elsevier Ltd. All rights reserved.

Wave functions for the methane molecule

If the potential field due to the nuclei in the methane molecule is expanded in terms of a set of spherical harmonics about the carbon nucleus, only the terms involving s, f, and higher harmonic functions differ from zero in the equilibrium configuration. Wave functions have been calculated for the equilibrium configuration, first including only the spherically symmetric s term in the potential, and secondly including both the s and the f terms. In the first calculation the complete Hartree-Fock S.C.F. wave functions were determined; in the second calculation a variation method was used to determine the best form of the wave function involving f harmonics. The resulting wave functions and electron density functions are presented and discussed

Vibrational intensities in methane

The absorption intensities of the two infra-red active vibrations in methane have been obtained from a perturbation calculation on the equilibrium wave functions derived in the preceding paper. The perturbation field is the change in the potential field due to the nuclei which results from moving the nuclei in the vibrational coordinate concerned, and a simplified form of second order perturbation theory, developed by Pople and Schofield, is used for the calculation. The main approximation involved is the neglect of f and higher harmonics in the spherical harmonic expansion of the nuclear field. The resulting dipole moment derivatives are approximately three times larger than the experimental values, but they show qualitative features and sign relationships which are significant.

Pollen ultrastructure in anther cultures of Datura innoxia. III. Incomplete microspore division

During the microspore division in Datura innoxia, the mitotic spindle is oriented in planes both perpendicular (PE) and oblique (OB) to the spore wall against which the nucleus is situated. However, irrespective of polarity, the usual type of hemispherical wall is laid down at cytokinesis and isolates the generative cell from the rest of the pollen grain (type A). In PE spores the vegetative nucleus initially occupies a central position in the pollen grain, whereas in OB spores the vegetative nucleus lies at the periphery of the grain close to the generative cell. In anther cultures initiated just before the microspore division is due to take place, no marked change can be observed in either orientation or symmetry of the mitotic spindle when the spores divide. In some, however, cytokinesis is disrupted and deposition of the hemispherical wall arrested. In the absence of a complete wall, differentiation of the generative cell cannot take place and binucleate pollen grains are formed having 2 vegetative-type nuclei (type B). The 2 nuclei in the B pollens are always situated against the pollen-grain wall, suggesting that the disruption phenomenon is related to the OB spores. The incomplete wall always makes contact with the intine on the intine-side of the spindle. Wall material may be represented merely as short stubs projecting out from the intine into the cytoplasm, in which event the 2 nuclei lie close to each other and are separated by only a narrow zone of cytoplasm. In other grains the wall is partially developed between the nuclei and terminates at varying distances from the tonoplast; in these, the nuclei are separated by a wider zone of cytoplasm. The significance of these binucleate grains in pollen embryogenesis is discussed.

Cell and nuclear degradation in sweetpotato [Ipomoea batatas (L.) Lam.] root meristems following exposure to salinity

Sodium chloride-induced cell and nuclear degradation in the root meristems of sweetpotato [Ipomoea batatas (L.) Lam.] were determined using fluorescent microscopy and flow cytometry analysis. Two sweetpotato cultivars were grown in liquid Murashige and Skoog medium and subjected to 0 mM and 500 mM NaCl, with or without 15 mM CaCl2, for periods up to 24 h. Changes to the nuclei of root meristematic cells showed a similar pattern of damage to the nuclei using both fluorescent microscopy and flow cytometry analysis. Damage occurring after only a few hours was followed by nuclear degradation at 24 h. Flow cytometry histograms showed a reduction in G1 and G2 nuclei and an increase in degraded nuclei in NaCl-stressed roots. Salinity-induced nuclear degradation was alleviated by the addition of CaCl2.

Microgemma vivaresi n. sp (Microsporidia, Tetramicridae), infecting liver and skeletal muscle of sea scorpions, Taurulus bubalis (Euphrasen 1786) (Osteichthyes, Cottidae), an inshore, littoral fish

The ultrastructure of a new microsporidian species Microgemmia vivaresi n. sp. causing liver cell xenoma formation in sea scorpions, Taurulus bubalis, is described. Stages of merogony, sporogony, and sporogenesis are mixed in the central cytoplasm of developing xenomas. All stages have unpaired nuclei. Uninucleate and multinucleate meronts lie within vacuoles formed from host endoplasmic reticulum and divide by binary or multiple fission. Sporonts, no longer in vacuoles, deposit plaques of surface coat on the plasma membrane that cause the surface to pucker. Division occurs at the Puckered stage into sporoblast mother cells, on which plaques join up to complete the surface coat. A final binary fission gives rise to sporoblasts. A dense globule, thought to be involved in polar tube synthesis, is gradually dispersed during spore maturation. Spores are broadly ovoid, have a large posterior vacuole, and measure 3.6 mu m x 2.1 pint (fresh). The polar tube has a short wide anterior section that constricts abruptly, then runs posteriad to coil about eight times around the posterior vacuole with granular contents. The polaroplast has up to 40 membranes arranged in pairs mostly attached to the wide region of the polar tube and directed posteriorty around a cytoplasm of a coarsely granular appearance. The species is placed alongside the type species Microgemmia hepaticus Ralphs and Matthews 1986 within the family Tetramicridae, which is transferred from the class Dihaplophasea to the class Haplophasea, as there is no evidence for the occurrence of a diplokaryotic phase.

Evidence for a reversal with age in the pattern of near-wall blood flow around aortic branches

The changes that occur with age in the distribution of atherosclerotic lesions around arterial branch points challenge accepted theories relating disease to haemodynamic stresses. We investigated whether flow near branch points changes with age in a way that can account for the different lesion distributions. Flow around 20 branches from immature and mature aortas was investigated by examining the length:width ratio and orientation of endothelial nuclei; these properties depend on the magnitude and direction of near-wall flows, respectively. There were significant changes in the pattern of nuclear shape with age, consistent with a reversal in the pattern of shear around branches. In control regions away from branches, there were no such changes. The role of haemodynamic stresses in atherogenesis may require re-evaluation in the light of these results. (C) 2003 Elsevier Ireland Ltd. All rights reserved.

Time course of gag protein assembly in HIV-1-infected cells: a study by immunoelectron microscopy

Recent biochemical studies have identified high molecular complexes of the HIV Gag precursor in the cytosol of infected cells. Using immunoelectron microscopy we studied the time course of the synthesis and assembly of a HIV Gag precursor protein (pr55gag) in Sf9 cells infected with recombinant baculovirus expressing the HIV gag gene. We also immunolabeled for pr55gag human T4 cells acutely or chronically infected with HIV-1. In Sf9 cells, the time course study showed that the first Gag protein appeared in the cytoplasm at 28-30 h p.i. and that budding started 6-8 h later. Colloidal gold particles, used to visualize the Gag protein, were first scattered randomly throughout the cytoplasm, but soon clusters representing 100 to 1000 copies of pr55gag were also observed. By contrast, in cells with budding or released virus-like particles the cytoplasm was virtually free of gold particles while the released virus-like particles were heavily labeled. Statistical analysis showed that between 80 and 90% of the gold particles in the cytoplasm were seen as singles, as doublets, or in small groups of up to five particles probably representing small oligomers. Clusters of gold particles were also observed in acutely infected lymphocytes as well as in multinuclear cells of chronically infected cultures of T4 cells. In a few cases small aggregates of gold particles were found in the nuclei of T4 lymphocytes. These observations suggest that the Gag polyprotein forms small oligomers in the cytoplasm of expressing cells but that assembly into multimeric complexes takes place predominantly at the plasma membrane. Large accumulations of Gag protein in the cytoplasm may represent misfolded molecules destined for degradation.