Modeling nitrogen adsorption in spherical pores of siliceous materials by density functional theory

Adsorption of nitrogen in spherical pores of FDU-1 silica at 77 K is considered by means of a nonlocal density functional theory (NLDFT) accounting for a disordered structure of pore walls. Pore size distribution analysis of various FDU-1 samples subject to different temperatures of calcination revealed three distinct groups of pores. The principal group of pores is identified as ordered spherical mesopores connected with each other by smaller interconnecting pores and irregular micropores present in the mesopore walls. To account for the entrances (connecting pores) into spherical mesopores, a concept of solid mass distribution with respect to the apparent density was introduced. It is shown that the introduction of the aforementioned distribution was sufficient to quantitatively describe experimental adsorption isotherms over the entire range of relative pressures spanning six decades.

Scaling and non-scaling of muscle activity, kinematics, and dynamics in sit-ups with different degrees of difficulty

The purpose of this study was to investigate how the CNS adjusts motor patterns for variants of a complex axial movement-the situp. Adjustments were induced by changing the support surface contact and mass distribution of the body. Healthy adults performed straight-legged sit-ups, 3 s in duration, with support added to or removed from the lumbar trunk, or with mass added to the head or to the legs. Each of these interventions either increased or decreased the difficulty of the task. The study addressed the extent to which changes in sit-up difficulty are compensated by scaling of muscle activity, kinematics, and dynamics versus the extent to which they are compensated by changing discretely the motor pattern. The analysis of muscle activity, kinematics, and dynamics focused on the first 30-40% of the sit-up-the trunk flexion phase-since this is the most critical part of the movement. Our results demonstrate that, in some respects, sit-up kinematics and dynamics scaled with difficulty, but in other respects, they did not. Muscle activity also scaled, in many respects, but in more difficult sit-ups, abdominal flexor activity decreased instead of increased. Non-scaling changes in these parameters suggest that complex movements, such as the sit-up, may require discrete changes in motor pattern in order to deal with large loads, which challenge the available leverage. (C) 2005 Elsevier Ltd. All rights reserved.

Computer simulation of radio-frequency methane/hydrogen plasmas and their interaction with GaAs Surfaces

The following thesis describes the computer modelling of radio frequency capacitively coupled methane/hydrogen plasmas and the consequences for the reactive ion etching of (100) GaAs surfaces. In addition a range of etching experiments was undertaken over a matrix of pressure, power and methane concentration. The resulting surfaces were investigated using X-ray photoelectron spectroscopy and the results were discussed in terms of physical and chemical models of particle/surface interactions in addition to the predictions for energies, angles and relative fluxes to the substrate of the various plasma species. The model consisted of a Monte Carlo code which followed electrons and ions through the plasma and sheath potentials whilst taking account of collisions with background neutral gas molecules. The ionisation profile output from the electron module was used as input for the ionic module. Momentum scattering interactions of ions with gas molecules were investigated via different models and compared against results given by quantum mechanical code. The interactions were treated as central potential scattering events and the resulting neutral cascades were followed. The resulting predictions for ion energies at the cathode compared well to experimental ion energy distributions and this verified the particular form of the electrical potentials used and their applicability in the particular geometry plasma cell used in the etching experiments. The final code was used to investigate the effect of external plasma parameters on the mass distribution, energy and angles of all species impingent on the electrodes. Comparisons of electron energies in the plasma also agreed favourably with measurements made using a Langmuir electric probe. The surface analysis showed the surfaces all to be depleted in arsenic due to its preferential removal and the resultant Ga:As ratio in the surface was found to be directly linked to the etch rate. The etch rate was determined by the methane flux which was predicted by the code.

Nonaqueous capillary zone electrophoresis of synthetic organic polypeptides

Poly(Nε-trifluoroacetyl-l-lysine) was used as a model solute to investigate the potential of nonaqueous capillary electrophoresis (NACE) for the characterization of synthetic organic polymers. The information obtained by NACE was compared to that derived from size exclusion chromatography (SEC) experiments, and the two techniques were found to be complimentary for polymer characterization. On one hand, NACE permitted (i) the separation of oligomers according to their molar mass and (ii) the separation of the polymers according to the nature of the end groups. On the other hand, SEC experiments were used for the characterization of the molar mass distribution for higher molar masses. Due to the tendency of the solutes (polypeptides) to adsorb onto the fused-silica capillary wall, careful attention was paid to the rinsing procedure of the capillary between runs in order to keep the capillary surface clean. For that purpose, the use of electrophoretic desorption under denaturating conditions was very effective. Optimization of the separation was performed by studying (i) the influence of the proportion of methanol in a methanol/acetonitrile mixture and (ii) the influence of acetic acid concentration in the background electrolyte. Highly resolved separation of the oligomers (up to a degree of polymerization n of ∼50) was obtained by adding trifluoroacetic acid to the electrolyte. Important information concerning the polymer conformations could be obtained from the mobility data. Two different plots relating the effective mobility data to the degree of polymerization were proposed for monitoring the changes in polymer conformations as a function of the number of monomers.

Database of the GLAMAP project

A uniform chronology for foraminifera-based sea surface temperature records has been established in more than 120 sediment cores obtained from the equatorial and eastern Atlantic up to the Arctic Ocean. The chronostratigraphy of the last 30,000 years is mainly based on published d18O records and 14C ages from accelerator mass spectrometry, converted into calendar-year ages. The high-precision age control provides the database necessary for the uniform reconstruction of the climate interval of the Last Glacial Maximum within the GLAMAP-2000 project.

Prioritizing tropical habitats for long-distance migratory songbirds: an assessment of habitat quality at a stopover site in Colombia

Long-distance migratory birds are declining globally and migration has been identified as the primary source of mortality in this group. Despite this, our lack of knowledge of habitat use and quality at stopovers, i.e., sites where the energy for migration is accumulated, remains a barrier to designing appropriate conservation measures, especially in tropical regions. There is therefore an urgent need to assess stopover habitat quality and concurrently identify efficient and cost-effective methods for doing so. Given that fuel deposition rates directly influence stopover duration, departure fuel load, and subsequent speed of migration, they are expected to provide a direct measure of habitat quality and have the advantage of being measurable through body-mass changes. Here, we examined seven potential indicators of quality, including body-mass change, for two ecologically distinct Neotropical migratory landbirds on stopover in shade-coffee plantations and tropical humid premontane forest during spring migration in Colombia: (1) rate of body-mass change; (2) foraging rate; (3) recapture rate; (4) density; (5) flock size; (6) age and sex ratios; and (7) body-mass distribution. We found higher rates of mass change in premontane forest than in shade-coffee in Tennessee Warbler Oreothlypis peregrina, a difference that was mirrored in higher densities and body masses in forest. In Gray-cheeked Thrush Catharus minimus, a lack of recaptures in shade-coffee and higher densities in forest, also suggested that forest provided superior fueling conditions. For a reliable assessment of habitat quality, we therefore recommend using a suite of indicators, taking into account each species’ ecology and methodological considerations. Our results also imply that birds stopping over in lower quality habitats may spend a longer time migrating and require more stopovers, potentially leading to important carryover effects on reproductive fitness. Evaluating habitat quality is therefore imperative prior to defining the conservation value of newly identified stopover regions.

A Search for a New Particle Decaying to Pairs of Weak Gauge Bosons with the ATLAS Detector

A search for new heavy resonances decaying to boson pairs (WZ, WW or ZZ) using 20.3 inverse femtobarns of proton-proton collision data at a center of mass energy of 8 TeV is presented. The data were recorded by the ATLAS detector at the Large Hadron Collider (LHC) in 2012. The analysis combines several search channels with the leptonic, semi-leptonic and fully hadronic final states. The diboson invariant mass spectrum is studied for local excesses above the Standard Model background prediction, and no significant excess is observed for the combined analysis. 95$\%$ confidence limits are set on the cross section times branching ratios for three signal models: an extended gauge model with a heavy W boson, a bulk Randall-Sundrum model with a spin-2 graviton, and a simplified model with a heavy vector triplet. Among the individual search channels, the fully-hadronic channel is predominantly presented where boson tagging technique and jet substructure cuts are used. Local excesses are found in the dijet mass distribution around 2 TeV, leading to a global significance of 2.5 standard deviations. This deviation from the Standard Model prediction results in many theory explanations, and the possibilities could be further explored using the LHC Run 2 data.

Description of an 8 MW reference wind turbine

An 8 MW wind turbine is described in terms of mass distribution, dimensions, power curve, thrust curve, maximum design load and tower configuration. This turbine has been described as part of the EU FP7 project LEANWIND in order to facilitate research into logistics and naval architecture efficiencies for future offshore wind installations. The design of this 8 MW reference wind turbine has been checked and validated by the design consultancy DNV-GL. This turbine description is intended to bridge the gap between the NREL 5 MW and DTU 10 MW reference turbines and thus contribute to the standardisation of research and development activities in the offshore wind energy industry.

Planktic foraminifera of sediment core GeoB1710-3

In this study, we test various parameters in deep-sea sediments (bulk sediment parameters and changes in microfossil abundances and preservation character) which are generally accepted as indicators of calcium carbonate dissolution. We investigate sediment material from station GeoB 1710-3 in the northern Cape Basin (eastern South Atlantic), 280 km away from the Namibian coast, well outside today's coastal upwelling. As northern Benguela upwelling cells were displaced westward and periodically preceded the core location during the past 245 kyr (Volbers et al., submitted), GeoB 1710-3 sediments reflect these changes in upwelling productivity. Results of the most commonly used calcium carbonate dissolution proxies do not only monitor dissolution within these calcareous sediments but also reflect changes in upwelling intensity. Accordingly, these conventional proxy parameters misrepresent, to some extent, the extent of calcium carbonate dissolution. These results were verified by an independent dissolution proxy, the Globigerina bulloides dissolution index (BDX') (Volbers and Henrich, 2002, doi:10.1016/S0025-3227(02)00333-X). The BDX' is based on scanning electronic microscope ultrastructural investigation of planktonic foraminiferal tests and indicates persistent good carbonate preservation throughout the past 245 kyr, with the exception of one pronounced dissolution event at early oxygen isotopic stage (OIS) 6. The early OIS 6 is characterized by calcium carbonate contents, sand contents, and planktonic foraminiferal concentrations all at their lowest levels for the last 245 kyr. At the same time, the ratio of radiolarian to planktonic foraminiferal abundances and the ratio of benthic to planktonic foraminiferal tests are strongly increased, as are the rain ratio, the fragmentation index, and the BDX'. The sedimentary calcite lysocline rose above the core position and GeoB 1710-3 sediments were heavily altered, as attested to by the unusual accumulation of pellets, aggregates, sponge spicules, radiolaria, benthic foraminifera, and planktonic foraminiferal assemblages. Solely the early OIS 6 dissolution event altered the coarse fraction intensely, and is therefore reflected by all conventional calcium carbonate preservation proxies and the BDX'. We attribute the more than 1000 m rise of the sedimentary calcite lysocline to the combination of two processes: (a) a prominent change in the deep-water mass distribution within the South Atlantic and (b) intense degradation of organic material within the sediment (preserved as maximum total organic carbon content) creating microenvironments favorable for calcium carbonate dissolution.

(Table T1) Abundance of radiolarian skeletons in scrape samples of ODP Hole 178-1098B

Palmer Deep is a series of three glacially overdeepened basins on the Antarctic Peninsula shelf, ~20 km southwest of Anvers Island. Site 1098 (64°51.72'S, 64°12.48'W) was drilled in the shallowest basin, Basin I, at 1012 m water depth. The sediment recovered was primarily laminated, siliceous, biogenic, pelagic muds alternating with siliciclastic hemipelagic sediments (Barker, Camerlenghi, Acton, et al., 1999). Sedimentation rates of 0.1725 cm/yr in the upper 25 m and 0.7-0.80 cm/yr in the lower 25 m of the core have been estimated from 14C (Domack et al., 2001). The oldest datable sediments have an age of ~13 ka and were underlain by diamicton sediments of the last glacial maximum (Domack et al., 2001). The large-scale water-mass distribution and circulation in the vicinity of Palmer Deep is dominated by Circumpolar Deep Water (CDW) below 200 m (Hofmann et al., 1996). Palmer Deep is too far from the coast to be influenced by glacial meltwater and cold-tongue generation associated with it (Domack and Williams, 1990; Dixon and Domack, 1991). Circulation patterns in the Palmer Deep area are not well understood, but evidence suggests southward flow across Palmer Deep from Anvers Island to Renaud Island (Kock and Stein, 1978). The water south of Anvers Island is nearly open with loose pack ice from February through May. The area is covered with sea ice beginning in June (Gloersen et al., 1992; Leventer et al., 1996). Micropaleontologic data from the work of Leventer et al. (1996) on a 9-m piston core has revealed circulation and climate patterns for the past 3700 yr in the Palmer Deep. The benthic foraminifer assemblage is dominated by two taxa, Bulimina aculeata and Bolivina pseudopunctata, which are inversely related. High relative abundances of B. aculeata occur cyclically over a period of ~230 yr. The assemblage associated with high abundance of B. aculeata in Palmer Deep resembles that from the Bellingshausen shelf, which is associated with CDW. In addition to the faunal evidence, hydrographic data indicate incursions of CDW into Palmer Deep (Leventer et al., 1996). A distinctive diatom assemblage dominated by a single genus was associated with peaks in B. aculeata, whereas a few different assemblages were associated with lows in B. aculeata. Leventer et al. (1996) interpreted the variability in diatom assemblages as an indication of changes in productivity associated with changes in water column stability. Abelmann and Gowing (1997) studied the horizontal and vertical distributions of radiolarians in the Atlantic sector of the Southern Ocean. They show that the spatial distribution of radiolarian assemblages reflects hydrographic boundaries. In a transect from the subtropical Atlantic to polar Antarctic zones, radiolarians in the upper 1000 m of the water column occurred in distinct surface and deep-living assemblages related to water depth, temperature, salinity, and nutrient content. Living assemblages resembled those preserved in underlying surface sediments (Abelmann and Gowing, 1997). Circumantarctic coastal sediments from neritic environments contained a distinctive assemblage dominated by the Phormacantha hystrix/Plectacantha oikiskos group and Rhizoplegma boreale (Nishimura et al., 1997). Low diversity and species compositions distinguished the coastal sediments from the typical pelagic Antarctic assemblages. Factors that controlled the assemblages were water depth, proximity to the coast, occurrence of sea ice, and steepness of topography, rather than temperature and salinity. Nishimura et al. (1997) found a gradient of sorts from deep-water sites containing diverse assemblages typical of pelagic environments to coastal sites with low diversity assemblages dominated by P. hystrix/P. oikiskos group and R. boreale. In general, sites between these two extremes had increased proportions of the coastal assemblage with decreasing water depth (Nishimura et al., 1997). At a site near Hole 1098 (GC905), they showed that the relative abundance of the coastal assemblage increased downcore (Nishimura et al., 1997). The purpose of the research presented here was to make a cursory investigation into the radiolarian assemblages as possible paleoenvironmental indicators.

Resolved Stars' Insights into Galaxy Physics

Thesis (Ph.D.)--University of Washington, 2016-08

Experimental and analytical investigation of refrigerant and lubricant migration

The off-cycle refrigerant mass migration has a direct influence on the on-cycle performance since compressor energy is necessary to redistribute the refrigerant mass. No studies, as of today, are available in the open literature which experimentally measured the lubricant migration within a refrigeration system during cycling or stop/start transients. Therefore, experimental procedures measuring the refrigerant and lubricant migration through the major components of a refrigeration system during stop/start transients were developed and implemented. Results identifying the underlying physics are presented. The refrigerant and lubricant migration of an R134a automotive A/C system-utilizing a fixed orifice tube, minichannel condenser, plate and fin evaporator, U-tube type accumulator and fixed displacement compressor-was measured across five sections divided by ball valves. Using the Quick-Closing Valve Technique (QCVT) combined with the Remove and Weigh Technique (RWT) using liquid nitrogen as the condensing agent resulted in a measurement uncertainty of 0.4 percent regarding the total refrigerant mass in the system. The determination of the lubricant mass distribution was achieved by employing three different techniques-Remove and Weigh, Mix and Sample, and Flushing. To employ the Mix and Sample Technique a device-called the Mix and Sample Device-was built. A method to separate the refrigerant and lubricant was developed with an accuracy-after separation-of 0.04 grams of refrigerant left in the lubricant. When applying the three techniques, the total amount of lubricant mass in the system was determined to within two percent. The combination of measurement results-infrared photography and high speed and real time videography-provide unprecedented insight into the mechanisms of refrigerant and lubricant migration during stop-start operation. During the compressor stop period, the primary refrigerant mass migration is caused by, and follows, the diminishing pressure difference across the expansion device. The secondary refrigerant migration is caused by a pressure gradient as a result of thermal nonequilibrium within the system and causes only vapor phase refrigerant migration. Lubricant migration is proportional to the refrigerant mass during the primary refrigerant mass migration. During the secondary refrigerant mass migration lubricant is not migrating. The start-up refrigerant mass migration is caused by an imbalance of the refrigerant mass flow rates across the compressor and expansion device. The higher compressor refrigerant mass flow rate was a result of the entrainment of foam into the U-tube of the accumulator. The lubricant mass migration during the start-up was not proportional to the refrigerant mass migration. The presence of water condensate on the evaporator affected the refrigerant mass migration during the compressor stop period. Caused by an evaporative cooling effect the evaporator held 56 percent of the total refrigerant mass in the system after three minutes of compressor stop time-compared to 25 percent when no water condensate was present on the evaporator coil. Foam entrainment led to a faster lubricant and refrigerant mass migration out of the accumulator than liquid entrainment through the hole at the bottom of the U-tube. The latter was observed for when water condensate was present on the evaporator coil because-as a result of the higher amount of refrigerant mass in the evaporator before start-up-the entrainment of foam into the U-tube of the accumulator ceased before the steady state refrigerant mass distribution was reached.

Search for Higgs Boson Production Beyond the Standard Model Using the Razor Kinematic Variables in pp Collisions at √s=8 TeV and Optimization of Higgs Boson Identification Using a Quantum Annealer

In the first part of this thesis we search for beyond the Standard Model physics through the search for anomalous production of the Higgs boson using the razor kinematic variables. We search for anomalous Higgs boson production using proton-proton collisions at center of mass energy √s=8 TeV collected by the Compact Muon Solenoid experiment at the Large Hadron Collider corresponding to an integrated luminosity of 19.8 fb^-1.

In the second part we present a novel method for using a quantum annealer to train a classifier to recognize events containing a Higgs boson decaying to two photons. We train that classifier using simulated proton-proton collisions at √s=8 TeV producing either a Standard Model Higgs boson decaying to two photons or a non-resonant Standard Model process that produces a two photon final state.

The production mechanisms of the Higgs boson are precisely predicted by the Standard Model based on its association with the mechanism of electroweak symmetry breaking. We measure the yield of Higgs bosons decaying to two photons in kinematic regions predicted to have very little contribution from a Standard Model Higgs boson and search for an excess of events, which would be evidence of either non-standard production or non-standard properties of the Higgs boson. We divide the events into disjoint categories based on kinematic properties and the presence of additional b-quarks produced in the collisions. In each of these disjoint categories, we use the razor kinematic variables to characterize events with topological configurations incompatible with typical configurations found from standard model production of the Higgs boson.

We observe an excess of events with di-photon invariant mass compatible with the Higgs boson mass and localized in a small region of the razor plane. We observe 5 events with a predicted background of 0.54 ± 0.28, which observation has a p-value of 10^-3 and a local significance of 3.35σ. This background prediction comes from 0.48 predicted non-resonant background events and 0.07 predicted SM higgs boson events. We proceed to investigate the properties of this excess, finding that it provides a very compelling peak in the di-photon invariant mass distribution and is physically separated in the razor plane from predicted background. Using another method of measuring the background and significance of the excess, we find a 2.5σ deviation from the Standard Model hypothesis over a broader range of the razor plane.

In the second part of the thesis we transform the problem of training a classifier to distinguish events with a Higgs boson decaying to two photons from events with other sources of photon pairs into the Hamiltonian of a spin system, the ground state of which is the best classifier. We then use a quantum annealer to find the ground state of this Hamiltonian and train the classifier. We find that we are able to do this successfully in less than 400 annealing runs for a problem of median difficulty at the largest problem size considered. The networks trained in this manner exhibit good classification performance, competitive with the more complicated machine learning techniques, and are highly resistant to overtraining. We also find that the nature of the training gives access to additional solutions that can be used to improve the classification performance by up to 1.2% in some regions.

Numeric model to estimate pesticide deposition and losses from aerial spraying.

Pesticides applications have been described by many researches as a very inefficient process. In some cases, there are reports that only 0.02% of the applied products are used for the effective control of the problem. The main factor that influences pesticides applications is the droplet size formed on spraying nozzles. Many parameters affects the dynamic of the droplets, like wind, temperature, relative humidity, and others. Small droplets are biologically more active, but they are affected by evaporation and drift. On the other hand, the great droplets do not promote a good distribution of the product on the target. In this sense, associated with the risk of non target areas contamination and with the high costs involved in applications, the knowledge of the droplet size is of fundamental importance in the application technology. When sophisticated technology for droplets analysis is unavailable, is common the use of artificial targets like water-sensitive paper to sample droplets. On field sampling, water-sensitive papers are placed on the trials where product will be applied. When droplets impinging on it, the yellow surface of this paper will be stained dark blue, making easy their recognition. Collected droplets on this papers have different kinds of sizes. In this sense, the determination of the droplet size distribution gives a mass distribution of the material and so, the efficience of the application of the product. The stains produced by droplets shows a spread factor proportional to their respectives initial sizes. One of methodologies to analyse the droplets is a counting and measure of the droplets made in microscope. The Porton N-G12 graticule, that shows equaly spaces class intervals on geometric progression of square 2, are coulpled to the lens of the microscope. The droplet size parameters frequently used are the Volumetric Median Diameter (VMD) and the Numeric Median Diameter. On VMD value, a representative droplets sample is divided in two equal parts of volume, in such away one part contains droplets of sizes smaller than VMD and the other part contains droplets of sizes greater that VMD. The same process is done to obtaining the NMD, which divide the sample in two equal parts in relation to the droplets size. The ratio between VMD and NMD allows the droplets uniformity evaluation. After that, the graphics of accumulated probability of the volume and size droplets are plotted on log scale paper (accumulated probability versus median diameter of each size class). The graphics provides the NMD on the x-axes point corresponding to the value of 50% founded on the y-axes. All this process is very slow and subjected to operator error. So, in order to decrease the difficulty envolved with droplets measuring it was developed a numeric model, implemented on easy and accessfull computational language, which allows approximate VMD and NMD values, with good precision. The inputs to this model are the frequences of the droplets sizes colected on the water-sensitive paper, observed on the Porton N-G12 graticule fitted on microscope. With these data, the accumulated distribution of the droplet medium volumes and sizes are evaluated. The graphics obtained by plotting this distributions allow to obtain the VMD and NMD using linear interpolation, seen that on the middle of the distributions the shape of the curves are linear. These values are essential to evaluate the uniformity of droplets and to estimate the volume deposited on the observed paper by the density (droplets/cm2). This methodology to estimate the droplets volume was developed by 11.0.94.224 Project of the CNPMA/EMBRAPA. Observed data of herbicides aerial spraying samples, realized by Project on Pelotas/RS county, were used to compare values obtained manual graphic method and with those obtained by model has shown, with great precision, the values of VMD and NMD on each sampled collector, allowing to estimate a quantities of deposited product and, by consequence, the quantities losses by drifty. The graphics of variability of VMD and NMD showed that the quantity of droplets that reachs the collectors had a short dispersion, while the deposited volume shows a great interval of variation, probably because the strong action of air turbulence on the droplets distribution, enfasizing the necessity of a deeper study to verify this influences on drift.

Gravitational waves as dark sirens: an astrophysical and cosmological analysis

Despite the success of the ΛCDM model in describing the Universe, a possible tension between early- and late-Universe cosmological measurements is calling for new independent cosmological probes. Amongst the most promising ones, gravitational waves (GWs) can provide a self-calibrated measurement of the luminosity distance. However, to obtain cosmological constraints, additional information is needed to break the degeneracy between parameters in the gravitational waveform. In this thesis, we exploit the latest LIGO-Virgo-KAGRA Gravitational Wave Transient Catalog (GWTC-3) of GW sources to constrain the background cosmological parameters together with the astrophysical properties of Binary Black Holes (BBHs), using information from their mass distribution. We expand the public code MGCosmoPop, previously used for the application of this technique, by implementing a state-of-the-art model for the mass distribution, needed to account for the presence of non-trivial features, i.e. a truncated power law with two additional Gaussian peaks, referred to as Multipeak. We then analyse GWTC-3 comparing this model with simpler and more commonly adopted ones, both in the case of fixed and varying cosmology, and assess their goodness-of-fit with different model selection criteria, and their constraining power on the cosmological and population parameters. We also start to explore different sampling methods, namely Markov Chain Monte Carlo and Nested Sampling, comparing their performances and evaluating the advantages of both. We find concurring evidence that the Multipeak model is favoured by the data, in line with previous results, and show that this conclusion is robust to the variation of the cosmological parameters. We find a constraint on the Hubble constant of H0 = 61.10+38.65−22.43 km/s/Mpc (68% C.L.), which shows the potential of this method in providing independent constraints on cosmological parameters. The results obtained in this work have been included in [1].