Analysis of the convergence of the 1/t and Wang-Landau algorithms in the calculation of multidimensional integrals

In this communication, the convergence of the 1/t and Wang - Landau algorithms in the calculation of multidimensional numerical integrals is analyzed. Both simulation methods are applied to a wide variety of integrals without restrictions in one, two and higher dimensions. The errors between the exact and the calculated values of the integral are obtained and the efficiency and accuracy of the methods are determined by their dynamical behavior. The comparison between both methods and the simple sampling Monte Carlo method is also reported. It is observed that the time dependence of the errors calculated with 1/t algorithm goes as N^{-1/2} (with N the MC trials) in quantitative agreement with the simple sampling Monte Carlo method. It is also showed that the error for the Wang - Landau algorithm saturates in time evidencing the non-convergence of the methods. The sources for the error are also determined.Comment: 8 pages, 5 figure

Inversion of electrical conductivity data with Tikhonov regularization approach: some considerations

Electromagnetic induction measurements, which are generally used to determine lateral variations of apparent electrical conductivity, can provide quantitative estimates of the subsurface conductivity at different depths. Quantitative inference about the Earth’s interior from experimental data is, however, an ill-posed problem. Using the generalised McNeill’s theory for the EM38 ground conductivity meter, we generated synthetic apparent conductivity curves (input data vector) simulating measurements at different heights above the soil surface. The electrical conductivity profile (the Earth model) was then estimated solving a least squares problem with Tikhonov regularization optimised with a projected conjugate gradient algorithm. Although the Tikhonov approach improves the conditioning of the resulting linear system, profile reconstruction can be surprisingly far from the desired true one. On the contrary, the projected conjugate gradient provided the best solution without any explicit regularization (a = 0) of the objective function of the least squares problem. Also, if the initial guess belongs to the image of the system matrix, Im(A), we found that it provides a unique solution in the same subspace Im(A)

The circle of Gánovce. Natural history of an endocast

The natural endocranial cast of Gánovce was found in Slovakia in 1926, and then stored in the National Museum (Národní Muzeum) of Prague. The endocast was extensively studied by Emanuel Vlček (1925-2006), mostly during the 50s and 60s of the past century (e.g., Vlček, 1949, 1953, 1955, 1969), with a large set of analytical tools including radiographic and biochemical sur- veys, and outline shape analysis. He recognized the Neanderthal morphology of the cast, which was dated to 105 ka and which has an estimated volume of 1320cc. In particular, Vlček noticed a similarity with specimens such as Krapina 3, Gibraltar 1 and Saccopastore 1. In fact, these three specimens all display a similar overall cra- nial anatomy, being possibly representatives of an “early and small-brained” Neanderthal morpho- type (Bruner & Manzi, 2006, 2008)

Inversion of electrical conductivity data with Tikhonov regularization approach: some considerations

Electromagnetic induction measurements, which are generally used to determine lateral variations of apparent electrical conductivity, can provide quantitative estimates of the subsurface conductivity at different depths. Quantitative inference about the Earth's interior from experimental data is, however, an ill-posed problem. Using the generalised McNeill's theory for the EM38 ground conductivity meter, we generated synthetic apparent conductivity curves (input data vector) simulating measurements at different heights above the soil surface. The electrical conductivity profile (the Earth model) was then estimated solving a least squares problem with Tikhonov regularization optimised with a projected conjugate gradient algorithm. Although the Tikhonov approach improves the conditioning of the resulting linear system, profile reconstruction can be surprisingly far from the desired true one. On the contrary, the projected conjugate gradient provided the best solution without any explicit regularization ( a= 0) of the objective function of the least squares problem. Also, if the initial guess belongs to the image of the system matrix, Im(A), we found that it provides a unique solution in the same subspace Im(A)

Daytime turbulent exchange between the Amazon forest and the atmosphere

Detailed observations of turbulence just above and below the crown of the Amazon rain forest during the wet season are presented. The forest canopy is shown to remove high frequency turbulent fluctuations while passing lower frequencies. Filter characteristics of turbulent transfer into the Amazon rain forest canopy are quantified. Simple empirical relations that relate observed turbulent heat fluxes to horizontal wind variance are presented. Changes in the amount of turbulent coupling between the forest and the boundary layer associated with deep convective clouds are presented both as statistical averages and as a series of case studies. These convective processes during the rainy season are shown to alter the diurnal course of turbulent fluxes. In wake of giant coastal systems, no significant heat or moisture fluxes occur for up to a day after the event. Radar data is used to demonstrate that even small raining clouds are capable of evacuating the canopy of substances normally trapped by persistent static stability near the forest floor. Recovery from these events can take more than an hour, even during mid-day. In spite of the ubiquitous presence of clouds and frequent rain during this season, the average horizontal wind speed spectrum is well described by dry CBL similarity hypotheses originally found to apply in flat terrain

Large K-exciton dynamics in GaN epilayers: the non-thermal and thermal regime

We present a detailed investigation concerning the exciton dynamics in GaN epilayers grown on c-plane sapphire substrates, focussing on the exciton formation and the transition from the nonthermal to the thermal regime. The time-resolved kinetics of LO-phonon replicas is used to address the energy relaxation in the excitonic band. From ps time-resolved spectra we bring evidence for a long lasting non-thermal excitonic distribution which accounts for the rst 50 ps. Such a behavior is con rmed in di erent experimental conditions, both when non-resonant and resonant excitation are used. At low excitation power density the exciton formation and their subsequent thermalization is dominated by impurity scattering rather than by acoustic phonon scattering. The estimate of the average energy of the excitons as a function of delay after the excitation pulse provides information on the relaxation time, which describes the evolution of the exciton population to the thermal regime.Comment: 9 pages,8 figure

Echocardiography combined with cardiopulmonary exercise testing for the prediction of outcome in idiopathic pulmonary arterial hypertension

BACKGROUND: Right ventricular (RV) function is a major determinant of exercise intolerance and outcome in idiopathic pulmonary arterial hypertension (IPAH). The aim of the study was to evaluate the incremental prognostic value of echocardiography of the RV and cardiopulmonary exercise testing (CPET) on long-term prognosis in these patients. METHODS: One hundred-thirty treatment-naïve IPAH patients were enrolled and prospectively followed. Clinical worsening (CW) was defined by a reduction in 6-minute walk distance plus an increase in functional class, or non elective hospitalization for PAH, or death. Baseline evaluation included clinical, hemodynamic, echocardiographic and CPET variables. Cox regression modeling with c-statistic and bootstrapping validation methods were done. RESULTS: During a mean period of 528 ± 304 days, 54 patients experienced CW (53%). Among demographic, clinical and hemodynamic variables at catheterization, functional class and cardiac index were independent predictors of CW (Model-1). With addition of echocardiographic and CPET variables (Model-2), peak O2 pulse (peak VO2/heart rate) and RV fractional area change (RVFAC) independently improved the power of the prognostic model (AUC: 0.81 vs 0.66, respectively; p=0.005). Patients with low RVFAC and low O2 pulse (low RVFAC + low O2 pulse) and high RVFAC+low O2 pulse showed 99.8 and 29.4 increase in the hazard ratio, respectively (relative risk -RR- of 41.1 and 25.3, respectively), compared with high RVFAC+high O2 pulse (p=0.0001). CONCLUSIONS: Echocardiography combined with CPET provides relevant clinical and prognostic information. A combination of low RVFAC and low O2 pulse identifies patients at a particularly high risk of clinical deterioration

Sexual violence in post-conflict Liberia: survivors and their care.

Using routine data from three clinics offering care to survivors of sexual violence (SV) in Monrovia, Liberia, we describe the characteristics of SV survivors and the pattern of SV and discuss how the current approach could be better adapted to meet survivors' needs. There were 1500 survivors seeking SV care between January 2008 and December 2009. Most survivors were women (98%) and median age was 13 years (Interquartile range: 9-17 years). Sexual aggression occurred during day-to-day activities in 822 (55%) cases and in the survivor's home in 552 (37%) cases. The perpetrator was a known civilian in 1037 (69%) SV events. Only 619 (41%) survivors sought care within 72 h. The current approach could be improved by: effectively addressing the psychosocial needs of child survivors, reaching male survivors, targeting the perpetrators in awareness and advocacy campaigns and reducing delays in seeking care

Total phallic reconstruction after penile amputation for donkey bite: Case report and review of the literature

There are very few reported cases of traumatic amputation of the male genitalia due to animal bite. The management involves thorough washout of the wounds, debridement, antibiotic prophylaxis, tetanus and rabies immunization followed by immediate reconstruction or primary wound closure with delayed reconstruction, when immediate reconstruction is not feasible. When immediate reconstruction is not feasible, long-term good functional and cosmetic results are still possible in the majority of cases by performing total phallic reconstruction. In particular, it is now possible to fashion a cosmetically acceptable sensate phallus with incorporated neourethra, to allow the patient to void while standing and to ejaculate, and with enough bulk to allow the insertion of a penile prosthesis to guarantee the rigidity necessary to engage in penetrative sexual intercourse