Search NASA⌕ Search

SEARCH · Search NASA

Results for “Image coding”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,351 records · Page 75

Practical jet flavour through NNLO

Abstract An infrared and collinear (IRC) safe definition of the partonic flavour of a jet is vital for precision predictions of quantum chromodynamics at colliders. Jet flavour definitions have been presented in the literature, but they are typically defined through modification of the jet algorithm to be sensitive to partonic flavour at every stage of the clustering. While this does ensure that the sum of flavours in a jet is IRC safe, a flavour-sensitive clustering procedure is difficult to apply to realistic data. We introduce a distinct and novel approach to jet flavour that can be applied to a collection of partons defined by any algorithm. Our definition of jet flavour is the sum of flavours of all partons that remain after Soft Drop grooming, reclustered with the "Image missing" <#comment/> algorithm. We prove that this prescription is IRC safe through next-to-next-to-leading order (NNLO), and so can interface with the most precise fixed-order calculations for jets available at present. We validate the IRC safety of this definition with numeric fixed-order codes and further show that jet flavour with Soft Drop reclustered with a generalised $$k_T$$ k T algorithm fails to be IRC safe at NNLO.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Shadow Knows: Using Shadows to Investigate the Structure of the Pretransitional Disk of HD 100453

We present Gemini Planet Imager polarized intensity imagery of HD 100453 in Y , J , and K 1 bands that reveals an inner gap (9–18 au), an outer disk (18–39 au) with two prominent spiral arms, and two azimuthally localized dark features that are also present in Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) total intensity images. Spectral energy distribution fitting further suggests that the radial gap extends to 1 au. The narrow, wedge-like shape of the dark features appears similar to predictions of shadows cast by an inner disk that is misaligned with respect to the outer disk. Using the Monte Carlo radiative transfer code HOCHUNCK3D, we construct a model of the disk that allows us to determine its physical properties in more detail. From the angular separation of the features, we measure the difference in inclination between the disks (45°) and their major axes, PA = 140° east of north for the outer disk, and 100° for the inner disk. We find an outer-disk inclination of 25° ± 10° from face-on, in broad agreement with the Wagner et al. measurement of 34°. SPHERE data in J and H bands indicate a reddish disk, which indicates that HD 100453 is evolving into a young debris disk.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of an Enhanced Radiation Physics Toolset for Modeling Spectral and Imaging Signatures in the Warm Dense Matter Experiments

Radiative and atomic processes in plasmas play a critical role in a wide variety of high energy density laboratory plasma (HEDLP) experiments. The emission, absorption, and transport of radiation can strongly affect the overall energetics and evolution of such plasmas. In addition, radiation-based diagnostics – including imaging, spectroscopy, and absolute flux measurements – are widely used to determine key features of HEDLPs. To advance our understanding of HEDLP science, it is vital to have high-fidelity computational physics tools that have well-tested radiation physics modeling, and that are readily accessible to researchers in the HEDLP community. Simulations play an extremely important role for planning and designing the experiments, as well as for post-experiment data analysis. Prism Computational Sciences develops software that is used by National Laboratories and universities (including five members of LaserNetUS network). Prominent examples of such research efforts include z-pinch and short-pulse laser experiments designed to study the basic physics of photoionized plasmas and photoionization fronts, as well as their application to astrophysical plasmas. The main effort was dedicated to the development of non-equilibrium equation-of-state (EOS) models within the HELIOS-CR code, a hydrodynamics code with inline collisional-radiative atomic kinetics. Gas cell experiments on Z and Omega demonstrated the importance of non-equilibrium effects on atomic kinetics in photoionized plasmas. Recent proof-of-principle experiments on Omega EP confirmed the advantages of using a short-pulse laser to create an intense radiation drive, leading to additional experiments being proposed. Photoionization front experiments at LLE also emphasizes the importance of radiation and atomic physics. In both studies, HELIOS-CR simulations played a crucial role in computing non-equilibrium opacities and ionization distributions. A newly developed non-LTE EOS model will help addressing possible non-equilibrium effects, on for example specific heat, and their influence in plasma evolution. Prism also implemented support for open-source atomic data generated by the Flexible Atomic Code. This allows researchers to generate custom atomic tables and use them within the complex framework of simulation tools developed by Prism. The ability to use open-source atomic data would be extremely valuable for hydrodynamics and spectroscopic simulations that include high-Z materials, e.g., picosecond x-ray pulse generation experiments. Support for new atomic structures was fully implemented, and the data can be used by all simulation tools developed at Prism: radiation-hydrodynamics, imaging and spectroscopy, EOS and opacity. The development resulted in a significant fidelity enhancement to the simulations tools developed by Prism that are currently used in other cutting-edge experiments including: opacity measurement experiments performed to both understand the basic radiative and atomic properties of plasmas as well as provide data for more accurately modeling the internal structure of the Sun and other stars, high-intensity short-pulse laser experiments performed to develop short-wavelength light sources for use as backlighters and to investigate fast ignition concepts for inertial fusion energy; capsule implosion experiments designed to develop inertial fusion as an energy source, etc.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental measurements of fluid flow in an 84-pin hexagonal rod bundle with spacer grid for a gas-cooled fast modular reactor

A 50 MW e helium-cooled fast modular reactor (FMR) is under development, and as part of the Department of Energy Integrated Research Project (IRP), Texas A&M University is conducting the thermal-hydraulic characterization of its baseline core configuration. Here, we experimentally investigated the axial and cross flow fields characteristics in the hexagonal fuel rod bundle composed of 84 rods, a central rod, and spacer grids at a Reynolds number of 12,000. A fully transparent experimental facility resembling of one unit of the fuel assembly was constructed. Time-resolved particle image velocimetry (TR-PIV) measurements were performed to characterize hydraulic behavior downstream the spacer grid. Velocity measurements were conducted in the axial and radial direction of the rod bundle. From the PIV velocity vector fields, the full-field flow statistics were computed for the mean velocity, vorticity, and Reynolds stresses. An analysis of the energy decay downstream of the spacer grid was conducted with calculations of secondary-flow intensities, turbulent kinetic energy, power spectrum, and spatial-temporal velocity cross correlations. The vorticity field was obtained and pairs of counter-rotating vortexes were identified using a 3D reconstruction of several measurement planes. The data from this experimental campaign will be used to validate numerical models based on Computational Fluid Dynamics and other codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

'Edge illumination' in X-ray Phase Contrast Imaging

In the late '90s, the concept of 'edge illumination' was developed at ELETTRA in Italy as an alternative method to increase the phase sensitivity of an imaging system. The main idea was to be able to reproduce the fine angular selection of 'analyzer' crystals without actually using a crystal, as this would allow employing the method with divergent and polychromatic (i.e. conventional) x-ray sources. It was observed that this could be achieved by illuminating only the edges of the detector pixels, and that the method's sensitivity could be progressively increased by illuminating smaller pixel fractions closer to its physical edge. A few years later the idea was adapted for use with a conventional source by means of two sets of x-ray masks ('coded aperture' masks), which enabled obtaining the same effect for each row (or column) of pixels of an area detector illuminated by a cone beam. This article reviews the method and presents recent examples of application.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Experimental observations of detached bow shock formation in the interaction of a laser-produced plasma with a magnetized obstacle

The magnetic field produced by planets with active dynamos, like the Earth, can exert sufficient pressure to oppose supersonic stellar wind plasmas, leading to the formation of a standing bow shock upstream of the magnetopause, or pressure-balance surface. Scaled laboratory experiments studying the interaction of an inflowing solar wind analog with a strong, external magnetic field are a promising new way to study magnetospheric physics and to complement existing models, although reaching regimes favorable for magnetized shock formation is experimentally challenging. This paper presents experimental evidence of the formation of a magnetized bow shock in the interaction of a supersonic, super-Alfvenic plasma with a strongly magnetized obstacle at the OMEGA laser facility. The solar wind analog is generated by the collision and subsequent expansion of two counterpropagating, laser-driven plasma plumes. The magnetized obstacle is a thin wire, driven with strong electrical currents. Hydrodynamic simulations using the FLASH code predict that the colliding plasma source meets the criteria for bow shock formation. Spatially resolved, optical Thomson scattering measures the electron number density, and optical emission lines provide a measurement of the plasma temperature, from which we infer the presence of a fast magnetosonic shock far upstream of the obstacle. Proton images provide a measure of large-scale features in the magnetic field topology, and reconstructed path-integrated magnetic field maps from these images suggest the formation of a bow shock upstream of the wire and as a transient magnetopause. Finally, we compare features in the reconstructed fields to two-dimensional MHD simulations of the system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Results from a synthetic model of the ITER XRCS-Core diagnostic based on high-fidelity x-ray ray tracing

A high-fidelity synthetic diagnostic has been developed for the ITER core x-ray crystal spectrometer diagnostic based on x-ray ray tracing. This synthetic diagnostic has been used to model expected performance of the diagnostic, to aid in diagnostic design, and to develop engineering tolerances. The synthetic model is based on x-ray ray tracing using the recently developed xicsrt ray tracing code and includes a fully three-dimensional representation of the diagnostic based on the computer aided design. The modeled components are: plasma geometry and emission profiles, highly oriented pyrolytic graphite pre-reflectors, spherically bent crystals, and pixelated x-ray detectors. Plasma emission profiles have been calculated for Xe 44+ , Xe 47+ , and Xe 51+ , based on an ITER operational scenario available through the Integrated Modelling & Analysis Suite database, and modeled within the ray tracing code as a volumetric x-ray source; the shape of the plasma source is determined by equilibrium geometry and an appropriate wavelength distribution to match the expected ion temperature profile. All individual components of the x-ray optical system have been modeled with high-fidelity producing a synthetic detector image that is expected to closely match what will be seen in the final as-built system. Particular care is taken to maintain preservation of photon statistics throughout the ray tracing allowing for quantitative estimates of diagnostic performance.

47 OTHER INSTRUMENTATION↗

MODELING MAGNETIC FIELD STRUCTURE OF A SOLAR ACTIVE REGION CORONA USING NONLINEAR FORCE-FREE FIELDS IN SPHERICAL GEOMETRY

We test a nonlinear force-free field (NLFFF) optimization code in spherical geometry using an analytical solution from Low and Lou. Several tests are run, ranging from idealized cases where exact vector field data are provided on all boundaries, to cases where noisy vector data are provided on only the lower boundary (approximating the solar problem). Analytical tests also show that the NLFFF code in the spherical geometry performs better than that in the Cartesian one when the field of view of the bottom boundary is large, say, 20 Degree-Sign Multiplication-Sign 20 Degree-Sign . Additionally, we apply the NLFFF model to an active region observed by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory (SDO) both before and after an M8.7 flare. For each observation time, we initialize the models using potential field source surface (PFSS) extrapolations based on either a synoptic chart or a flux-dispersal model, and compare the resulting NLFFF models. The results show that NLFFF extrapolations using the flux-dispersal model as the boundary condition have slightly lower, therefore better, force-free, and divergence-free metrics, and contain larger free magnetic energy. By comparing the extrapolated magnetic field lines with the extreme ultraviolet (EUV) observations by the Atmospheric Imaging Assembly on board SDO, we find that the NLFFF performs better than the PFSS not only for the core field of the flare productive region, but also for large EUV loops higher than 50 Mm.

79 ASTRONOMY AND ASTROPHYSICS↗

Global energetics of solar flares. I. Magnetic energies

We present the first part of a project on the global energetics of solar flares and coronal mass ejections that includes about 400 M- and X-class flares observed with Atmospheric Imaging Assembly (AIA) and Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). We calculate the potential (E{sub p} ), the nonpotential (E {sub np}) or free energies (E {sub free} = E {sub np} – E{sub p} ), and the flare-dissipated magnetic energies (E {sub diss}). We calculate these magnetic parameters using two different NLFFF codes: the COR-NLFFF code uses the line-of-sight magnetic field component B{sub z} from HMI to define the potential field, and the two-dimensional (2D) coordinates of automatically detected coronal loops in six coronal wavelengths from AIA to measure the helical twist of coronal loops caused by vertical currents, while the PHOT-NLFFF code extrapolates the photospheric three-dimensional (3D) vector fields. We find agreement between the two codes in the measurement of free energies and dissipated energies within a factor of ≲ 3. The size distributions of magnetic parameters exhibit powerlaw slopes that are approximately consistent with the fractal-diffusive self-organized criticality model. The magnetic parameters exhibit scaling laws for the nonpotential energy, E{sub np}∝E{sub p}{sup 1.02}, for the free energy, E{sub free}∝E{sub p}{sup 1.7} and E{sub free}∝B{sub φ}{sup 1.0}L{sup 1.5}, for the dissipated energy, E{sub diss}∝E{sub p}{sup 1.6} and E{sub diss}∝E{sub free}{sup 0.9}, and the energy dissipation volume, V∝E{sub diss}{sup 1.2}. The potential energies vary in the range of E{sub p} = 1 × 10{sup 31}-4 × 10{sup 33} erg, while the free energy has a ratio of E {sub free}/E{sub p} ≈ 1%-25%. The Poynting flux amounts to F {sub flare} ≈ 5 × 10{sup 8}-10{sup 10} erg cm{sup –2} s{sup –1} during flares, which averages to F {sub AR} ≈ 6 × 10{sup 6} erg cm{sup –2} s{sup –1} during the entire observation period and is comparable with the coronal heating rate requirement in active regions.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-J-197: Investigation of Microsoft Kinect 2.0 Depth Resolution for Patient Motion Tracking

Purpose: Investigate the use of the Kinect 2.0 for patient motion tracking during radiotherapy by studying spatial and depth resolution capabilities. Methods: Using code written in C#, depth map data was abstracted from the Kinect to create an initial depth map template indicative of the initial position of an object to be compared to the depth map of the object over time. To test this process, simple setup was created in which two objects were imaged: a 40 cm × 40 cm board covered in non reflective material and a 15 cm × 26 cm textbook with a slightly reflective, glossy cover. Each object, imaged and measured separately, was placed on a movable platform with object to camera distance measured. The object was then moved a specified amount to ascertain whether the Kinect’s depth camera would visualize the difference in position of the object. Results: Initial investigations have shown the Kinect depth resolution is dependent on the object to camera distance. Measurements indicate that movements as small as 1 mm can be visualized for objects as close as 50 cm away. This depth resolution decreases linearly with object to camera distance. At 4 m, the depth resolution had decreased to observe a minimum movement of 1 cm. Conclusion: The improved resolution and advanced hardware of the Kinect 2.0 allows for increase of depth resolution over the Kinect 1.0. Although obvious that the depth resolution should decrease with increasing distance from an object given the decrease in number of pixels representing said object, the depth resolution at large distances indicates its usefulness in a clinical setting.

60 APPLIED LIFE SCIENCES↗

A miniaturized long-life acoustic transmitter for animal tracking in marine environments

Marine and hydrokinetic (MHK) energy devices are early-stage technologies and their impacts to ocean life and environment are largely unknown, which consequently limits the progress of deploying these devices at commercial scale. Optical and acoustic imaging have been the primary technologies for monitoring marine animals’ behavior at existing or potential MHK sites. However, these technologies lack the ability to identify individual animals, which makes them less useful in identifying protected species as well as studying the animals’ long-term behavior. Implantable acoustic transmitters with unique identification codes could be used as a complementary technology. To minimize the tag burden on the tagged animal (and hence the resulting bias) and to satisfy the long-term monitoring needs, it is desired that the implanted transmitter be as small and light as possible and can last for as long as possible. To date, the smallest acoustic transmitters suitable for marine animal tracking have relatively short service life and detection range (up to 200 m). In this report, we introduce a new miniaturized acoustic transmitter which operates at 208.3 kHz. It offers a smaller dry weight (0.45 g), a significantly longer service life (140 days at a 5-second ping rate) as well as a greater detection range than the existing commercial counterparts, which will help obtain more comprehensive movement and behavioral information of animals around MHK energy devices.

Li, Huidong↗

MO-AB-206-00: Nuclear Medicine Physics and Testing

This education session will cover the physics and operation principles of gamma cameras and PET scanners. The first talk will focus on PET imaging. An overview of the principles of PET imaging will be provided, including positron decay physics, and the transition from 2D to 3D imaging. More recent advances in hardware and software will be discussed, such as time-of-flight imaging, and improvements in reconstruction algorithms that provide for options such as depth-of-interaction corrections. Quantitative applications of PET will be discussed, as well as the requirements for doing accurate quantitation. Relevant performance tests will also be described. Learning Objectives: Be able to describe basic physics principles of PET and operation of PET scanners. Learn about recent advances in PET scanner hardware technology. Be able to describe advances in reconstruction techniques and improvements Be able to list relevant performance tests. The second talk will focus on gamma cameras. The Nuclear Medicine subcommittee has charged a task group (TG177) to develop a report on the current state of physics testing of gamma cameras, SPECT, and SPECT/CT systems. The report makes recommendations for performance tests to be done for routine quality assurance, annual physics testing, and acceptance tests, and identifies those needed satisfy the ACR accreditation program and The Joint Commission imaging standards. The report is also intended to be used as a manual with detailed instructions on how to perform tests under widely varying conditions. Learning Objectives: At the end of the presentation members of the audience will: Be familiar with the tests recommended for routine quality assurance, annual physics testing, and acceptance tests of gamma cameras for planar imaging. Be familiar with the tests recommended for routine quality assurance, annual physics testing, and acceptance tests of SPECT systems. Be familiar with the tests of a SPECT/CT system that include the CT images for SPECT reconstructions. Become knowledgeable of items to be included in annual acceptance testing reports including CT dosimetry and PACS monitor measurements. T. Turkington, GE Healthcare.

61 RADIATION PROTECTION AND DOSIMETRY↗

MPF X-Ray Module Setup for 2018

The Mobile Plutonium Facility (MPF) is a portable module design facility that can be deployed anywhere in the world for the safe recovery, packaging, and shipment of a nations legacy, or excess plutonium stockpile. All of the modules of MPF have specific purposes for operating in a mission if one were to ever come up. For my past four internships here at SRNL I have been working on the support of the X-Ray module for MPF. This year I installed new equipment I had purchased last year for the X-Ray Module and made sure the module was mission ready for the chance that the team deploys in the future for plutonium recovery. Using digital radiography (DR) and computed tomography (CT) provides the MPF team the ability to verify the contents of containers that would be holding material. This also ensures the safety of the team members in the glovebox facilities by verifying that there are no sharps, liquids, or other items that may harm the workers or glovebox operations when working with radioactive material. X-Ray imaging can also pick up low density objects like plastic bags and tape. Note how the DR image hides the puncture in the inner can while the CT image shows the damage. New DR and CT computers with upgraded software, New imaging panel, NOVO Portable X-Ray System. Because of the new upgrades to the X-Ray module system containers can be scanned and evaluated quicker than before. With the last system the 3D model reconstruction alone would take anywhere from 10-12 minutes, now with the new system that time has been brought down to just shy of 1 minute 30 seconds. The new portable system also can 1,2 or 4 images of an object and stitch them together into a mosaic for scanning of larger containers.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Poster - 15: Air kerma rate measurements using a cylindrical free-air ionization chamber on the 05B1-1 beamline at the Canadian Light Source

Purpose: The BioMedical Imaging and Therapy (BMIT) beamlines at the Canadian Light Source offer the opportunity for investigating novel imaging and therapy applications of synchrotron radiation. A necessary component in advancing this research, and in progressing toward clinical applications, is the availability of accurate dosimetry that is traceable to a standards institution. However, dosimetry in this setting is challenging. These beams are typically small, non-uniform, and highly intense. This work describes air kerma rate measurements on a BMIT beamline using a free-air ionization chamber (FAC). Methods: The measurements were taken at the 05B1-1 beamline (∼8 – 100 keV) for several beam qualities with mean energies between 20.0 and 84.0 keV. The Victoreen Model 480 cylindrical FAC, with a specially fabricated 0.52 mm diameter aperture, was used to measure air kerma rates. The required correction factors were determined using a variety of methods: tabulated data, measurements, theoretical calculations and Monte Carlo simulations (EGSnrc user code egs-fac). Results: The experimental air kerma rates measured between 0.270 ± 13.6% and 312 ± 2.7% Gy/min. At lower energies (low filtration), the most impactful correction factors were those for ion recombination and for x-ray attenuation. Conclusions: These measurements marked the first absolute dosimetry performed at the BMIT beamlines. The experimental and Monte Carlo methods developed will allow air kerma rates to be measured under other experimental conditions, provide a benchmark to which other dosimeters will be compared, and provide a reference for imaging and therapy research programs on this beamline.

60 APPLIED LIFE SCIENCES↗

Investigation of radiative bow-shocks in magnetically accelerated plasma flows

We present a study of the formation of bow shocks in radiatively cooled plasma flows. This work uses an inverse wire array to provide a quasi-uniform, large scale hydrodynamic flow accelerated by Lorentz forces to supersonic velocities. This flow impacts a stationary object placed in its path, forming a well-defined Mach cone. Interferogram data are used to determine a Mach number of ∼6, which may increase with radial position suggesting a strongly cooling flow. Self-emission imaging shows the formation of a thin (<60 μm) strongly emitting shock region, where T{sub e} ∼ 40–50 eV, and rapid cooling behind the shock. Emission is observed upstream of the shock position which appears consistent with a radiation driven phenomenon. Data are compared to 2-dimensional simulations using the Gorgon MHD code, which show good agreement with the experiments. The simulations are also used to investigate the effect of magnetic field in the target, demonstrating that the bow-shocks have a high plasma β, and the influence of B-field at the shock is small. This consistent with experimental measurement with micro bdot probes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cathodoluminescence spectrum imaging analysis of CdTe thin-film bevels

We conducted T = 6 K cathodoluminescence (CL) spectrum imaging with a nanoscale electron beam on beveled surfaces of CdTe thin films at the critical stages of standard CdTe solar cell fabrication. We find that the through-thickness CL total intensity profiles are consistent with a reduction in grain-boundary recombination due to the CdCl 2 treatment. The color-coded CL maps of the near-band-edge transitions indicate significant variations in the defect recombination activity at the micron and sub-micron scales within grains, from grain to grain, throughout the film depth, and between films with different processing histories. We estimated the grain-interior sulfur-alloying fraction in the interdiffused CdTe/CdS region of the CdCl 2 -treated films from a sample of 35 grains and found that it is not strongly correlated with CL intensity. A kinetic rate-equation model was used to simulate grain-boundary (GB) and grain-interior CL spectra. Simulations indicate that the large reduction in the exciton band intensity and relatively small decrease in the lower-energy band intensity at CdTe GBs or dislocations can be explained by an enhanced electron-hole non-radiative recombination rate at the deep GB or dislocation defects. Additionally, simulations also show that higher GB concentrations of donors and/or acceptors can increase the lower-energy band intensity, while slightly decreasing the exciton band intensity.

14 SOLAR ENERGY↗

Study of the spatial resolution of laser thermochemical technology for recording diffraction microstructures

The thermochemical method for recording data, which is based on local laser oxidation of a thin metal film with subsequent etching of the unirradiated region, is an alternative to laser photolithography and direct laser removal of the film material. This recording technology is characterised by the absence of thermal and hydrodynamic image distortions, as in the case of laser ablation, and the number of necessary technological operations is much smaller as compared with the photomask preparation in classical photolithography. The main field of application of the thermochemical technology is the fabrication of diffraction optical elements (DOEs), which are widely used in printers, bar-code readers, CD and DVD laser players, etc. The purpose of this study is to increase the resolution of thermochemical data recording on thin chromium films. (interaction of laser radiation with matter)

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

GLOBAL OPTIMIZATION METHODS FOR GRAVITATIONAL LENS SYSTEMS WITH REGULARIZED SOURCES

Several approaches exist to model gravitational lens systems. In this study, we apply global optimization methods to find the optimal set of lens parameters using a genetic algorithm. We treat the full optimization procedure as a two-step process: an analytical description of the source plane intensity distribution is used to find an initial approximation to the optimal lens parameters; the second stage of the optimization uses a pixelated source plane with the semilinear method to determine an optimal source. Regularization is handled by means of an iterative method and the generalized cross validation (GCV) and unbiased predictive risk estimator (UPRE) functions that are commonly used in standard image deconvolution problems. This approach simultaneously estimates the optimal regularization parameter and the number of degrees of freedom in the source. Using the GCV and UPRE functions, we are able to justify an estimation of the number of source degrees of freedom found in previous work. We test our approach by applying our code to a subset of the lens systems included in the SLACS survey.

79 ASTRONOMY AND ASTROPHYSICS↗