Low energy elastic and electronically inelastic electron scattering from biomolecules.

Reactions initiated by collisions with low-energy secondary electrons has been found to be the prominent mechanism toward the radiation damage on living tissues through DNA strand breaks. Now it is widely accepted that during the interaction with these secondary species the selective breaking of chemical bonds is triggered by dissociative electron attachment (DEA), that is, the capture of the incident electron and the formation of temporary negative ion states [1,2,3]. One of the approaches largely used toward a deeper understanding of the radiation damage to DNA is through modeling of DEA with its basic constituents (nucleotide bases, sugar and other subunits). We have tried to simplify this approach and attempt to make it comprehensible at a more fundamental level by looking at even simple molecules. Studies involving organic systems such as carboxylic acids, alcohols and simple ¯ve-membered heterocyclic compounds are taken as starting points for these understanding. In the present study we investigate the role played by elastic scattering and electronic excitation of molecules on electron-driven chemical processes. Special attention is focused on the analysis of the in°uence of polarization and multichannel coupling e®ects on the magnitude of elastic and electronically inelastic cross-sections. Our aim is also to investigate the existence of resonances in the elastic and electronically inelastic channels as well as to characterize them with respect to its type (shape, core-excited or Feshbach), symmetry and position. The relevance of these issues is evaluated within the context of possible applications for the modeling of discharge environments and implications in the understanding of mutagenic rupture of DNA chains. The scattering calculations were carried out with the Schwinger multichannel method (SMC) [4] and its implementation with pseudopotentials (SMCPP) [5] at di®erent levels of approximation for impact energies ranging from 0.5 eV to 30 eV. References [1] B. Boudai®a, P. Cloutier, D. Hunting, M. A. Huels and L. Sanche, Science 287, 1658 (2000). [2] X. Pan, P. Cloutier, D. Hunting and L. Sanche, Phys. Rev. Lett. 90, 208102 (2003). [3] F. Martin, P. D. Burrow, Z. Cai, P. Cloutier, D. Hunting and L. Sanche, Phys. Rev. Lett. 93, 068101 (2004). [4] K. Takatsuka and V. McKoy, Phys. Rev. A 24, 2437 (1981); ibid. Phys. Rev. A 30, 1734 (1984). [5] M. H. F. Bettega, L. G. Ferreira and M. A. P. Lima, Phys. Rev. A 47, 1111 (1993).

(18O,16O) Two-neutron transfer reactions for spectroscopic studies.

A systematic study of the response of different nuclei to the (18O, 16O) two-neutron transfer reaction at 84 MeV incident energy was pursued at the INFN-LNS in Catania (Italy). The experiments were performed using several solid targets from light (9Bc, 11 B, 12,13C, 16O, 28Si) to heavier ones (58,64Ni, 120Sn, 208Pb). The 16O ejectiles were detected at forward angles by the MAGNEX magnetic spectrometer and identified without the need of time of flight measurements. Exploiting the large momentum (≈ 25%) and angular (50 msr) acceptance of the spectrometer, energy spectra were obtained with a relevant yield up to about 20 MeV excitation energy. A common feature of the light nuclei spectra is the strong population of states with well known configuration of two-particle over a core and the appearance of unknown resonant structures in the continuum. These latter can reveal the excitation of a collective mode connected with the transfer of a pair. For the heavier nuclei as 66Ni a completely different behaviour is observed indicating the presence of more dissipative processes in the reaction mechanisms that hide the spectroscopic information.

The submarine volcano eruption off El Hierro Island: Effects on the scattering migrant biota and the evolution of the pelagic communities

[EN] The submarine volcano eruption off El Hierro Island (Canary Islands) on 10 October 2011 promoted dramatic perturbation of the water column leading to changes in the distribution of pelagic fauna. To study the response of the scattering biota, we combined acoustic data with hydrographic profiles and concurrent sea surface turbidity indexes from satellite imagery. We also monitored changes in the plankton and nekton communities through the eruptive and post-eruptive phases. Decrease of oxygen, acidification, rising temperature and deposition of chemicals in shallow waters resulted in a reduction of epipelagic stocks and a disruption of diel vertical migration (nocturnal ascent) of mesopelagic organisms. Furthermore, decreased light levels at depth caused by extinction in the volcanic plume resulted in a significant shallowing of the deep acoustic scattering layer. Once the eruption ceased, the distribution and abundances of the pelagic biota returned to baseline levels. There was no evidence of a volcano-induced bloom in the plankton community.

The submarine volcano eruption off El Hierro island: effects on the scattering biota and the evolution of the pelagic community

Trabajo realizado por Ariza, A. V., Kaartvedt, S. Rostad, A. Garijo, J. C., Arístegui, J. Fraile-Nuez, E., Hernández-León, S.

4H silicon carbide particle detectors: study of the defects induced by high energy neutron irradiation

During the last decade advances in the field of sensor design and improved base materials have pushed the radiation hardness of the current silicon detector technology to impressive performance. It should allow operation of the tracking systems of the Large Hadron Collider (LHC) experiments at nominal luminosity (1034 cm-2s-1) for about 10 years. The current silicon detectors are unable to cope with such an environment. Silicon carbide (SiC), which has recently been recognized as potentially radiation hard, is now studied. In this work it was analyzed the effect of high energy neutron irradiation on 4H-SiC particle detectors. Schottky and junction particle detectors were irradiated with 1 MeV neutrons up to fluence of 1016 cm-2. It is well known that the degradation of the detectors with irradiation, independently of the structure used for their realization, is caused by lattice defects, like creation of point-like defect, dopant deactivation and dead layer formation and that a crucial aspect for the understanding of the defect kinetics at a microscopic level is the correct identification of the crystal defects in terms of their electrical activity. In order to clarify the defect kinetic it were carried out a thermal transient spectroscopy (DLTS and PICTS) analysis of different samples irradiated at increasing fluences. The defect evolution was correlated with the transport properties of the irradiated detector, always comparing with the un-irradiated one. The charge collection efficiency degradation of Schottky detectors induced by neutron irradiation was related to the increasing concentration of defects as function of the neutron fluence.

Calculation of the eigenfunctions of the two-group neutron diffusion equation and application to modal decomposition of BWR instabilities

In this thesis, numerical methods aiming at determining the eigenfunctions, their adjoint and the corresponding eigenvalues of the two-group neutron diffusion equations representing any heterogeneous system are investigated. First, the classical power iteration method is modified so that the calculation of modes higher than the fundamental mode is possible. Thereafter, the Explicitly-Restarted Arnoldi method, belonging to the class of Krylov subspace methods, is touched upon. Although the modified power iteration method is a computationally-expensive algorithm, its main advantage is its robustness, i.e. the method always converges to the desired eigenfunctions without any need from the user to set up any parameter in the algorithm. On the other hand, the Arnoldi method, which requires some parameters to be defined by the user, is a very efficient method for calculating eigenfunctions of large sparse system of equations with a minimum computational effort. These methods are thereafter used for off-line analysis of the stability of Boiling Water Reactors. Since several oscillation modes are usually excited (global and regional oscillations) when unstable conditions are encountered, the characterization of the stability of the reactor using for instance the Decay Ratio as a stability indicator might be difficult if the contribution from each of the modes are not separated from each other. Such a modal decomposition is applied to a stability test performed at the Swedish Ringhals-1 unit in September 2002, after the use of the Arnoldi method for pre-calculating the different eigenmodes of the neutron flux throughout the reactor. The modal decomposition clearly demonstrates the excitation of both the global and regional oscillations. Furthermore, such oscillations are found to be intermittent with a time-varying phase shift between the first and second azimuthal modes.

Atmospheric radiative transfer in multiple scattering conditions. Application to NWP models.

High spectral resolution radiative transfer (RT) codes are essential tools in the study of the radiative energy transfer in the Earth atmosphere and a support for the development of parameterizations for fast RT codes used in climate and weather prediction models. Cirrus clouds cover permanently 30% of the Earth's surface, representing an important contribution to the Earth-atmosphere radiation balance. The work has been focussed on the development of the RT model LBLMS. The model, widely tested in the infra-red spectral range, has been extended to the short wave spectrum and it has been used in comparison with airborne and satellite measurements to study the optical properties of cirrus clouds. A new database of single scattering properties has been developed for mid latitude cirrus clouds. Ice clouds are treated as a mixture of ice crystals with various habits. The optical properties of the mixture are tested in comparison to radiometric measurements in selected case studies. Finally, a parameterization of the mixture for application to weather prediction and global circulation models has been developed. The bulk optical properties of ice crystals are parameterized as functions of the effective dimension of measured particle size distributions that are representative of mid latitude cirrus clouds. Tests with the Limited Area Weather Prediction model COSMO have shown the impact of the new parameterization with respect to cirrus cloud optical properties based on ice spheres.

Bestimmung des elektrischen Formfaktors des Neutrons G en in der Reaktion 3 He (e,e'n) bei einem Impulsübertrag Q 2 =0.67 (GeV/c) 2

In der Überprüfung von Kernmodellen stellt der elektrische Formfaktor des Neutrons Gen eine wichtige Eingangsgröße dar. Im Rahmen dieser Arbeit wurde am polarisierten Elektronenstrahl des Mainzer Mikrotrons MAMI der elektrische Formfaktor des Neutrons in der Reaktion 3He(e,e'n) bestimmt. Aus dem Verhältnis der Asymmetrien respektive der Helizitätsumkehr der Elektronenpolarisation mit Targetspin senkrecht und parallel zum Impulsübertrag konnte in einer integralen Analyse Gen bestimmt werden. Zusammengefasst ergab sich Gen(Q2=0.67 (GeV/c)2) = 0.0468 +- 0.0064(stat} +- 0.0027(syst}. Um den Einfluß von Kernbindungseffekten, die z. Zt. von der Gruppe um Prof. Glöckle in Bochum theoretisch gerechnet werden, auch experimentell abzusichern, wurden parallel zur Gen-Messung die Targetasymmetrie Ay0 3He(e,e'n) und Protonenasymmetrie 3He(e,e'p) bestimmt. Die Empfindlichkeit dieser Observablen auf FSI und D-Wellenbeiträge schafft Redundanzen, aus denen auf die Eigenschaften des freien Neutrons im gebundenen System geschlossen werden kann.

Dispersion approach to real and virtual Compton scattering

In recent years, new precision experiments have become possible withthe high luminosity accelerator facilities at MAMIand JLab, supplyingphysicists with precision data sets for different hadronic reactions inthe intermediate energy region, such as pion photo- andelectroproduction and real and virtual Compton scattering.By means of the low energy theorem (LET), the global properties of thenucleon (its mass, charge, and magnetic moment) can be separated fromthe effects of the internal structure of the nucleon, which areeffectively described by polarizabilities. Thepolarizabilities quantify the deformation of the charge andmagnetization densities inside the nucleon in an applied quasistaticelectromagnetic field. The present work is dedicated to develop atool for theextraction of the polarizabilities from these precise Compton data withminimum model dependence, making use of the detailed knowledge of pionphotoproduction by means of dispersion relations (DR). Due to thepresence of t-channel poles, the dispersion integrals for two ofthe six Compton amplitudes diverge. Therefore, we have suggested to subtract the s-channel dispersion integrals at zero photon energy($nu=0$). The subtraction functions at $nu=0$ are calculated through DRin the momentum transfer t at fixed $nu=0$, subtracted at t=0. For this calculation, we use the information about the t-channel process, $gammagammatopipito Nbar{N}$. In this way, four of thepolarizabilities can be predicted using the unsubtracted DR in the $s$-channel. The other two, $alpha-beta$ and $gamma_pi$, are free parameters in ourformalism and can be obtained from a fit to the Compton data.We present the results for unpolarized and polarized RCS observables,%in the kinematics of the most recent experiments, and indicate anenhanced sensitivity to the nucleon polarizabilities in theenergy range between pion production threshold and the $Delta(1232)$-resonance.newlineindentFurthermore,we extend the DR formalism to virtual Compton scattering (radiativeelectron scattering off the nucleon), in which the concept of thepolarizabilities is generalized to the case of avirtual initial photon by introducing six generalizedpolarizabilities (GPs). Our formalism provides predictions for the fourspin GPs, while the two scalar GPs $alpha(Q^2)$ and $beta(Q^2)$ have to befitted to the experimental data at each value of $Q^2$.We show that at energies betweenpion threshold and the $Delta(1232)$-resonance position, thesensitivity to the GPs can be increased significantly, as compared tolow energies, where the LEX is applicable. Our DR formalism can be used for analysing VCS experiments over a widerange of energy and virtuality $Q^2$, which allows one to extract theGPs from VCS data in different kinematics with a minimum of model dependence.

Confocal microscopy applied to the study of single entity fluorescence and light scattering

A sample scanning confocal optical microscope (SCOM) was designed and constructed in order to perform local measurements of fluorescence, light scattering and Raman scattering. This instrument allows to measure time resolved fluorescence, Raman scattering and light scattering from the same diffraction limited spot. Fluorescence from single molecules and light scattering from metallic nanoparticles can be studied. First, the electric field distribution in the focus of the SCOM was modelled. This enables the design of illumination modes for different purposes, such as the determination of the three-dimensional orientation of single chromophores. Second, a method for the calculation of the de-excitation rates of a chromophore was presented. This permits to compare different detection schemes and experimental geometries in order to optimize the collection of fluorescence photons. Both methods were combined to calculate the SCOM fluorescence signal of a chromophore in a general layered system. The fluorescence excitation and emission of single molecules through a thin gold film was investigated experimentally and modelled. It was demonstrated that, due to the mediation of surface plasmons, single molecule fluorescence near a thin gold film can be excited and detected with an epi-illumination scheme through the film. Single molecule fluorescence as close as 15nm to the gold film was studied in this manner. The fluorescence dynamics (fluorescence blinking and excited state lifetime) of single molecules was studied in the presence and in the absence of a nearby gold film in order to investigate the influence of the metal on the electronic transition rates. The trace-histogram and the autocorrelation methods for the analysis of single molecule fluorescence blinking were presented and compared via the analysis of Monte-Carlo simulated data. The nearby gold influences the total decay rate in agreement to theory. The gold presence produced no influence on the ISC rate from the excited state to the triplet but increased by a factor of 2 the transition rate from the triplet to the singlet ground state. The photoluminescence blinking of Zn0.42Cd0.58Se QDs on glass and ITO substrates was investigated experimentally as a function of the excitation power (P) and modelled via Monte-Carlo simulations. At low P, it was observed that the probability of a certain on- or off-time follows a negative power-law with exponent near to 1.6. As P increased, the on-time fraction reduced on both substrates whereas the off-times did not change. A weak residual memory effect between consecutive on-times and consecutive off-times was observed but not between an on-time and the adjacent off-time. All of this suggests the presence of two independent mechanisms governing the lifetimes of the on- and off-states. The simulated data showed Poisson-distributed off- and on-intensities, demonstrating that the observed non-Poissonian on-intensity distribution of the QDs is not a product of the underlying power-law probability and that the blinking of QDs occurs between a non-emitting off-state and a distribution of emitting on-states with different intensities. All the experimentally observed photo-induced effects could be accounted for by introducing a characteristic lifetime tPI of the on-state in the simulations. The QDs on glass presented a tPI proportional to P-1 suggesting the presence of a one-photon process. Light scattering images and spectra of colloidal and C-shaped gold nano-particles were acquired. The minimum size of a metallic scatterer detectable with the SCOM lies around 20 nm.

Strahlungsbegleitete Pion-Photoproduktion am Proton

Im Jahre 1997 wurden von Tatischeff et al. bei der Reaktion p p -> X p pi+ resonanzartige Zustände im Spektrum der invarianten Masse des fehlenden Nukleons X bei M = 1004, 1044 und 1094 MeV gefunden. In einem zweiten Experiment von Filkov et al. beobachtete man bei der Reaktion p d -> p p X Resonanzstrukturen bei M = 966, 986 und 1003 MeV. Solche exotischen Resonanzen widersprechen etablierten Nukleonenmodellen, die die Delta(1232)-Resonanz als ersten Anregungszustand beschreiben. Zur Deutung der beobachteten Strukturen wurden Quarkcluster-Modelle mit und ohne Farb-Magnet-Wechselwirkungen entwickelt. Lvov et al. zweifelten die experimentellen Ergebnisse an, da keine Strukturen in den Daten zur reellen Comptonstreuung gefunden wurden. Als Gegenargument wurde von Kobushkin vorgeschlagen, dass diese Resonanzen eine total-antisymmetrische Spin-Flavour-Wellenfunktion haben und nur der N-2Gamma-Zerfall erlaubt wäre. In dieser Arbeit wurde die Reaktion g p -> X pi+ -> n g g pi+ zur Suche nach diesen exotischen Resonanzen verwendet. Die Daten wurden parallel zur Messung der Pion-Polarisierbarkeiten am Mainzer Beschleuniger MAMI genommen. Durch Bremsstrahlung der Elektronen an einer Radiatorfolie wurden reelle Photonen erzeugt, deren Energie von der A2-Photonenmarkierungsanlage (Glasgow-Tagger) bestimmt wurde. Als Protontarget wurde ein 10 cm langes Flüssigwasserstoff-Target verwendet. Geladene Reaktionsprodukte wurden unter Vorwärtswinkeln Theta < 20 Grad bezüglich der Strahlachse in einer Vieldraht-Proportionalkammer nachgewiesen, während Photonen im Spektrometer TAPS mit 526 BaF2-Kristallen unter Polarwinkeln Theta > 60 Grad detektiert wurden. Zum Nachweis von Neutronen stand ein Flugzeitdetektor mit insgesamt 111 Einzelmodulen zur Verfügung. Zum Test der Analysesoftware und des experimentellen Aufbaus wurden zusätzlich die Reaktionskanäle g p -> p pi0 und g p -> n pi0 pi+ ausgewertet. Für die Ein-Pion-Produktion wurden differentielle Wirkungsquerschnitte unter Rückwärtswinkeln bestimmt und mit theoretischen Modellen und experimentellen Werten verglichen. Für den Kanal g p -> n pi0 pi+ wurden Spektren invarianter Massen für verschiedene Teilchenkombinationen ermittelt und mit einer Simulation verglichen. Die Daten legen nahe, dass die Reaktion hauptsächlich über eine Anregung der Delta0(1232)-Resonanz verläuft. Bei der Suche nach exotischen Resonanzen wurden keine statistisch signifikanten Strukturen gefunden. Es wurden Obergrenzen für den differentiellen Wirkungsquerschnitt ermittelt.

108 Sn studied with intermediate-energy Coulomb excitation

In dieser Arbeit wurde der instabile, Neutronenarme Kern 108Sn mit Hilfe der Coulomb-Anregung bei intermediaeren Energien in inverser Kinematik studiert. Diese Methode wurde bisher zur Untersuchung der ersten angeregten 2+ Zustaende und deren E2 Zerfallsraten in Kernen mit Kernladungszahl Z< 30 angewendet. 108Sn ist somit der Kern mit der groeßten Kernladungszahl, bei dem diese Studien bisher stattfanden. Das Ziel dieses Experiments war die Messung der unbekannten reduzierten Uebergangswahrscheinlichkeit B(E2,0+ -> 2+). Der B(E2)-Wert von 0.230(57) e2b2 wurde relativ zu dem bekannten Wert des Isotops 112Sn bestimmt. Das Experiment wurde an der GSI Darmstadt mit Hilfe des RISING Detektors und des Fragmentseperators (FRS) durchgefuehrt. Sekundaere Strahlen (108Sn, 112Sn) mit einer Energie von ca. 150 MeV pro Nukleon wurden auf ein 386 mg/cm2 dickes 197Au Target geschossen. Die Projektilfragmente wurden mit Hilfe des Fragmentseparators selektiert und identifiziert. Zur Selektion des Reaktionskanals und zur Bestimmung des Winkels der gestreuten Fragmente wurde das Teilchenteleskop CATE, das sich hinter dem Target befand, verwendet. Gammastrahlung, die in Koinzidenz mit den Projektilrestkernen emittiert wurde, wurde in den Germanium-Cluster Detektoren des RISING Detektors nachgewiesen. Der gemessene B(E2,0+ -> 2+)-Wert von 108Sn ist in Uebereinstimmung mit neueren Schalenmodellrechnungen, die auf realistischen effektiven Wechselwirkungen basieren und im Rahmen eines verallgemeinerten Seniorit¨ats-Schemas erklaert werden.

Study of electromagnetic reactions on light nuclei with the Lorentz Integral Transform method

In dieser Arbeit aus dem Bereich der Wenig-Nukleonen-Physik wird die neu entwickelte Methode der Lorentz Integral Transformation (LIT) auf die Untersuchung von Kernphotoabsorption und Elektronenstreuung an leichten Kernen angewendet. Die LIT-Methode ermoeglicht exakte Rechnungen durchzufuehren, ohne explizite Bestimmung der Endzustaende im Kontinuum. Das Problem wird auf die Loesung einer bindungzustandsaehnlichen Gleichung reduziert, bei der die Endzustandswechselwirkung vollstaendig beruecksichtigt wird. Die Loesung der LIT-Gleichung wird mit Hilfe einer Entwicklung nach hypersphaerischen harmonischen Funktionen durchgefuehrt, deren Konvergenz durch Anwendung einer effektiven Wechselwirkung im Rahmem des hypersphaerischen Formalismus (EIHH) beschleunigt wird. In dieser Arbeit wird die erste mikroskopische Berechnung des totalen Wirkungsquerschnittes fuer Photoabsorption unterhalb der Pionproduktionsschwelle an 6Li, 6He und 7Li vorgestellt. Die Rechnungen werden mit zentralen semirealistischen NN-Wechselwirkungen durchgefuehrt, die die Tensor Kraft teilweise simulieren, da die Bindungsenergien von Deuteron und von Drei-Teilchen-Kernen richtig reproduziert werden. Der Wirkungsquerschnitt fur Photoabsorption an 6Li zeigt nur eine Dipol-Riesenresonanz, waehrend 6He zwei unterschiedliche Piks aufweist, die dem Aufbruch vom Halo und vom Alpha-Core entsprechen. Der Vergleich mit experimentellen Daten zeigt, dass die Addition einer P-Wellen-Wechselwirkung die Uebereinstimmung wesentlich verbessert. Bei 7Li wird nur eine Dipol-Riesenresonanz gefunden, die gut mit den verfuegbaren experimentellen Daten uebereinstimmt. Bezueglich der Elektronenstreuung wird die Berechnung der longitudinalen und transversalen Antwortfunktionen von 4He im quasi-elastischen Bereich fuer mittlere Werte des Impulsuebertrages dargestellt. Fuer die Ladungs- und Stromoperatoren wird ein nichtrelativistisches Modell verwendet. Die Rechnungen sind mit semirealistischen Wechselwirkungen durchgefuert und ein eichinvarianter Strom wird durch die Einfuehrung eines Mesonaustauschstroms gewonnen. Die Wirkung des Zweiteilchenstroms auf die transversalen Antwortfunktionen wird untersucht. Vorlaeufige Ergebnisse werden gezeigt und mit den verfuegbaren experimentellen Daten verglichen.

Monte Carlo simulation of the WENDI-2 neutron dosimeter

Il presente lavoro di tesi, sviluppato nell’arco di sei mesi presso l’Institut Supérieur Industriel de Bruxelles (ISIB) in collaborazione con Ion Beam Application Group (IBA, Louvain la Neuve), ha come principale soggetto lo studio della risposta del rem meter WENDI-2 commercializzato da Thermo Scientific. Lo studio si è basato principalmente sull’uso del codice Monte Carlo MCNPX 2.5.0, simulando la risposta del detector sia in caso di campi di radiazione neutronica monoenergetici sia in corrispondenza di spettri neutronici continui. La prima fase è stata dedicata alla modellizzazione MCNPX del rem counter, consentendo così la valutazione della sua funzione risposta. Questa è stata ricostruita interpolando 93 punti, ciascuno calcolato in corrispondenza di un singolo valore di energia di una sorgente puntiforme, compreso tra 1 meV e 5 GeV. In tal caso è stata rilevata un’ottima corrispondenza tra i risultati ottenuti e quelli riportati nella letteratura scientifica esistente. In una seconda fase, al fine di ottenere informazioni sulla risposta di WENDI II in corrispondenza di campi complessi di radiazione, simulazioni MCNPX sono state realizzate riproducendo un ambiente di lavoro esistente presso la sede IBA di Louvain la Neuve: la risposta del detector è stata valutata in corrispondenza di 9 diverse posizioni all’interno di un bunker contenente un ciclotrone PET (18 MeV H-), implicando la rilevazione di campi di radiazione neutronica continui ed estesi dalle energie termiche fino a 18 MeV. I risultati ottenuti sono stati infine comparati con i valori di dose ambiente equivalente calcolata nelle stesse condizioni di irraggiamento.

Measurement of the muon-induced neutron flux at LNGS with the LVD experiment

In this thesis we describe in detail the Monte Carlo simulation (LVDG4) built to interpret the experimental data collected by LVD and to measure the muon-induced neutron yield in iron and liquid scintillator. A full Monte Carlo simulation, based on the Geant4 (v 9.3) toolkit, has been developed and validation tests have been performed. We used the LVDG4 to determine the active vetoing and the shielding power of LVD. The idea was to evaluate the feasibility to host a dark matter detector in the most internal part, called Core Facility (LVD-CF). The first conclusion is that LVD is a good moderator, but the iron supporting structure produce a great number of neutrons near the core. The second conclusions is that if LVD is used as an active veto for muons, the neutron flux in the LVD-CF is reduced by a factor 50, of the same order of magnitude of the neutron flux in the deepest laboratory of the world, Sudbury. Finally, the muon-induced neutron yield has been measured. In liquid scintillator we found $(3.2 \pm 0.2) \times 10^{-4}$ n/g/cm$^2$, in agreement with previous measurements performed at different depths and with the general trend predicted by theoretical calculations and Monte Carlo simulations. Moreover we present the first measurement, in our knowledge, of the neutron yield in iron: $(1.9 \pm 0.1) \times 10^{-3}$ n/g/cm$^2$. That measurement provides an important check for the MC of neutron production in heavy materials that are often used as shield in low background experiments.