A hybrid molecular dynamics/fluctuating hydrodynamics method for modelling liquids at multiple scales in space and time

A new 3D implementation of a hybrid model based on the analogy with two-phase hydrodynamics has been developed for the simulation of liquids at microscale. The idea of the method is to smoothly combine the atomistic description in the molecular dynamics zone with the Landau-Lifshitz fluctuating hydrodynamics representation in the rest of the system in the framework of macroscopic conservation laws through the use of a single "zoom-in" user-defined function s that has the meaning of a partial concentration in the two-phase analogy model. In comparison with our previous works, the implementation has been extended to full 3D simulations for a range of atomistic models in GROMACS from argon to water in equilibrium conditions with a constant or a spatially variable function s. Preliminary results of simulating the diffusion of a small peptide in water are also reported.

On the undamped free vibration of a mass interacting with a Hertzian contact stiffness

Three methods are used to determine the natural frequency of undamped free vibration of a mass interacting with a Hertzian contact stiffness. The exact value is determined using the first integral of motion. The harmonic balance method is used on a transformed equation for an approximate solution, and the multiple scales method is used on an approximate equation. The maximum initial displacement avoiding contact loss is also determined, and the corresponding exact natural frequency is also obtained analytically. The methods are evaluated by studying the free vibration of an elastic sphere on a flat rigid surface. © 2011 Elsevier Ltd © 2011 Elsevier Ltd. All rights reserved.

Kinetics of martensitic transitions in Cu-Al-Mn under thermal cycling: Analysis at multiple length scales

In this paper we study the evolution of the kinetic features of the martensitic transition in a Cu-Al-Mn single crystal under thermal cycling. The use of several experimental techniques including optical microscopy, calorimetry, and acoustic emission, has enabled us to perform an analysis at multiple scales. In particular, we have focused on the analysis of avalanche events (associated with the nucleation and growth of martensitic domains), which occur during the transition. There are significant differences between the kinetics at large and small length scales. On the one hand, at small length scales, small avalanche events tend to sum to give new larger events in subsequent loops. On the other hand, at large length scales the large domains tend to split into smaller ones on thermal cycling. We suggest that such different behavior is the necessary ingredient that leads the system to the final critical state corresponding to a power-law distribution of avalanches.

Kinetics of martensitic transitions in Cu-Al-Mn under thermal cycling: Analysis at multiple length scales

In this paper we study the evolution of the kinetic features of the martensitic transition in a Cu-Al-Mn single crystal under thermal cycling. The use of several experimental techniques including optical microscopy, calorimetry, and acoustic emission, has enabled us to perform an analysis at multiple scales. In particular, we have focused on the analysis of avalanche events (associated with the nucleation and growth of martensitic domains), which occur during the transition. There are significant differences between the kinetics at large and small length scales. On the one hand, at small length scales, small avalanche events tend to sum to give new larger events in subsequent loops. On the other hand, at large length scales the large domains tend to split into smaller ones on thermal cycling. We suggest that such different behavior is the necessary ingredient that leads the system to the final critical state corresponding to a power-law distribution of avalanches.

uvby-beta photometry of solar twins The solar colors, model atmospheres, and the T(eff) and metallicity scales

Aims. Solar colors have been determined on the uvby-beta photometric system to test absolute solar fluxes, to examine colors predicted by model atmospheres as a function of stellar parameters (T(eff), log g, [Fe/H]), and to probe zero-points of T(eff) and metallicity scales. Methods. New uvby-beta photometry is presented for 73 solar-twin candidates. Most stars of our sample have also been observed spectroscopically to obtain accurate stellar parameters. Using the stars that most closely resemble the Sun, and complementing our data with photometry available in the literature, the solar colors on the uvby-beta system have been inferred. Our solar colors are compared with synthetic solar colors computed from absolute solar spectra and from the latest Kurucz (ATLAS9) and MARCS model atmospheres. The zero-points of different T(eff) and metallicity scales are verified and corrections are proposed. Results. Our solar colors are (b - y)(circle dot) = 0.4105 +/- 0.0015, m(1,circle dot) = 0.2122 +/- 0.0018, c(1,circle dot) = 0.3319 +/- 0.0054, and beta(circle dot) = 2.5915 +/- 0.0024. The (b - y)(circle dot) and m(1,circle dot) colors obtained from absolute spectrophotometry of the Sun agree within 3-sigma with the solar colors derived here when the photometric zero-points are determined from either the STIS HST observations of Vega or an ATLAS9 Vega model, but the c(1,circle dot) and beta(circle dot) synthetic colors inferred from absolute solar spectra agree with our solar colors only when the zero-points based on the ATLAS9 model are adopted. The Kurucz solar model provides a better fit to our observations than the MARCS model. For photometric values computed from the Kurucz models, (b - y)(circle dot) and m(1,circle dot) are in excellent agreement with our solar colors independently of the adopted zero-points, but for c(1,circle dot) and beta circle dot agreement is found only when adopting the ATLAS9 zero-points. The c(1,circle dot) color computed from both the Kurucz and MARCS models is the most discrepant, probably revealing problems either with the models or observations in the u band. The T(eff) calibration of Alonso and collaborators has the poorest performance (similar to 140 K off), while the relation of Casagrande and collaborators is the most accurate (within 10 K). We confirm that the Ramirez & Melendez uvby metallicity calibration, recommended by Arnadottir and collaborators to obtain [Fe/H] in F, G, and K dwarfs, needs a small (similar to 10%) zero-point correction to place the stars and the Sun on the same metallicity scale. Finally, we confirm that the c(1) index in solar analogs has a strong metallicity sensitivity.

Buckling and post-buckling of extensible rods revisited: A multiple-scale solution

An exact non-linear formulation of the equilibrium of elastic prismatic rods subjected to compression and planar bending is presented, electing as primary displacement variable the cross-section rotations and taking into account the axis extensibility. Such a formulation proves to be sufficiently general to encompass any boundary condition. The evaluation of critical loads for the five classical Euler buckling cases is pursued, allowing for the assessment of the axis extensibility effect. From the quantitative viewpoint, it is seen that such an influence is negligible for very slender bars, but it dramatically increases as the slenderness ratio decreases. From the qualitative viewpoint, its effect is that there are not infinite critical loads, as foreseen by the classical inextensible theory. The method of multiple (spatial) scales is used to survey the post-buckling regime for the five classical Euler buckling cases, with remarkable success, since very small deviations were observed with respect to results obtained via numerical integration of the exact equation of equilibrium, even when loads much higher than the critical ones were considered. Although known beforehand that such classical Euler buckling cases are imperfection insensitive, the effect of load offsets were also looked at, thus showing that the formulation is sufficiently general to accommodate this sort of analysis. (c) 2008 Elsevier Ltd. All rights reserved.

Statin use and peripheral blood progenitor cells mobilization in patients with multiple myeloma.

PURPOSE: Statins have beneficial effects in patients after myocardial infarction and at least part of the benefit results from mobilization of marrow endothelial progenitors to repopulate damaged myocardial tissues. This study examines if statins may have the same effect in mobilizing marrow progenitors to be harvested and subsequently used in high-dose chemotherapy with progenitor cell rescue in multiple myeloma. METHODS: From 2006 to 2012, 86 patients with multiple myeloma were mobilized with the use of G-CSF and were retrospectively analyzed. Patients with other malignancies or mobilized with the use of chemotherapy or with plerixafor were excluded. RESULTS: The median age of the patients was 60 years. 72 patients had received one line of chemotherapy and 14 patients two or more lines of chemotherapy. Twenty patients were taking statins at the time of the harvest while 66 patients were not. In the group of patients taking statins the success rate of first leukapheresis (obtaining the target number of 4 × 10(6) CD34+ cells/kg) was 85 % while in the group not taking statins this rate was 63.6 %. Despite the comparatively small number of patients this difference approached statistical significance (χ (2) = 0.07). CONCLUSION: This retrospective analysis of 86 patients shows for the first time a possible benefit of statins for peripheral blood progenitor cells mobilization in patients with multiple myeloma. Larger studies would be required to clarify the issue. If their effectiveness is confirmed, statins could be a safe and cheaper addition to chemotherapy and plerixafor for peripheral hematopoietic stem cell mobilization.

99mTc-sestamibi imaging and bone marrow karyotyping in the assessment of multiple myeloma and MGUS

AIM: The first pathogenetic step in multiple myeloma is the emergence of a limited number of clonal plasma cells, clinically known as monoclonal gammopathy of undetermined significance (MGUS). Patients with MGUS do not have symptoms or end-organ damage but they do have a 1% annual risk of progression to multiple myeloma or related malignant disorders. With progression of MGUS to multiple myeloma, complex genetic events occur in the neoplastic plasma cell. Karyotyping and fluorescence in-situ hybridization (FISH) were shown to be of prognostic value in patients with multiple myeloma. Tc-sestamibi imaging reflects myeloma disease activity in bone marrow with very high sensitivity and specificity predicting disease evolution. This study was undertaken to evaluate the role of Tc-sestamibi imaging and cytogenetic analysis in prognosis prediction of MGUS and multiple myeloma. METHODS: We enrolled 30 consecutive patients with a confirmed diagnosis of multiple myeloma or MGUS. Bone marrow biopsy and biochemical staging according to the International Staging System (ISS) were performed in all cases. Karyotype analysis and FISH were performed in 11 of 12 patients with MGUS and in 17 of 18 patients with multiple myeloma having adequate metaphases. RESULTS: The karyotype was abnormal in four of 11 MGUS and in six of 17 multiple myeloma. Abnormalities of chromosome 13 were present in one case of MGUS and in six cases of multiple myeloma whereas the involvement of immunoglobulin was observed in one case of multiple myeloma. An abnormal FISH panel was found in four MGUS and nine multiple myeloma patients. All patients with MGUS showed a normal MIBI scan (score 0). Among patients with multiple myeloma only three, all with ISS stage I, showed a normal scan while a positive scan was obtained in others (score range, 1-7). The MIBI uptake was strongly related to the bone marrow plasma cell infiltration and to cytogenetic abnormalities. Particularly, a MIBI uptake score above 5 identified patients with poor prognosis encompassing all stage III multiple myeloma and three of seven stage II multiple myeloma. On the other hand all stage I and II patients having a MIBI score less than 5 showed a good prognosis. CONCLUSION: Both cytogenetic analysis and a MIBI scan add no relevant prognostic information to the ISS in patients with stage I and III multiple myeloma. The MIBI scan was of prognostic value in stage II multiple myeloma patients. Additionally, MIBI imaging may be useful to guide bone marrow biopsy in order to obtain adequate samples for cytogenetic analysis.

Multiple imputations applied to the DREAM3 phosphoproteomics challenge: a winning strategy.

DREAM is an initiative that allows researchers to assess how well their methods or approaches can describe and predict networks of interacting molecules [1]. Each year, recently acquired datasets are released to predictors ahead of publication. Researchers typically have about three months to predict the masked data or network of interactions, using any predictive method. Predictions are assessed prior to an annual conference where the best predictions are unveiled and discussed. Here we present the strategy we used to make a winning prediction for the DREAM3 phosphoproteomics challenge. We used Amelia II, a multiple imputation software method developed by Gary King, James Honaker and Matthew Blackwell[2] in the context of social sciences to predict the 476 out of 4624 measurements that had been masked for the challenge. To chose the best possible multiple imputation parameters to apply for the challenge, we evaluated how transforming the data and varying the imputation parameters affected the ability to predict additionally masked data. We discuss the accuracy of our findings and show that multiple imputations applied to this dataset is a powerful method to accurately estimate the missing data. We postulate that multiple imputations methods might become an integral part of experimental design as a mean to achieve cost savings in experimental design or to increase the quantity of samples that could be handled for a given cost.

Trends in the prevalence, risk and pregnancy outcome of multiple births with congenital anomaly: a registry-based study in 14 European countries 1984-2007.

OBJECTIVE: To assess the public health consequences of the rise in multiple births with respect to congenital anomalies. DESIGN: Descriptive epidemiological analysis of data from population-based congenital anomaly registries. SETTING: Fourteen European countries. POPULATION: A total of 5.4 million births 1984-2007, of which 3% were multiple births. METHODS: Cases of congenital anomaly included live births, fetal deaths from 20 weeks of gestation and terminations of pregnancy for fetal anomaly. MAIN OUTCOME MEASURES: Prevalence rates per 10,000 births and relative risk of congenital anomaly in multiple versus singleton births (1984-2007); proportion prenatally diagnosed, proportion by pregnancy outcome (2000-07). Proportion of pairs where both co-twins were cases. RESULTS: Prevalence of congenital anomalies from multiple births increased from 5.9 (1984-87) to 10.7 per 10,000 births (2004-07). Relative risk of nonchromosomal anomaly in multiple births was 1.35 (95% CI 1.31-1.39), increasing over time, and of chromosomal anomalies was 0.72 (95% CI 0.65-0.80), decreasing over time. In 11.4% of affected twin pairs both babies had congenital anomalies (2000-07). The prenatal diagnosis rate was similar for multiple and singleton pregnancies. Cases from multiple pregnancies were less likely to be terminations of pregnancy for fetal anomaly, odds ratio 0.41 (95% CI 0.35-0.48) and more likely to be stillbirths and neonatal deaths. CONCLUSIONS: The increase in babies who are both from a multiple pregnancy and affected by a congenital anomaly has implications for prenatal and postnatal service provision. The contribution of assisted reproductive technologies to the increase in risk needs further research. The deficit of chromosomal anomalies among multiple births has relevance for prenatal risk counselling.

Analysis and spatial modeling of the indoor radon Swiss data

The present research deals with an important public health threat, which is the pollution created by radon gas accumulation inside dwellings. The spatial modeling of indoor radon in Switzerland is particularly complex and challenging because of many influencing factors that should be taken into account. Indoor radon data analysis must be addressed from both a statistical and a spatial point of view. As a multivariate process, it was important at first to define the influence of each factor. In particular, it was important to define the influence of geology as being closely associated to indoor radon. This association was indeed observed for the Swiss data but not probed to be the sole determinant for the spatial modeling. The statistical analysis of data, both at univariate and multivariate level, was followed by an exploratory spatial analysis. Many tools proposed in the literature were tested and adapted, including fractality, declustering and moving windows methods. The use of Quan-tité Morisita Index (QMI) as a procedure to evaluate data clustering in function of the radon level was proposed. The existing methods of declustering were revised and applied in an attempt to approach the global histogram parameters. The exploratory phase comes along with the definition of multiple scales of interest for indoor radon mapping in Switzerland. The analysis was done with a top-to-down resolution approach, from regional to local lev¬els in order to find the appropriate scales for modeling. In this sense, data partition was optimized in order to cope with stationary conditions of geostatistical models. Common methods of spatial modeling such as Κ Nearest Neighbors (KNN), variography and General Regression Neural Networks (GRNN) were proposed as exploratory tools. In the following section, different spatial interpolation methods were applied for a par-ticular dataset. A bottom to top method complexity approach was adopted and the results were analyzed together in order to find common definitions of continuity and neighborhood parameters. Additionally, a data filter based on cross-validation was tested with the purpose of reducing noise at local scale (the CVMF). At the end of the chapter, a series of test for data consistency and methods robustness were performed. This lead to conclude about the importance of data splitting and the limitation of generalization methods for reproducing statistical distributions. The last section was dedicated to modeling methods with probabilistic interpretations. Data transformation and simulations thus allowed the use of multigaussian models and helped take the indoor radon pollution data uncertainty into consideration. The catego-rization transform was presented as a solution for extreme values modeling through clas-sification. Simulation scenarios were proposed, including an alternative proposal for the reproduction of the global histogram based on the sampling domain. The sequential Gaussian simulation (SGS) was presented as the method giving the most complete information, while classification performed in a more robust way. An error measure was defined in relation to the decision function for data classification hardening. Within the classification methods, probabilistic neural networks (PNN) show to be better adapted for modeling of high threshold categorization and for automation. Support vector machines (SVM) on the contrary performed well under balanced category conditions. In general, it was concluded that a particular prediction or estimation method is not better under all conditions of scale and neighborhood definitions. Simulations should be the basis, while other methods can provide complementary information to accomplish an efficient indoor radon decision making.

Prevalence and risk of Down syndrome in monozygotic and dizygotic multiple pregnancies in Europe: implications for prenatal screening.

OBJECTIVE: To determine risk of Down syndrome (DS) in multiple relative to singleton pregnancies, and compare prenatal diagnosis rates and pregnancy outcome. DESIGN: Population-based prevalence study based on EUROCAT congenital anomaly registries. SETTING: Eight European countries. POPULATION: 14.8 million births 1990-2009; 2.89% multiple births. METHODS: DS cases included livebirths, fetal deaths from 20 weeks, and terminations of pregnancy for fetal anomaly (TOPFA). Zygosity is inferred from like/unlike sex for birth denominators, and from concordance for DS cases. MAIN OUTCOME MEASURES: Relative risk (RR) of DS per fetus/baby from multiple versus singleton pregnancies and per pregnancy in monozygotic/dizygotic versus singleton pregnancies. Proportion of prenatally diagnosed and pregnancy outcome. STATISTICAL ANALYSIS: Poisson and logistic regression stratified for maternal age, country and time. RESULTS: Overall, the adjusted (adj) RR of DS for fetus/babies from multiple versus singleton pregnancies was 0.58 (95% CI 0.53-0.62), similar for all maternal ages except for mothers over 44, for whom it was considerably lower. In 8.7% of twin pairs affected by DS, both co-twins were diagnosed with the condition. The adjRR of DS for monozygotic versus singleton pregnancies was 0.34 (95% CI 0.25-0.44) and for dizygotic versus singleton pregnancies 1.34 (95% CI 1.23-1.46). DS fetuses from multiple births were less likely to be prenatally diagnosed than singletons (adjOR 0.62 [95% CI 0.50-0.78]) and following diagnosis less likely to be TOPFA (adjOR 0.40 [95% CI 0.27-0.59]). CONCLUSIONS: The risk of DS per fetus/baby is lower in multiple than singleton pregnancies. These estimates can be used for genetic counselling and prenatal screening.

Fiedler trees for multiscale surface analysis

In this work we introduce a new hierarchical surface decomposition method for multiscale analysis of surface meshes. In contrast to other multiresolution methods, our approach relies on spectral properties of the surface to build a binary hierarchical decomposition. Namely, we utilize the first nontrivial eigenfunction of the Laplace-Beltrami operator to recursively decompose the surface. For this reason we coin our surface decomposition the Fiedler tree. Using the Fiedler tree ensures a number of attractive properties, including: mesh-independent decomposition, well-formed and nearly equi-areal surface patches, and noise robustness. We show how the evenly distributed patches can be exploited for generating multiresolution high quality uniform meshes. Additionally, our decomposition permits a natural means for carrying out wavelet methods, resulting in an intuitive method for producing feature-sensitive meshes at multiple scales. Published by Elsevier Ltd.

Evaluation of survey methods for sampling anuran species richness in the neotropics.

Species richness is central to ecological theory, with practical applications in conservation, environmental management and monitoring. Several techniques are available for measuring species richness and composition of amphibians in breeding pools, but the relative efficacy of these methods for sampling high-diversity Neotropical amphibian fauna is poorly understood. I evaluated seven studies from south and south-eastern Brazil to compare the relative and combined effectiveness of two methods for measuring species richness at anuran breeding pools: acoustic surveys with visual encounter of adults and dipnet surveys of larvae. I also compared the relative efficacy of each survey method in detecting species with different reproductive modes. Results showed that both survey methods underestimated the number of species when used separately; however, a close approximation of the actual number of species in each breeding pool was obtained when the methods were combined. There was no difference between survey methods in detecting species with different reproductive modes. These results indicate that researchers should employ multiple survey methods that target both adult and larval life history stages in order to accurately assess anuran species richness at breeding pools in the Neotropics.

Comportamento dinâmico não linear e controle de sistemas eletromecânicos em macro e micro escalas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)