Monte Carlo simulation of interface alloying

An Ising-like model, with interactions ranging up to next-nearest-neighbor pairs, is used to simulate the process of interface alloying. Interactions are chosen to stabilize an intermediate "antiferromagnetic" ordered structure. The dynamics proceeds exclusively by atom-vacancy exchanges. In order to characterize the process, the time evolution of the width of the intermediate ordered region and the diffusion length is studied. Both lengths are found to follow a power-law evolution with exponents depending on the characteristic features of the model.

Modeling premartensitic effects in Ni2MnGa: A mean-field and Monte Carlo simulation study

ic first-order transition line ending in a critical point. This critical point is responsible for the existence of large premartensitic fluctuations which manifest as broad peaks in the specific heat, not always associated with a true phase transition. The main conclusion is that premartensitic effects result from the interplay between the softness of the anomalous phonon driving the modulation and the magnetoelastic coupling. In particular, the premartensitic transition occurs when such coupling is strong enough to freeze the involved mode phonon. The implication of the results in relation to the available experimental data is discussed.

Finite-size scaling analysis of the avalanches in the three-dimensional Gaussian random-field Ising model with metastable dynamics

A numerical study is presented of the third-dimensional Gaussian random-field Ising model at T=0 driven by an external field. Standard synchronous relaxation dynamics is employed to obtain the magnetization versus field hysteresis loops. The focus is on the analysis of the number and size distribution of the magnetization avalanches. They are classified as being nonspanning, one-dimensional-spanning, two-dimensional-spanning, or three-dimensional-spanning depending on whether or not they span the whole lattice in different space directions. Moreover, finite-size scaling analysis enables identification of two different types of nonspanning avalanches (critical and noncritical) and two different types of three-dimensional-spanning avalanches (critical and subcritical), whose numbers increase with L as a power law with different exponents. We conclude by giving a scenario for avalanche behavior in the thermodynamic limit.

Non-invasive detection of coronary endothelial response to sequential handgrip exercise in coronary artery disease patients and healthy adults.

OBJECTIVES: Our objective is to test the hypothesis that coronary endothelial function (CorEndoFx) does not change with repeated isometric handgrip (IHG) stress in CAD patients or healthy subjects. BACKGROUND: Coronary responses to endothelial-dependent stressors are important measures of vascular risk that can change in response to environmental stimuli or pharmacologic interventions. The evaluation of the effect of an acute intervention on endothelial response is only valid if the measurement does not change significantly in the short term under normal conditions. Using 3.0 Tesla (T) MRI, we non-invasively compared two coronary artery endothelial function measurements separated by a ten minute interval in healthy subjects and patients with coronary artery disease (CAD). METHODS: Twenty healthy adult subjects and 12 CAD patients were studied on a commercial 3.0 T whole-body MR imaging system. Coronary cross-sectional area (CSA), peak diastolic coronary flow velocity (PDFV) and blood-flow were quantified before and during continuous IHG stress, an endothelial-dependent stressor. The IHG exercise with imaging was repeated after a 10 minute recovery period. RESULTS: In healthy adults, coronary artery CSA changes and blood-flow increases did not differ between the first and second stresses (mean % change ±SEM, first vs. second stress CSA: 14.8%±3.3% vs. 17.8%±3.6%, p = 0.24; PDFV: 27.5%±4.9% vs. 24.2%±4.5%, p = 0.54; blood-flow: 44.3%±8.3 vs. 44.8%±8.1, p = 0.84). The coronary vasoreactive responses in the CAD patients also did not differ between the first and second stresses (mean % change ±SEM, first stress vs. second stress: CSA: -6.4%±2.0% vs. -5.0%±2.4%, p = 0.22; PDFV: -4.0%±4.6% vs. -4.2%±5.3%, p = 0.83; blood-flow: -9.7%±5.1% vs. -8.7%±6.3%, p = 0.38). CONCLUSION: MRI measures of CorEndoFx are unchanged during repeated isometric handgrip exercise tests in CAD patients and healthy adults. These findings demonstrate the repeatability of noninvasive 3T MRI assessment of CorEndoFx and support its use in future studies designed to determine the effects of acute interventions on coronary vasoreactivity.

Pure amnesia after unilateral left polar thalamic infarct: topographic and sequential neuropsychological and metabolic (PET) correlations.

A 54-year-old patient who had an isolated small polar thalamic infarct and acute global amnesia with slight frontal type dysfunction but without other neurological dysfunction was studied. Memory improved partially within 8 months. At all stages the impairment was more severe for verbal than non-verbal memory. Autobiographic recollections and newly acquired information tended to be disorganised with respect to temporal order. Procedural memory was unaffected. Both emotional involvement and pleasure in reading were lost. On MRI, the infarct was limited to the left anterior thalamic nuclei and the adjacent mamillothalamic tract. The regional cerebral metabolic rate of glucose (measured with PET) was decreased on the left in the thalamus, amygdala, and posterior cingulate cortex 2 weeks after the infarct, and in the thalamus and posterior cingulate cortex 9 months later. These findings stress the specific role of the left anterior thalamic region in memory and confirm that longlasting amnesia from a thalamic lesion can occur without significant structural damage to the dorsomedial nucleus. Furthermore, they suggest that the anterior thalamic nuclei and possibly their connections with the posterior cingulate cortex play a role in emotional involvement linked to ipsilateral hemispheric functions.

Application of Visualization and Simulation Program to Improve Work Zone Safety and Mobility, January 2010

A previous study sponsored by the Smart Work Zone Deployment Initiative, “Feasibility of Visualization and Simulation Applications to Improve Work Zone Safety and Mobility,” demonstrated the feasibility of combining readily available, inexpensive software programs, such as SketchUp and Google Earth, with standard two-dimensional civil engineering design programs, such as MicroStation, to create animations of construction work zones. The animations reflect changes in work zone configurations as the project progresses, representing an opportunity to visually present complex information to drivers, construction workers, agency personnel, and the general public. The purpose of this study is to continue the work from the previous study to determine the added value and resource demands created by including more complex data, specifically traffic volume, movement, and vehicle type. This report describes the changes that were made to the simulation, including incorporating additional data and converting the simulation from a desktop application to a web application.

Concurrent or Sequential Adjuvant Hormonoradiotherapy after Conservative Surgery for Early-Stage Breast Cancer: Clinical Results of the CO-HO-RT Phase II Randomized Trial.

Purpose: Letrozole (LET) has recently been shown to be superior to tamoxifen for postmenopausal patients (pts). In addition, LET radiosensitizes breast cancer cells in vitro. We conducted a phase II randomized study to evaluate concurrent and sequential radiotherapy (RT)-LET in the adjuvant setting. We present here clinical results with a minimum follow-up of 24 months. Patients and Methods: Postmenopausal pts with early-stage breast cancer were randomized after conservative surgery to either: A) concurrent RT-LET (LET started 3 weeks before the first day of RT) or B) sequential RT-LET (LET started 3 weeks after the end of RT). Whole breast RT was delivered to a total dose of 50 Gy. A 10-16 Gy boost was allowed according to age and pathological prognostic factors. Pts were stratified by center, adjuvant chemotherapy, boost, and radiation-induced CD8 apoptosis (RILA). RILA was performed before RT as previously published (Ozsahin et al. Clin Cancer Res, 2005). An independent monitoring committee reviewed individual safety data. Skin toxicities were evaluated by two different clinicians at each medical visit (CTCAE v3.0). Lung CT-scan and functional pulmonary tests were performed regularly. DNA samples were screened for SNPs in candidate genes as recently published (Azria et al., Clin Cancer Res, 2008). Results: A total of 150 pts were randomized between 01/05 and 02/07. Median follow-up is 26 months (range, 3-40 months). No statistical differences were identified between the two arms in terms of mean age; initial TNM; median surgical bed volume; post surgical breast volume. Chemotherapy and RT boost were delivered in 19% and 38% of pts, respectively. Nodes received 50 Gy in 23% of patients without differences between both arms. During RT and within the first 6 weeks after RT, 10 patients (6.7%) presented grade 3 acute skin dermatitis during RT but no differences were observed between both arms (4 and 6 patients in arm A and B, respectively). At 26 month of follow-up, grade 2 and more radiation-induced subcutaneous fibrosis (RISCF) was present in 4 patients (3%) without any difference between arm A (n = 2) and B (n = 2), p=0.93. In both arms, all patients that presented a RICSF had a RILA lower than 16%. Sensitivity and specificity were 100% and 39%, respectively.No acute lung toxicities were observed and quality of life was good to excellent for all patients.SNPs analyses are still on-going (Pr Rosenstein, NY). Conclusion: Acute and early late grade 2 dermatitis were similar in both arms. The only factor that influenced RISCF was a low radiation-induced lymphocyte apoptosis yield. We confirmed prospectively the capacity of RILA for identifying hypersensitive patients to radiation. Indeed, patients with RILA superior to 16% did not present late effects to radiation and confirmed the first prospective trial we published in 2005 (Ozsahin et al., Clin Cancer Res).

Use of the swat model for hydro-sedimentologic simulation in a small rural watershed

Mathematical models have great potential to support land use planning, with the goal of improving water and land quality. Before using a model, however, the model must demonstrate that it can correctly simulate the hydrological and erosive processes of a given site. The SWAT model (Soil and Water Assessment Tool) was developed in the United States to evaluate the effects of conservation agriculture on hydrological processes and water quality at the watershed scale. This model was initially proposed for use without calibration, which would eliminate the need for measured hydro-sedimentologic data. In this study, the SWAT model was evaluated in a small rural watershed (1.19 km²) located on the basalt slopes of the state of Rio Grande do Sul in southern Brazil, where farmers have been using cover crops associated with minimum tillage to control soil erosion. Values simulated by the model were compared with measured hydro-sedimentological data. Results for surface and total runoff on a daily basis were considered unsatisfactory (Nash-Sutcliffe efficiency coefficient - NSE < 0.5). However simulation results on monthly and annual scales were significantly better. With regard to the erosion process, the simulated sediment yields for all years of the study were unsatisfactory in comparison with the observed values on a daily and monthly basis (NSE values < -6), and overestimated the annual sediment yield by more than 100 %.

Sex differences in urinary levels of several biological indicators of exposure: a simulation study using a compartmental based toxicokinetic model

Toxicokinetic modeling is a useful tool to describe or predict the behavior of a chemical agent in the human or animal organism. A general model based on four compartments was developed in a previous study in order to quantify the effect of human variability on a wide range of biological exposure indicators. The aim of this study was to adapt this existing general toxicokinetic model to three organic solvents, which were methyl ethyl ketone, 1-methoxy-2-propanol and 1,1,1,-trichloroethane, and to take into account sex differences. We assessed in a previous human volunteer study the impact of sex on different biomarkers of exposure corresponding to the three organic solvents mentioned above. Results from that study suggested that not only physiological differences between men and women but also differences due to sex hormones levels could influence the toxicokinetics of the solvents. In fact the use of hormonal contraceptive had an effect on the urinary levels of several biomarkers, suggesting that exogenous sex hormones could influence CYP2E1 enzyme activity. These experimental data were used to calibrate the toxicokinetic models developed in this study. Our results showed that it was possible to use an existing general toxicokinetic model for other compounds. In fact, most of the simulation results showed good agreement with the experimental data obtained for the studied solvents, with a percentage of model predictions that lies within the 95% confidence interval varying from 44.4 to 90%. Results pointed out that for same exposure conditions, men and women can show important differences in urinary levels of biological indicators of exposure. Moreover, when running the models by simulating industrial working conditions, these differences could even be more pronounced. In conclusion, a general and simple toxicokinetic model, adapted for three well known organic solvents, allowed us to show that metabolic parameters can have an important impact on the urinary levels of the corresponding biomarkers. These observations give evidence of an interindividual variablity, an aspect that should have its place in the approaches for setting limits of occupational exposure.

Numerical simulation of gel electrophoresis of DNA knots in weak and strong electric fields.

Gel electrophoresis allows one to separate knotted DNA (nicked circular) of equal length according to the knot type. At low electric fields, complex knots, being more compact, drift faster than simpler knots. Recent experiments have shown that the drift velocity dependence on the knot type is inverted when changing from low to high electric fields. We present a computer simulation on a lattice of a closed, knotted, charged DNA chain drifting in an external electric field in a topologically restricted medium. Using a Monte Carlo algorithm, the dependence of the electrophoretic migration of the DNA molecules on the knot type and on the electric field intensity is investigated. The results are in qualitative and quantitative agreement with electrophoretic experiments done under conditions of low and high electric fields.

Statistical mechanics in the extended Gaussian ensemble

The extended Gaussian ensemble (EGE) is introduced as a generalization of the canonical ensemble. This ensemble is a further extension of the Gaussian ensemble introduced by Hetherington [J. Low Temp. Phys. 66, 145 (1987)]. The statistical mechanical formalism is derived both from the analysis of the system attached to a finite reservoir and from the maximum statistical entropy principle. The probability of each microstate depends on two parameters ß and ¿ which allow one to fix, independently, the mean energy of the system and the energy fluctuations, respectively. We establish the Legendre transform structure for the generalized thermodynamic potential and propose a stability criterion. We also compare the EGE probability distribution with the q-exponential distribution. As an example, an application to a system with few independent spins is presented.

Sequential partitioning: an alternative to understanding size distributions of avalanches in first-order phase transitions

We study the problem of the partition of a system of initial size V into a sequence of fragments s1,s2,s3 . . . . By assuming a scaling hypothesis for the probability p(s;V) of obtaining a fragment of a given size, we deduce that the final distribution of fragment sizes exhibits power-law behavior. This minimal model is useful to understanding the distribution of avalanche sizes in first-order phase transitions at low temperatures.