Combining Legendre's Polynomials and Genetic Algorithm in the Solution of Nonlinear Initial-Value Problems

This work develops a method for solving ordinary differential equations, that is, initial-value problems, with solutions approximated by using Legendre's polynomials. An iterative procedure for the adjustment of the polynomial coefficients is developed, based on the genetic algorithm. This procedure is applied to several examples providing comparisons between its results and the best polynomial fitting when numerical solutions by the traditional Runge-Kutta or Adams methods are available. The resulting algorithm provides reliable solutions even if the numerical solutions are not available, that is, when the mass matrix is singular or the equation produces unstable running processes.

THE GEMINI NICI PLANET-FINDING CAMPAIGN: DISCOVERY OF A CLOSE SUBSTELLAR COMPANION TO THE YOUNG DEBRIS DISK STAR PZ Tel

We report the discovery of a tight substellar companion to the young solar analog PZ Tel, a member of the beta Pic moving group observed with high-contrast adaptive optics imaging as part of the Gemini Near-Infrared Coronagraphic Imager Planet-Finding Campaign. The companion was detected at a projected separation of 16.4 +/- 1.0 AU (0.'' 33 +/- 0.'' 01) in 2009 April. Second-epoch observations in 2010 May demonstrate that the companion is physically associated and shows significant orbital motion. Monte Carlo modeling constrains the orbit of PZ Tel B to eccentricities >0.6. The near-IR colors of PZ Tel B indicate a spectral type of M7 +/- 2 and thus this object will be a new benchmark companion for studies of ultracool, low-gravity photospheres. Adopting an age of 12(-4)(+8) Myr for the system, we estimate a mass of 36 +/- 6 M(Jup) based on the Lyon/DUSTY evolutionary models. PZ Tel B is one of the few young substellar companions directly imaged at orbital separations similar to those of giant planets in our own solar system. Additionally, the primary star PZ Tel A shows a 70 mu m emission excess, evidence for a significant quantity of circumstellar dust that has not been disrupted by the orbital motion of the companion.

Detection of high-velocity material from the wind-wind collision zone of Eta Carinae across the 2009.0 periastron passage

We report near-infrared spectroscopic observations of the Eta Carinae massive binary system during 2008-2009 using the CRIRES spectrograph mounted on the 8m UT 1 Very Large Telescope (VLT Antu). We detect a strong, broad absorption wing in He I lambda 10833 extending up to -1900 km s(-1) across the 2009.0 spectroscopic event. Analysis of archival Hubble Space Telescope/Space Telescope Imaging Spectrograph ultraviolet and optical data identifies a similar high-velocity absorption (up to -2100 km s(-1)) in the ultraviolet resonance lines of Si IV lambda lambda 1394, 1403 across the 2003.5 event. Ultraviolet resonance lines from low-ionization species, such as Si II lambda lambda 1527, 1533 and CII lambda lambda 1334, 1335, show absorption only up to -1200 km s(-1), indicating that the absorption with velocities -1200 to -2100 km s(-1) originates in a region markedly more rapidly moving and more ionized than the nominal wind of the primary star. Seeing-limited observations obtained at the 1.6m OPD/LNA telescope during the last four spectroscopic cycles of Eta Carinae (1989-2009) also show high-velocity absorption in He I lambda 10833 during periastron. Based on the large OPD/LNA dataset, we determine that material with velocities more negative than -900 km s(-1) is present in the phase range 0.976 <= phi <= 1.023 of the spectroscopic cycle, but absent in spectra taken at phi <= 0.947 and phi >= 1.049. Therefore, we constrain the duration of the high-velocity absorption to be 95 to 206 days (or 0.047 to 0.102 in phase). We propose that the high-velocity absorption component originates in shocked gas in the wind-wind collision zone, at distances of 15 to 45 AU in the line-of-sight to the primary star. With the aid of three-dimensional hydrodynamical simulations of the wind-wind collision zone, we find that the dense high-velocity gas is along the line-of-sight to the primary star only if the binary system is oriented in the sky such that the companion is behind the primary star during periastron, corresponding to a longitude of periastron of omega similar to 240 degrees-270 degrees. We study a possible tilt of the orbital plane relative to the Homunculus equatorial plane and conclude that our data are broadly consistent with orbital inclinations in the range i = 40 degrees-60 degrees.

Nature of the X(3872) from QCD

We have studied some possible four-quark and molecule configurations of the X(3872) using double ratios of sum rules, which are more accurate than the usual simple ratios often used in the literature to obtain hadron masses. We found that the different structures ((3) over bar - (3) over bar and (6) over bar - 6 tetraquarks and D - D(*) molecule) lead to the same prediction for the mass (within the accuracy of the method), indicating that the alone prediction of the X mass may not be sufficient to reveal its nature. In doing these analyses, we also find that (within our approximation) the use of the (MS) over bar running (m) over bar (c)(m(c)(2)), rather than the on-shell mass, is more appropriate to obtain the J/psi and X meson masses. Using vertex sum rules to roughly estimate the X(3872) hadronic and radiative widths, we found that the available experimental data does not exclude a lambda - J/psi-like molecule current.

Radiative decay of the X(3872) as a mixed molecule-charmonium state in QCD sum rules

We use QCD sum rules (QCDSR) to calculate the width of the radiative decay of the meson X(3872), assumed to be a mixture between charmonium and exotic molecular [c (q) over bar][q (c) over bar] states with J(PC) = 1(++). We find that in a small range for the values of the mixing angle, 5 degrees <= theta <= 13 degrees, we get the branching ratio Gamma(X -> J/psi gamma)/Gamma(X -> J/psi pi(+)pi(-)) = 0.19 +/- 0.13, which is in agreement, with the experimental value. This result is compatible with the analysis of the mass and decay width of the mode J/psi(n pi) performed in the same approach.

Bs0 meson and the Bs0BK coupling from QCD sum rules

We evaluate the mass of the B(s0) scalar meson and the coupling constant in the B(s0)BK vertex in the framework of QCD sum rules. We consider the B(s0) as a tetraquark state to evaluate its mass. We get m(Bs0) = (5.85 +/- 0.13) GeV, which is in agreement, considering the uncertainties, with predictions supposing it as a b (s) over bar state or a B (K) over bar bound state with J(P) = 0(+). To evaluate the g(Bs0BK) coupling, we use the three-point correlation functions of the vertex, considering B(s0) as a normal b (s) over bar state. The obtained coupling constant is: g(Bs0BK) = (16.3 +/- 3.2) GeV. This number is in agreement with light-cone QCD sum rules calculation. We have also compared the decay width of the B(s0) -> BK process considering the B(s0) to be a b (s) over bar state and a BK molecular state. The width obtained for the BK molecular state is twice as big as the width obtained for the b (s) over bar state. Therefore, we conclude that with the knowledge of the mass and the decay width of the B(s0) meson, one can discriminate between the different theoretical proposals for its structure.

QCD sum rules for the X(3872) as a mixed molecule-charmonium state

We use QCD sum rules to test the nature of the meson X(3872), assumed to be a mixture between charmonium and exotic molecular [c (q) over bar][q (c) over bar] states with J(PC) = 1(++). We find that there is only a small range for the values of the mixing angle theta that can provide simultaneously good agreement with the experimental value of the mass and the decay width, and this range is 5(0) <= theta <= 3(0). In this range we get m(X) = (3.77 +/- 0.18) GeV and Gamma(X -> J/psi pi(+)pi(-)) = (9.3 +/- 6.9) MeV, which are compatible, within the errors, with the experimental values. We, therefore, conclude that the X(3872) is approximately 97% a charmonium state with 3% admixture of similar to 88% D(0)D*(0) molecule and similar to 12% D(+)D*(-) molecule.

Model Independent Tests of Skyrmions and Their Holographic Cousins

We describe a new exact relation for large N(c) QCD for the long-distance behavior of baryon form factors in the chiral limit. This model-independent relation is used to test the consistency of the structure of several baryon models. All 4D semiclassical chiral soliton models satisfy the relation, as does the Pomarol-Wulzer holographic model of baryons as 5D Skyrmions. However, remarkably, we find that the holographic model treating baryons as instantons in the Sakai-Sugimoto model does not satisfy the relation.

Width of exotics from QCD sum rules : Tetraquarks or molecules?

We investigate the widths of the recently observed charmonium like resonances X(3872), Z(4430), and Z(2)(4250) using QCD sum rules. Extending previous analyses regarding these states as diquark-antiquark states or molecules of D mesons, we introduce the Breit-Wigner function in the pole term. We find that introducing the width increases the mass at the small Borel window region. Using the operator-product expansion up to dimension 8, we find that the sum rules based on interpolating current with molecular components give a stable Borel curve from which both the masses and widths of these resonances can be well obtained. Thus the QCD sum rule approach strongly favors the molecular description of these states.

Identifying Multiquark Hadrons from Heavy Ion Collisions

Identifying hadronic molecular states and/or hadrons with multiquark components either with or without exotic quantum numbers is a long-standing challenge in hadronic physics. We suggest that studying the production of these hadrons in relativistic heavy ion collisions offers a promising resolution to this problem as yields of exotic hadrons are expected to be strongly affected by their structures. Using the coalescence model for hadron production, we find that, compared to the case of a nonexotic hadron with normal quark numbers, the yield of an exotic hadron is typically an order of magnitude smaller when it is a compact multiquark state and a factor of 2 or more larger when it is a loosely bound hadronic molecule. We further find that some of the newly proposed heavy exotic states could be produced and realistically measured in these experiments.

A Geometric Approach for Turbo Decoding of Reed-Solomon Codes in QAM Modulation Schemes

This work studies the turbo decoding of Reed-Solomon codes in QAM modulation schemes for additive white Gaussian noise channels (AWGN) by using a geometric approach. Considering the relations between the Galois field elements of the Reed-Solomon code and the symbols combined with their geometric dispositions in the QAM constellation, a turbo decoding algorithm, based on the work of Chase and Pyndiah, is developed. Simulation results show that the performance achieved is similar to the one obtained with the pragmatic approach with binary decomposition and analysis.

Extension of a Recent Method for Obtaining Exact Solutions of the Bruce and Klute Equation

A method based on a specific power-law relationship between the hydraulic head and the Boltzmann variable was recently presented. We generalized this relationship to a range of powers and extended the solution to include the saturated zone. As a result, the new solution satisfies the Bruce and Klute equation exactly.

New Analytic Solution Related to the Richards, Philip, and Green-Ampt Equations for Infiltration

Based on physical laws of similarity, an analytic solution of the soil water potential form of the Richards equation was derived for water infiltration into a homogeneous sand. The derivation assumes a similarity between the soil water retention function and that of the soil water content profiles taken at fixed times. The new solution successfully described soil water content profiles experimentally measured for water infiltrating downward, upward, and horizontally into a homogeneous sand and agrees with that presented by Philip in 1957. The utility of this analysis is still to be verified, but it is expected to hold for soils that have a narrow pore-size distribution before wetting and that manifest a sharp increase of water content at the wetting front during infiltration. The effect of van Genuchten`s parameters alpha and n on the application of the solution to other porous media was investigated. The solution also improves and provides a more realistic description of the infiltration process than that pioneered by Green and Ampt in 1911.

New Analytic Solution of Boltzmann Transform for Horizontal Water Infiltration into Sand

The use of the Boltzmann transform function, lambda(theta), to solve the Richards equation when the diffusivity, D, is a function of only soil water content,., is now commonplace in the literature. Nevertheless, a new analytic solution of the Boltzmann transform lambda(h) as a function of matric potential for horizontal water infiltration into a sand was derived without invoking the concept or use of D(theta). The derivation assumes that a similarity exists between the soil water retention function and the Boltzmann transform lambda(theta). The solution successfully described soil water content profiles experimentally measured for different infiltration times into a homogeneous sand and agrees with those presented by Philip in 1955 and 1957. The applicability of this solution for all soils remains open, but it is anticipated to hold for soils whose air-filled pore-size distribution before wetting is sufficiently narrow to yield a sharp increase of water content at the wetting front during infiltration. It also improves and provides a versatile alternative to the well-known analysis pioneered by Green and Ampt in 1911.

Alternative Analytical Expressions for the General van Genuchten-Mualem and van Genuchten-Burdine Hydraulic Conductivity Models

The van Genuchten expressions for the unsaturated soil hydraulic properties, first published in 1980, are used frequently in various vadose zone flow and transport applications assuming a specific relationship between the m and n soil hydraulic parameters. By comparison, probably because of the complexity of the hydraulic conductivity equations, the more general solutions with independent m and n values are rarely used. We expressed the general van Genuchten-Mualem and van Genuchten-Burdine hydraulic conductivity equations in terms of hypergeometric functions, which can be approximated by infinite series that converge rapidly for relatively large values of the van Genuchten-Mualem parameter n but only very slowly when n is close to one. Alternative equations were derived that provide very close approximations of the analytical results. The newly proposed equations allow the use of independent values of the parameters m and n in the soil water retention model of van Genuchten for subsequent prediction of the van Genuchten-Mualem and van Genuchten-Burdine hydraulic conductivity models, thus providing more flexibility in fitting experimental pressure-head-dependent water content, theta(h), and hydraulic conductivity, K(h), or K(theta) data.