Savani: A variable resolution whole-mantle model of anisotropic shear velocity variations based on multiple data sets

International audience

On The Consistency Of Monte Carlo Track Structure Dna Damage Simulations.

Monte Carlo track structures (MCTS) simulations have been recognized as useful tools for radiobiological modeling. However, the authors noticed several issues regarding the consistency of reported data. Therefore, in this work, they analyze the impact of various user defined parameters on simulated direct DNA damage yields. In addition, they draw attention to discrepancies in published literature in DNA strand break (SB) yields and selected methodologies. The MCTS code Geant4-DNA was used to compare radial dose profiles in a nanometer-scale region of interest (ROI) for photon sources of varying sizes and energies. Then, electron tracks of 0.28 keV-220 keV were superimposed on a geometric DNA model composed of 2.7 × 10(6) nucleosomes, and SBs were simulated according to four definitions based on energy deposits or energy transfers in DNA strand targets compared to a threshold energy ETH. The SB frequencies and complexities in nucleosomes as a function of incident electron energies were obtained. SBs were classified into higher order clusters such as single and double strand breaks (SSBs and DSBs) based on inter-SB distances and on the number of affected strands. Comparisons of different nonuniform dose distributions lacking charged particle equilibrium may lead to erroneous conclusions regarding the effect of energy on relative biological effectiveness. The energy transfer-based SB definitions give similar SB yields as the one based on energy deposit when ETH ≈ 10.79 eV, but deviate significantly for higher ETH values. Between 30 and 40 nucleosomes/Gy show at least one SB in the ROI. The number of nucleosomes that present a complex damage pattern of more than 2 SBs and the degree of complexity of the damage in these nucleosomes diminish as the incident electron energy increases. DNA damage classification into SSB and DSB is highly dependent on the definitions of these higher order structures and their implementations. The authors' show that, for the four studied models, different yields are expected by up to 54% for SSBs and by up to 32% for DSBs, as a function of the incident electrons energy and of the models being compared. MCTS simulations allow to compare direct DNA damage types and complexities induced by ionizing radiation. However, simulation results depend to a large degree on user-defined parameters, definitions, and algorithms such as: DNA model, dose distribution, SB definition, and the DNA damage clustering algorithm. These interdependencies should be well controlled during the simulations and explicitly reported when comparing results to experiments or calculations.

Descriptor-and fragment-based QSAR models for a series of Schistosoma mansoni purine nucleoside inhibitors

The enzyme purine nucleoside phosphorylase from Schistosoma mansoni (SmPNP) is an attractive molecular target for the treatment of major parasitic infectious diseases, with special emphasis on its role in the discovery of new drugs against schistosomiasis, a tropical disease that affects millions of people worldwide. In the present work, we have determined the inhibitory potency and developed descriptor- and fragment-based quantitative structure-activity relationships (QSAR) for a series of 9-deazaguanine analogs as inhibitors of SmPNP. Significant statistical parameters (descriptor-based model: r² = 0.79, q² = 0.62, r²pred = 0.52; and fragment-based model: r² = 0.95, q² = 0.81, r²pred = 0.80) were obtained, indicating the potential of the models for untested compounds. The fragment-based model was then used to predict the inhibitory potency of a test set of compounds, and the predicted values are in good agreement with the experimental results

A novel polar-based human face recognition computational model

Motivated by a recently proposed biologically inspired face recognition approach, we investigated the relation between human behavior and a computational model based on Fourier-Bessel (FB) spatial patterns. We measured human recognition performance of FB filtered face images using an 8-alternative forced-choice method. Test stimuli were generated by converting the images from the spatial to the FB domain, filtering the resulting coefficients with a band-pass filter, and finally taking the inverse FB transformation of the filtered coefficients. The performance of the computational models was tested using a simulation of the psychophysical experiment. In the FB model, face images were first filtered by simulated V1- type neurons and later analyzed globally for their content of FB components. In general, there was a higher human contrast sensitivity to radially than to angularly filtered images, but both functions peaked at the 11.3-16 frequency interval. The FB-based model presented similar behavior with regard to peak position and relative sensitivity, but had a wider frequency band width and a narrower response range. The response pattern of two alternative models, based on local FB analysis and on raw luminance, strongly diverged from the human behavior patterns. These results suggest that human performance can be constrained by the type of information conveyed by polar patterns, and consequently that humans might use FB-like spatial patterns in face processing.

Structure of disordered gold-polymer thin films using small angle x-ray scattering

We have investigated the structure of disordered gold-polymer thin films using small angle x-ray scattering and compared the results with the predictions of a theoretical model based on two approaches-a structure form factor approach and the generalized Porod law. The films are formed of polymer-embedded gold nanoclusters and were fabricated by very low energy gold ion implantation into polymethylmethacrylate (PMMA). The composite films span (with dose variation) the transition from electrically insulating to electrically conducting regimes, a range of interest fundamentally and technologically. We find excellent agreement with theory and show that the PMMA-Au films have monodispersive or polydispersive characteristics depending on the implanted ion dose. (C) 2010 American Institute of Physics. [doi:10.1063/1.3493241]

Local atomic structure in tetragonal pure ZrO(2) nanopowders

The local atomic structures around the Zr atom of pure (undoped) ZrO(2) nanopowders with different average crystallite sizes, ranging from 7 to 40 nm, have been investigated. The nanopowders were synthesized by different wet-chemical routes, but all exhibit the high-temperature tetragonal phase stabilized at room temperature, as established by synchrotron radiation X-ray diffraction. The extended X-ray absorption fine structure (EXAFS) technique was applied to analyze the local structure around the Zr atoms. Several authors have studied this system using the EXAFS technique without obtaining a good agreement between crystallographic and EXAFS data. In this work, it is shown that the local structure of ZrO(2) nanopowders can be described by a model consisting of two oxygen subshells (4 + 4 atoms) with different Zr-O distances, in agreement with those independently determined by X-ray diffraction. However, the EXAFS study shows that the second oxygen subshell exhibits a Debye-Waller (DW) parameter much higher than that of the first oxygen subshell, a result that cannot be explained by the crystallographic model accepted for the tetragonal phase of zirconia-based materials. However, as proposed by other authors, the difference in the DW parameters between the two oxygen subshells around the Zr atoms can be explained by the existence of oxygen displacements perpendicular to the z direction; these mainly affect the second oxygen subshell because of the directional character of the EXAFS DW parameter, in contradiction to the crystallographic value. It is also established that this model is similar to another model having three oxygen subshells, with a 4 + 2 + 2 distribution of atoms, with only one DW parameter for all oxygen subshells. Both models are in good agreement with the crystal structure determined by X-ray diffraction experiments.

Charm and longitudinal structure functions within the Kharzeev-Levin-Nardi model

We use the Kharzeev-Levin-Nardi (KLN) model of the low x gluon distributions to fit recent HERA data on F(L) and F(2)(c)(F(2)(b)). Having checked that this model gives a good description of the data, we use it to predict F(L) and F(2)(c) to be measured in a future electron-ion collider. The results are similar to those obtained with the de Florian-Sassot and Eskola-Paukkunen-Salgado nuclear gluon distributions. The conclusion of this exercise is that the KLN model, simple as it is, may still be used as an auxiliary tool to make estimates for both heavy-ion and electron-ion collisions.

Communication structure of cortical networks

Large-scale cortical networks exhibit characteristic topological properties that shape communication between brain regions and global cortical dynamics. Analysis of complex networks allows the description of connectedness, distance, clustering, and centrality that reveal different aspects of how the network's nodes communicate. Here, we focus on a novel analysis of complex walks in a series of mammalian cortical networks that model potential dynamics of information flow between individual brain regions. We introduce two new measures called absorption and driftness. Absorption is the average length of random walks between any two nodes, and takes into account all paths that may diffuse activity throughout the network. Driftness is the ratio between absorption and the corresponding shortest path length. For a given node of the network, we also define four related measurements, namely in-and out-absorption as well as in-and out-driftness, as the averages of the corresponding measures from all nodes to that node, and from that node to all nodes, respectively. We find that the cat thalamo-cortical system incorporates features of two classic network topologies, Erdos-Renyi graphs with respect to in-absorption and in-driftness, and configuration models with respect to out-absorption and out-driftness. Moreover, taken together these four measures separate the network nodes based on broad functional roles (visual, auditory, somatomotor, and frontolimbic).

Dynamics of carbon, biomass, and structure in two Amazonian forests

Amazon forests are potentially globally significant sources or sinks for atmospheric carbon dioxide. In this study, we characterize the spatial trends in carbon storage and fluxes in both live and dead biomass (necromass) in two Amazonian forests, the Biological Dynamic of Forest Fragments Project (BDFFP), near Manaus, Amazonas, and the Tapajos National Forest (TNF) near Santarem, Para. We assessed coarse woody debris (CWD) stocks, tree growth, mortality, and recruitment in ground-based plots distributed across the terra firme forest at both sites. Carbon dynamics were similar within each site, but differed significantly between the sites. The BDFFP and the TNF held comparable live biomass (167 +/- 7.6 MgC.ha(-1) versus 149 +/- 6.0 MgC.ha(-1), respectively), but stocks of CWD were 2.5 times larger at TNF (16.2 +/- 1.5 MgC.ha(-1) at BDFFP, versus 40.1 +/- 3.9 MgC.ha(-1) at TNF). A model of current forest dynamics suggests that the BDFFP was close to carbon balance, and its size class structure approximated a steady state. The TNF, by contrast, showed rapid carbon accrual to live biomass (3.24 +/- 0.22 MgC.ha(-1).a(-1) in TNF, 2.59 +/- 0.16 MgC.ha(-1).a(-1) in BDFFP), which was more than offset by losses from large stocks of CWD, as well as ongoing shifts of biomass among size classes. This pattern in the TNF suggests recovery from a significant disturbance. The net loss of carbon from the TNF will likely last 10 - 15 years after the initial disturbance (controlled by the rate of decay of coarse woody debris), followed by uptake of carbon as the forest size class structure and composition continue to shift. The frequency and longevity of forests showing such disequilibruim dynamics within the larger matrix of the Amazon remains an essential question to understanding Amazonian carbon balance.

The Influence of Different Land Uses on the Structure of Archaeal Communities in Amazonian Anthrosols Based on 16S rRNA and amoA Genes

Soil from the Amazonian region is usually regarded as unsuitable for agriculture because of its low organic matter content and low pH; however, this region also contains extremely rich soil, the Terra Preta Anthrosol. A diverse archaeal community usually inhabits acidic soils, such as those found in the Amazon. Therefore, we hypothesized that this community should be sensitive to changes in the environment. Here, the archaeal community composition of Terra Preta and adjacent soil was examined in four different sites in the Brazilian Amazon under different anthropic activities. The canonical correspondence analysis of terminal restriction fragment length polymorphisms has shown that the archaeal community structure was mostly influenced by soil attributes that differentiate the Terra Preta from the adjacent soil (i.e., pH, sulfur, and organic matter). Archaeal 16S rRNA gene clone libraries indicated that the two most abundant genera in both soils were Candidatus nitrosphaera and Canditatus nitrosocaldus. An ammonia monoxygenase gene (amoA) clone library analysis indicated that, within each site, there was no significant difference between the clone libraries of Terra Preta and adjacent soils. However, these clone libraries indicated there were significant differences between sites. Quantitative PCR has shown that Terra Preta soils subjected to agriculture displayed a higher number of amoA gene copy numbers than in adjacent soils. On the other hand, soils that were not subjected to agriculture did not display significant differences on amoA gene copy numbers between Terra Preta and adjacent soils. Taken together, our findings indicate that the overall archaeal community structure in these Amazonian soils is determined by the soil type and the current land use.

Integrating structural and input design of a 2-DOF high-speed parallel manipulator: A flexible model-based approach

This paper discusses the integrated design of parallel manipulators, which exhibit varying dynamics. This characteristic affects the machine stability and performance. The design methodology consists of four main steps: (i) the system modeling using flexible multibody technique, (ii) the synthesis of reduced-order models suitable for control design, (iii) the systematic flexible model-based input signal design, and (iv) the evaluation of some possible machine designs. The novelty in this methodology is to take structural flexibilities into consideration during the input signal design; therefore, enhancing the standard design process which mainly considers rigid bodies dynamics. The potential of the proposed strategy is exploited for the design evaluation of a two degree-of-freedom high-speed parallel manipulator. The results are experimentally validated. (C) 2010 Elsevier Ltd. All rights reserved.

New prediction methods for CO2 evaporation inside tubes: Part II - An updated general flow boiling heat transfer model based on flow patterns

Corresponding to the updated flow pattern map presented in Part I of this study, an updated general flow pattern based flow boiling heat transfer model was developed for CO2 using the Cheng-Ribatski-Wojtan-Thome [L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside horizontal tubes, Int. J. Heat Mass Transfer 49 (2006) 4082-4094; L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, Erratum to: ""New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside tubes"" [Heat Mass Transfer 49 (21-22) (2006) 4082-4094], Int. J. Heat Mass Transfer 50 (2007) 391] flow boiling heat transfer model as the starting basis. The flow boiling heat transfer correlation in the dryout region was updated. In addition, a new mist flow heat transfer correlation for CO2 was developed based on the CO2 data and a heat transfer method for bubbly flow was proposed for completeness sake. The updated general flow boiling heat transfer model for CO2 covers all flow regimes and is applicable to a wider range of conditions for horizontal tubes: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m(2) s, heat fluxes from 1.8 to 46 kW/m(2) and saturation temperatures from -28 to 25 degrees C (reduced pressures from 0.21 to 0.87). The updated general flow boiling heat transfer model was compared to a new experimental database which contains 1124 data points (790 more than that in the previous model [Cheng et al., 2006, 2007]) in this study. Good agreement between the predicted and experimental data was found in general with 71.4% of the entire database and 83.2% of the database without the dryout and mist flow data predicted within +/-30%. However, the predictions for the dryout and mist flow regions were less satisfactory due to the limited number of data points, the higher inaccuracy in such data, scatter in some data sets ranging up to 40%, significant discrepancies from one experimental study to another and the difficulties associated with predicting the inception and completion of dryout around the perimeter of the horizontal tubes. (C) 2007 Elsevier Ltd. All rights reserved.

Deterministic mathematical model for the prediction of the as-cast macrostructure

A multiphase deterministic mathematical model was implemented to predict the formation of the grain macrostructure during unidirectional solidification. The model consists of macroscopic equations of energy, mass, and species conservation coupled with dendritic growth models. A grain nucleation model based on a Gaussian distribution of nucleation undercoolings was also adopted. At some solidification conditions, the cooling curves calculated with the model showed oscillations (""wiggles""), which prevented the correct prediction of the average grain size along the structure. Numerous simulations were carried out at nucleation conditions where the oscillations are absent, enabling an assessment of the effect of the heat transfer coefficient on the average grain size and columnar-to-equiaxed transition.

Stable Model Predictive Control for Integrating Systems with Optimizing Targets

This paper concern the development of a stable model predictive controller (MPC) to be integrated with real time optimization (RTO) in the control structure of a process system with stable and integrating outputs. The real time process optimizer produces Optimal targets for the system inputs and for Outputs that Should be dynamically implemented by the MPC controller. This paper is based oil a previous work (Comput. Chem. Eng. 2005, 29, 1089) where a nominally stable MPC was proposed for systems with the conventional control approach where only the outputs have set points. This work is also based oil the work of Gonzalez et at. (J. Process Control 2009, 19, 110) where the zone control of stable systems is studied. The new control for is obtained by defining ail extended control objective that includes input targets and zone controller the outputs. Additional decision variables are also defined to increase the set of feasible solutions to the control problem. The hard constraints resulting from the cancellation of the integrating modes Lit the end of the control horizon are softened,, and the resulting control problem is made feasible to a large class of unknown disturbances and changes of the optimizing targets. The methods are illustrated with the simulated application of the proposed,approaches to a distillation column of the oil refining industry.

Valve friction and nonlinear process model closed-loop identification

Among several process variability sources, valve friction and inadequate controller tuning are supposed to be two of the most prevalent. Friction quantification methods can be applied to the development of model-based compensators or to diagnose valves that need repair, whereas accurate process models can be used in controller retuning. This paper extends existing methods that jointly estimate the friction and process parameters, so that a nonlinear structure is adopted to represent the process model. The developed estimation algorithm is tested with three different data sources: a simulated first order plus dead time process, a hybrid setup (composed of a real valve and a simulated pH neutralization process) and from three industrial datasets corresponding to real control loops. The results demonstrate that the friction is accurately quantified, as well as ""good"" process models are estimated in several situations. Furthermore, when a nonlinear process model is considered, the proposed extension presents significant advantages: (i) greater accuracy for friction quantification and (ii) reasonable estimates of the nonlinear steady-state characteristics of the process. (C) 2010 Elsevier Ltd. All rights reserved.