A framework for implementation of organ effect models in TOPAS with benchmarks extended to proton therapy

The aim of this work was to develop a framework for modeling organ effects within TOPAS (TOol for PArticle Simulation), a wrapper of the Geant4 Monte Carlo toolkit that facilitates particle therapy simulation. The DICOM interface for TOPAS was extended to permit contour input, used to assign voxels to organs. The following dose response models were implemented: The Lyman-Kutcher-Burman model, the critical element model, the population based critical volume model, the parallel-serial model, a sigmoid-based model of Niemierko for normal tissue complication probability and tumor control probability (TCP), and a Poisson-based model for TCP. The framework allows easy manipulation of the parameters of these models and the implementation of other models. As part of the verification, results for the parallel-serial and Poisson model for x-ray irradiation of a water phantom were compared to data from the AAPM Task Group 166. When using the task group dose-volume histograms (DVHs), results were found to be sensitive to the number of points in the DVH, with differences up to 2.4%, some of which are attributable to differences between the implemented models. New results are given with the point spacing specified. When using Monte Carlo calculations with TOPAS, despite the relatively good match to the published DVH's, differences up to 9% were found for the parallel-serial model (for a maximum DVH difference of 2%) and up to 0.5% for the Poisson model (for a maximum DVH difference of 0.5%). However, differences of 74.5% (in Rectangle1), 34.8% (in PTV) and 52.1% (in Triangle) for the critical element, critical volume and the sigmoid-based models were found respectively. We propose a new benchmark for verification of organ effect models in proton therapy. The benchmark consists of customized structures in the spread out Bragg peak plateau, normal tissue, tumor, penumbra and in the distal region. The DVH's, DVH point spacing, and results of the organ effect models are provided. The models were used to calculate dose response for a Head and Neck patient to demonstrate functionality of the new framework and indicate the degree of variability between the models in proton therapy.

Improved efficiency in Monte Carlo simulation for passive-scattering proton therapy

The aim of this work was to improve the computational efficiency of Monte Carlo simulations when tracking protons through a proton therapy treatment head. Two proton therapy facilities were considered, the Francis H Burr Proton Therapy Center (FHBPTC) at the Massachusetts General Hospital and the Crocker Lab eye treatment facility used by University of California at San Francisco (UCSFETF). The computational efficiency was evaluated for phase space files scored at the exit of the treatment head to determine optimal parameters to improve efficiency while maintaining accuracy in the dose calculation. For FHBPTC, particles were split by a factor of 8 upstream of the second scatterer and upstream of the aperture. The radius of the region for Russian roulette was set to 2.5 or 1.5 times the radius of the aperture and a secondary particle production cut (PC) of 50 mm was applied. For UCSFETF, particles were split a factor of 16 upstream of a water absorber column and upstream of the aperture. Here, the radius of the region for Russian roulette was set to 4 times the radius of the aperture and a PC of 0.05 mm was applied. In both setups, the cylindrical symmetry of the proton beam was exploited to position the split particles randomly spaced around the beam axis. When simulating a phase space for subsequent water phantom simulations, efficiency gains between a factor of 19.9  ±  0.1 and 52.21  ±  0.04 for the FHTPC setups and 57.3  ±  0.5 for the UCSFETF setups were obtained. For a phase space used as input for simulations in a patient geometry, the gain was a factor of 78.6  ±  7.5. Lateral-dose curves in water were within the accepted clinical tolerance of 2%, with statistical uncertainties of 0.5% for the two facilities. For the patient geometry and by considering the 2% and 2mm criteria, 98.4% of the voxels showed a gamma index lower than unity. An analysis of the dose distribution resulted in systematic deviations below of 0.88% for 20% of the voxels with dose of 20% of the maximum or more.

Experimental validation of the TOPAS Monte Carlo system for passive scattering proton therapy

PURPOSE

TOPAS (TOol for PArticle Simulation) is a particle simulation code recently developed with the specific aim of making Monte Carlo simulations user-friendly for research and clinical physicists in the particle therapy community. The authors present a thorough and extensive experimental validation of Monte Carlo simulations performed with TOPAS in a variety of setups relevant for proton therapy applications. The set of validation measurements performed in this work represents an overall end-to-end testing strategy recommended for all clinical centers planning to rely on TOPAS for quality assurance or patient dose calculation and, more generally, for all the institutions using passive-scattering proton therapy systems.

METHODS

The authors systematically compared TOPAS simulations with measurements that are performed routinely within the quality assurance (QA) program in our institution as well as experiments specifically designed for this validation study. First, the authors compared TOPAS simulations with measurements of depth-dose curves for spread-out Bragg peak (SOBP) fields. Second, absolute dosimetry simulations were benchmarked against measured machine output factors (OFs). Third, the authors simulated and measured 2D dose profiles and analyzed the differences in terms of field flatness and symmetry and usable field size. Fourth, the authors designed a simple experiment using a half-beam shifter to assess the effects of multiple Coulomb scattering, beam divergence, and inverse square attenuation on lateral and longitudinal dose profiles measured and simulated in a water phantom. Fifth, TOPAS' capabilities to simulate time dependent beam delivery was benchmarked against dose rate functions (i.e., dose per unit time vs time) measured at different depths inside an SOBP field. Sixth, simulations of the charge deposited by protons fully stopping in two different types of multilayer Faraday cups (MLFCs) were compared with measurements to benchmark the nuclear interaction models used in the simulations.

RESULTS

SOBPs' range and modulation width were reproduced, on average, with an accuracy of +1, -2 and ±3 mm, respectively. OF simulations reproduced measured data within ±3%. Simulated 2D dose-profiles show field flatness and average field radius within ±3% of measured profiles. The field symmetry resulted, on average in ±3% agreement with commissioned profiles. TOPAS accuracy in reproducing measured dose profiles downstream the half beam shifter is better than 2%. Dose rate function simulation reproduced the measurements within ∼2% showing that the four-dimensional modeling of the passively modulation system was implement correctly and millimeter accuracy can be achieved in reproducing measured data. For MLFCs simulations, 2% agreement was found between TOPAS and both sets of experimental measurements. The overall results show that TOPAS simulations are within the clinical accepted tolerances for all QA measurements performed at our institution.

CONCLUSIONS

Our Monte Carlo simulations reproduced accurately the experimental data acquired through all the measurements performed in this study. Thus, TOPAS can reliably be applied to quality assurance for proton therapy and also as an input for commissioning of commercial treatment planning systems. This work also provides the basis for routine clinical dose calculations in patients for all passive scattering proton therapy centers using TOPAS.

Range verification of passively scattered proton beams based on prompt gamma time patterns

We propose a proton range verification technique for passive scattering proton therapy systems where spread out Bragg peak (SOBP) fields are produced with rotating range modulator wheels. The technique is based on the correlation of time patterns of the prompt gamma ray emission with the range of protons delivering the SOBP. The main feature of the technique is the ability to verify the proton range with a single point of measurement and a simple detector configuration. We performed four-dimensional (time-dependent) Monte Carlo simulations using TOPAS to show the validity and accuracy of the technique. First, we validated the hadronic models used in TOPAS by comparing simulations and prompt gamma spectrometry measurements published in the literature. Second, prompt gamma simulations for proton range verification were performed for the case of a water phantom and a prostate cancer patient. In the water phantom, the proton range was determined with 2 mm accuracy with a full ring detector configuration for a dose of ~2.5 cGy. For the prostate cancer patient, 4 mm accuracy on range determination was achieved for a dose of ~15 cGy. The results presented in this paper are encouraging in view of a potential clinical application of the technique.

Geometrical splitting technique to improve the computational efficiency in Monte Carlo calculations for proton therapy

PURPOSE

To present the implementation and validation of a geometrical based variance reduction technique for the calculation of phase space data for proton therapy dose calculation.

METHODS

The treatment heads at the Francis H Burr Proton Therapy Center were modeled with a new Monte Carlo tool (TOPAS based on Geant4). For variance reduction purposes, two particle-splitting planes were implemented. First, the particles were split upstream of the second scatterer or at the second ionization chamber. Then, particles reaching another plane immediately upstream of the field specific aperture were split again. In each case, particles were split by a factor of 8. At the second ionization chamber and at the latter plane, the cylindrical symmetry of the proton beam was exploited to position the split particles at randomly spaced locations rotated around the beam axis. Phase space data in IAEA format were recorded at the treatment head exit and the computational efficiency was calculated. Depth-dose curves and beam profiles were analyzed. Dose distributions were compared for a voxelized water phantom for different treatment fields for both the reference and optimized simulations. In addition, dose in two patients was simulated with and without particle splitting to compare the efficiency and accuracy of the technique.

RESULTS

A normalized computational efficiency gain of a factor of 10-20.3 was reached for phase space calculations for the different treatment head options simulated. Depth-dose curves and beam profiles were in reasonable agreement with the simulation done without splitting: within 1% for depth-dose with an average difference of (0.2 ± 0.4)%, 1 standard deviation, and a 0.3% statistical uncertainty of the simulations in the high dose region; 1.6% for planar fluence with an average difference of (0.4 ± 0.5)% and a statistical uncertainty of 0.3% in the high fluence region. The percentage differences between dose distributions in water for simulations done with and without particle splitting were within the accepted clinical tolerance of 2%, with a 0.4% statistical uncertainty. For the two patient geometries considered, head and prostate, the efficiency gain was 20.9 and 14.7, respectively, with the percentages of voxels with gamma indices lower than unity 98.9% and 99.7%, respectively, using 2% and 2 mm criteria.

CONCLUSIONS

The authors have implemented an efficient variance reduction technique with significant speed improvements for proton Monte Carlo simulations. The method can be transferred to other codes and other treatment heads.

GPU-based fast Monte Carlo dose calculation for proton therapy

Accurate radiation dose calculation is essential for successful proton radiotherapy. Monte Carlo (MC) simulation is considered to be the most accurate method. However, the long computation time limits it from routine clinical applications. Recently, graphics processing units (GPUs) have been widely used to accelerate computationally intensive tasks in radiotherapy. We have developed a fast MC dose calculation package, gPMC, for proton dose calculation on a GPU. In gPMC, proton transport is modeled by the class II condensed history simulation scheme with a continuous slowing down approximation. Ionization, elastic and inelastic proton nucleus interactions are considered. Energy straggling and multiple scattering are modeled. Secondary electrons are not transported and their energies are locally deposited. After an inelastic nuclear interaction event, a variety of products are generated using an empirical model. Among them, charged nuclear fragments are terminated with energy locally deposited. Secondary protons are stored in a stack and transported after finishing transport of the primary protons, while secondary neutral particles are neglected. gPMC is implemented on the GPU under the CUDA platform. We have validated gPMC using the TOPAS/Geant4 MC code as the gold standard. For various cases including homogeneous and inhomogeneous phantoms as well as a patient case, good agreements between gPMC and TOPAS/Geant4 are observed. The gamma passing rate for the 2%/2 mm criterion is over 98.7% in the region with dose greater than 10% maximum dose in all cases, excluding low-density air regions. With gPMC it takes only 6-22 s to simulate 10 million source protons to achieve ∼1% relative statistical uncertainty, depending on the phantoms and energy. This is an extremely high efficiency compared to the computational time of tens of CPU hours for TOPAS/Geant4. Our fast GPU-based code can thus facilitate the routine use of MC dose calculation in proton therapy.

A modular method to handle multiple time-dependent quantities in Monte Carlo simulations

A general method for handling time-dependent quantities in Monte Carlo simulations was developed to make such simulations more accessible to the medical community for a wide range of applications in radiotherapy, including fluence and dose calculation. To describe time-dependent changes in the most general way, we developed a grammar of functions that we call 'Time Features'. When a simulation quantity, such as the position of a geometrical object, an angle, a magnetic field, a current, etc, takes its value from a Time Feature, that quantity varies over time. The operation of time-dependent simulation was separated into distinct parts: the Sequence samples time values either sequentially at equal increments or randomly from a uniform distribution (allowing quantities to vary continuously in time), and then each time-dependent quantity is calculated according to its Time Feature. Due to this modular structure, time-dependent simulations, even in the presence of multiple time-dependent quantities, can be efficiently performed in a single simulation with any given time resolution. This approach has been implemented in TOPAS (TOol for PArticle Simulation), designed to make Monte Carlo simulations with Geant4 more accessible to both clinical and research physicists. To demonstrate the method, three clinical situations were simulated: a variable water column used to verify constancy of the Bragg peak of the Crocker Lab eye treatment facility of the University of California, the double-scattering treatment mode of the passive beam scattering system at Massachusetts General Hospital (MGH), where a spinning range modulator wheel accompanied by beam current modulation produces a spread-out Bragg peak, and the scanning mode at MGH, where time-dependent pulse shape, energy distribution and magnetic fields control Bragg peak positions. Results confirm the clinical applicability of the method.

Efficient voxel navigation for proton therapy dose calculation in TOPAS and Geant4

A key task within all Monte Carlo particle transport codes is 'navigation', the calculation to determine at each particle step what volume the particle may be leaving and what volume the particle may be entering. Navigation should be optimized to the specific geometry at hand. For patient dose calculation, this geometry generally involves voxelized computed tomography (CT) data. We investigated the efficiency of navigation algorithms on currently available voxel geometry parameterizations in the Monte Carlo simulation package Geant4: G4VPVParameterisation, G4VNestedParameterisation and G4PhantomParameterisation, the last with and without boundary skipping, a method where neighboring voxels with the same Hounsfield unit are combined into one larger voxel. A fourth parameterization approach (MGHParameterization), developed in-house before the latter two parameterizations became available in Geant4, was also included in this study. All simulations were performed using TOPAS, a tool for particle simulations layered on top of Geant4. Runtime comparisons were made on three distinct patient CT data sets: a head and neck, a liver and a prostate patient. We included an additional version of these three patients where all voxels, including the air voxels outside of the patient, were uniformly set to water in the runtime study. The G4VPVParameterisation offers two optimization options. One option has a 60-150 times slower simulation speed. The other is compatible in speed but requires 15-19 times more memory compared to the other parameterizations. We found the average CPU time used for the simulation relative to G4VNestedParameterisation to be 1.014 for G4PhantomParameterisation without boundary skipping and 1.015 for MGHParameterization. The average runtime ratio for G4PhantomParameterisation with and without boundary skipping for our heterogeneous data was equal to 0.97: 1. The calculated dose distributions agreed with the reference distribution for all but the G4PhantomParameterisation with boundary skipping for the head and neck patient. The maximum memory usage ranged from 0.8 to 1.8 GB depending on the CT volume independent of parameterizations, except for the 15-19 times greater memory usage with the G4VPVParameterisation when using the option with a higher simulation speed. The G4VNestedParameterisation was selected as the preferred choice for the patient geometries and treatment plans studied.

Measurement of inclusive radiative B-meson decays with a photon energy threshold of 1.7 GeV

Using 605 fb(-1) of data collected at the Upsilon(4S) resonance we present a measurement of the inclusive radiative B-meson decay channel, B-->X(s)gamma. For the lower photon energy thresholds of 1.7, 1.8, 1.9, and 2.0 GeV, as defined in the rest frame of the B meson, we measure the partial branching fraction and the mean and variance of the photon energy spectrum. At the 1.7 GeV threshold we obtain the partial branching fraction BF(B-->X(s)}gamma)=(3.45+/-0.15+/-0.40)x10(-4), where the errors are statistical and systematic.

Measurement of the Decay B_(s);(0)-->D_(s);(-)pi;(+) and Evidence for B_(s);(0)-->D_(s);(-/+)K;(+/-) in e;(+)e;(-) Annihilation at sqrt[s] approximately 10.87 GeV

We have studied B_(s);(0)-->D_(s);(-)pi;(+) and B_(s);(0)-->D_(s);(-/+)K;(+/-) decays using 23.6 fb;(-1) of data collected at the Upsilon(5S) resonance with the Belle detector at the KEKB e;(+)e;(-) collider. This highly pure B_(s);(0)-->D_(s);(-)pi;(+) sample is used to measure the branching fraction, B(B_(s);(0)-->D_(s);(-)pi;(+))=[3.67_(-0.33);(+0.35)(stat)-0.42+0.43(syst)+/-0.49(f_(s))]x10;(-3) (f_(s)=N_(B_(s);((*))B[over ]_(s);((*)))/N_(bb[over ])) and the fractions of B_(s);(0) event types at the Upsilon(5S) energy, in particular N_(B_(s);(*)B[over ]_(s);(*))/N_(B_(s);((*))B[over ]_(s);((*)))=(90.1_(-4.0);(+3.8)+/-0.2)%. We also determine the masses M(B_(s);(0))=(5364.4+/-1.3+/-0.7) MeV/c;(2) and M(B_(s);(*))=(5416.4+/-0.4+/-0.5) MeV/c;(2). In addition, we observe B_(s);(0)-->D_(s);(-/+)K;(+/-) decays with a significance of 3.5sigma and measure B(B_(s);(0)-->D_(s);(-/+)K;(+/-))=[2.4_(-1.0);(+1.2)(stat)+/-0.3(syst)+/-0.3(f_(s))]x10;(-4).

Evidence for neutral B meson decays to omegaK*0

We present the results of a study of the charmless vector-vector decay B0-->omegaK*0 with 657 x 10(6) BB pairs collected with the Belle detector at the KEKB e+e- collider. We measure the branching fraction to be B(B0-->omegaK*0) = [1.8+/-0.7(stat)+/-0.3(syst)]x 10(-6) with 3.0sigma significance. We also perform a helicity analysis of the omega and K*0 vector mesons, and obtain the longitudinal polarization fraction fL(B0-->omegaK*0) = 0.56+/-0.29(stat) -0.08(+0.18)(syst). Finally, we measure a large nonresonant branching fraction B[B0-->omegaK+pi(-);M Kpi in(0.755,1.250) GeV/c2] = [5.1+/-0.7(stat)+/-0.7(syst)]x 10(-6) with a significance of 9.5sigma.

Measurement of branching fractions, isospin, and CP-violating asymmetries for exclusive b-->dgamma modes

We report new measurements of the decays B+-->rho+gamma, B0-->rho0gamma, and B0-->omegagamma using a data sample of 657x10(6) B meson pairs accumulated with the Belle detector at the KEKB e+e- collider. We measure branching fractions B(B+-->rho+gamma)=(8.7_-2.7-1.1;+2.9+0.9)x10(-7), B(B0-->rho0gamma)=(7.8_-1.6-1.0;+1.7+0.9)x10(-7), and B(B0-->omegagamma)=(4.0_-1.7;+1.9+/-1.3)x10(-7). We also report the isospin asymmetry Delta(rhogamma)=-0.48_-0.19-0.09;+0.21+0.08 and the first measurement of the direct CP-violating asymmetry ACP(B+-->rho+gamma)=-0.11+/-0.32+/-0.09, where the first and second errors are statistical and systematic, respectively.

Search for B+-->D*+pi0 decay

We report on a search for the doubly Cabibbo suppressed decay B+-->D*+pi0, based on a data sample of 657x10(6) BB pairs collected at the Upsilon(4S) resonance with the Belle detector at the KEKB asymmetric-energy e+e- collider. We find no significant signal and set an upper limit of B(B+-->D*+pi0)<3.6x10(-6) at the 90% confidence level. This limit can be used to constrain the ratio between suppressed and favored B-->D*pi decay amplitudes, r<0.051, at the 90% confidence level.

Observation of B_(s)(0)-->phigamma and search for B_(s)(0)-->gammagamma decays at Belle

We search for the radiative penguin decays B_{s}{0}-->varphigamma and B_{s}{0}-->gammagamma in a 23.6 fb{-1} data sample collected at the Upsilon(5S) resonance with the Belle detector at the KEKB e{+}e{-} asymmetric-energy collider. We observe for the first time a radiative penguin decay of the B_{s}{0} meson in the B_{s}{0}-->varphigamma mode and we measure B(B_{s}{0}-->varphigamma)=(57_{-15}{+18}(stat)-11+12(syst))x10{-6}. No significant B_{s}{0}-->gammagamma signal is observed and we set a 90% confidence level upper limit of B(B_{s}{0}-->gammagamma)<8.7x10{-6}.

Observation of anomalous upsilon(1S)pi+pi(-) and upsilon(2S)pi+pi(-) production near the upsilon(5S) resonance

We report the first observation of e;{+}e;{-}-->Upsilon(1S)pi;{+}pi;{-}, Upsilon(2S)pi;{+}pi;{-}, and first evidence for e;{+}e;{-}-->Upsilon(3S)pi;{+}pi;{-}, Upsilon(1S)K+K-, near the peak of the Upsilon(5S) resonance at sqrt[s] approximately 10.87 GeV. The results are based on a data sample of 21.7 fb;{-1} collected with the Belle detector at the KEKB e;{+}e;{-} collider. Attributing the signals to the Upsilon(5S) resonance, the partial widths Gamma(Upsilon(5S)-->Upsilon(1S)pi;{+}pi;{-})=0.59+/-0.04(stat)+/-0.09(syst) MeV and Gamma(Upsilon(5S)-->Upsilon(2S)pi;{+}pi;{-})=0.85+/-0.07(stat)+/-0.16(syst) MeV are obtained from the observed cross sections. These values exceed by more than 2 orders of magnitude the previously measured partial widths for dipion transitions between lower Upsilon resonances.

Time-dependent CP-violating asymmetry in B0-->rho0gamma decays

We report the first measurement of CP-violation parameters in B0-->rho0gamma decays based on a data sample of 657x10(6)BB pairs collected with the Belle detector at the KEKB asymmetric-energy e+e- collider. We obtain the time-dependent and direct CP-violating parameters, Srho0gamma=-0.83+/-0.65(stat)+/-0.18(syst) and Arho0gamma=-0.44+/-0.49(stat)+/-0.14(syst), respectively.

Search for B --> h(*)nunu[over ] decays at Belle

We present a search for the rare decays B --> h(*)nunu[over ], where h(*) stands for a light meson. A data sample of 535 x 10{6} BB[over ] pairs collected with the Belle detector at the KEKB e{+}e{-} collider is used. Signal candidates are required to have an accompanying B meson fully reconstructed in a hadronic mode and signal side particles consistent with a single h(*) meson. No significant signal is observed and we set upper limits on the branching fractions at 90% confidence level. The limits on B{0} --> K{*0}nunu[over ] and B{+} --> K{+}nunu[over ] decays are more stringent than the previous constraints, while the first searches for B{0} --> K{0}nunu[over ], pi{0}nunu[over ], rho{0}nunu[over ], varphinunu[over ] and B{+} --> K{*+}nunu[over ], rho{+}nunu[over ] are reported.

Observation of B(0)-->D(*-)tau(+)nu(tau) decay at Belle

We report an observation of the decay B{0}-->D{*-}tau{+}nu{tau} in a data sample containing 535x10{6} BB pairs collected with the Belle detector at the KEKB asymmetric-energy e{+}e{-} collider. We find a signal with a significance of 5.2sigma and measure the branching fraction B(B{0}-->D{*-}tau{+}nu{tau})=(2.02{-0.37}{+0.40}(stat)+/-0.37(syst))%. This is the first observation of an exclusive B decay with a b-->ctaunu{tau} transition.

Measurement of the e+e- -->pi+pi- J/psi cross section via initial-state radiation at Belle

The cross section for e(+)e(-)-->pi(+)pi(-)J/psi between 3.8 and 5.5 GeV/c(2) is measured using a 548 fb(-1) data sample collected on or near the Upsilon(4S) resonance with the Belle detector at KEKB. A peak near 4.25 GeV/c(2), corresponding to the so called Y(4260), is observed. In addition, there is another cluster of events at around 4.05 GeV/c(2). A fit using two interfering Breit-Wigner shapes describes the data better than one that uses only the Y(4260), especially for the lower-mass side of the 4.25 GeV enhancement.

Measurement of Einstein-Podolsky-Rosen-type flavor entanglement in Upsilon(4S) --> B0 B0 decays

The neutral B meson pair produced at the Upsilon(4S) should exhibit a nonlocal correlation of the type discussed by Einstein, Podolsky, and Rosen. We measure this correlation using the time-dependent flavor asymmetry of semileptonic B(0) decays, which we compare with predictions from quantum mechanics and two local realistic models. The data are consistent with quantum mechanics, and inconsistent with the other models. Assuming that some B pairs disentangle to produce B(0) and B(0) with definite flavor, we find a decoherent fraction of 0.029 +/ -0.057, consistent with no decoherence.

Measurement of D0-D0 mixing parameters in D0 --> Ks pi+ pi- decays

We report a measurement of D0-D(0) mixing parameters in D(0) --> K(s)(0) pi(+) pi(-) decays using a time-dependent Dalitz-plot analysis. We first assume CP conservation and subsequently allow for CP violation. The results are based on 540 fb(-1) of data accumulated with the Belle detector at the KEKB e(+)e(-) collider. Assuming negligible CP violation, we measure the mixing parameters x = (0.80 +/- 0.29(-0.07-0.14)(+0.09+0.10))% and y = (0.33+/-0.24(-0.12-0.08)(+0.08+0.06))%, where the errors are statistical, experimental systematic, and systematic due to the Dalitz decay model, respectively. Allowing for CP violation, we obtain the CP-violating parameters |q / p| = 0.86(-0.29-0.03)(+0.30+0.06) +/- 0.08 and arg(q/p) = (-14(-18-3-4)(+16+5+2)) degrees .

Measurements of branching fractions for B --> K pi and B --> pi pi decays

We report measurements of branching fractions for B --> K pi and B --> pi pi decays based on a data sample of 449 x 10(6) BB[over] pairs collected at the Upsilon(4S) resonance with the Belle detector at the KEKB asymmetric-energy e(+)e(-) collider. We also measure the ratios of partial widths for B-->Kpi decays, namely R(c) identical with 2Gamma(B(+) --> K(+) pi(0))/Gamma(B(+) --> K(0) pi(+)) = 1.08+/-0.06+/-0.08 and R(n) identical with Gamma(B(0) --> K(+) pi(-))/2 Gamma(B(0) --> K(0) pi(0)) = 1.08+/-0.08+/-0.08, where the first and the second errors are statistical and systematic, respectively. These ratios are sensitive to enhanced electroweak penguin contributions from new physics; the new measurements are, however, consistent with standard model expectations.

Measurement of the tau lepton mass and an upper limit on the mass difference between tau+ and tau-

The mass of the tau lepton has been measured in the decay mode tau-->3pinutau using a pseudomass technique. The result obtained from 414 fb-1 of data collected with the Belle detector is Mtau=[1776.61+/-0.13(stat)+/-0.35(sys)] MeV/c2. The upper limit on the relative mass difference between positive and negative tau leptons is |Mtau+-Mtau-|/Mtau<2.8 x 10-4 at 90% confidence level.

Experimental constraints on the spin and parity of the Lambdac(2880)+

We report the results of several studies of the Lambda(c)(+)pi(+)pi(-)X final state in continuum e(+)e(-) annihilation data collected by the Belle detector. An analysis of angular distributions in Lambda(c)(2880)(+)-->Sigma(c)(2455)(0,++)pi(+,-) decays strongly favors a Lambda(c)(2880)(+) spin assignment of 5/2 over 3/2 or 1/2. We find evidence for Lambda(c)(2880)(+)-->Sigma(c)(2520)(0,++)pi(+,-) decay and measure the ratio of Lambda(c)(2880)(+) partial widths Gamma(Sigma(c)(2520)pi)/Gamma(Sigma(c)(2455)pi)=0.225+/-0.062+/-0.025. This value favors the Lambda(c)(2880)(+) spin-parity assignment of 5/2(+) over 5/2(-). We also report the first observation of Lambda(c)(2940)(+)-->Sigma(c)(2455)(0,++)pi(+,-) decay and measure Lambda(c)(2880)(+) and Lambda(c)(2940)(+) mass and width parameters. These studies are based on a 553 fb(-1) data sample collected at or near the Upsilon(4S) resonance at the KEKB collider.

Evidence for CP violation in B0-->D+D- decays

We report measurements of the branching fraction and CP violation parameters in B(0)-->D+ D- decays. The results are based on a data sample that contains 535 x 10(6) BB pairs collected at the Upsilon(4S) resonance, with the Belle detector at the KEKB asymmetric-energy e+e- collider. We obtain [1.97+/-0.20(stat) +/- 0.20(syst)] x 10(-4) for the branching fraction of B0-->D+D-. The measured values of the CP violation parameters are S=-1.13+/-0.37+/-0.09, A=0.91+/-0.23+/-0.06, where the first error is statistical and the second is systematic. We find evidence of CP violation in B0-->D+D- at the 4.1sigma confidence level. While the value of S is consistent with expectations from other measurements, the value of the parameter A favors large direct CP violation at the 3.2sigma confidence level, in contradiction to standard model expectations.

Measurement of CP asymmetry in a time-dependent Dalitz analysis of B0-->(rhopi)0 and a constraint on the quark mixing matrix angle phi2

We present a measurement of CP asymmetry using a time-dependent Dalitz plot analysis of B0-->pi+pi-pi0 decays based on a 414 fb(-1) data sample containing 449 x 10(6) BB pairs. The data was collected on the Upsilon(4S) resonance with the Belle detector at the KEKB asymmetric energy e+ e- collider. Combining our analysis with information on charged B decay modes, we perform a full Dalitz and isospin analysis and obtain a constraint on the CKM angle phi2, 68 degrees < phi2 < 95 degrees as the 68.3% confidence interval for the phi2 solution consistent with the standard model (SM). A large SM-disfavored region also remains.

Evidence for D0-D0 mixing

We observe evidence for D(0)-D(0) mixing by measuring the difference in the apparent lifetime when a D(0) meson decays to the CP eigenstates K(+)K(-) and pi(+)pi(-) and when it decays to the final state K(-)pi(+). We find the relative difference of the lifetimes y(CP) to be [1.31+/-0.32(stat)+/-0.25(syst)]%, 3.2 standard deviations from zero. We also search for a CP asymmetry between D(0) and D(0) decays; no evidence for CP violation is found. These results are based on 540 fb(-1) of data recorded by the Belle detector at the KEKB e(+)e(-) collider.

Observation of direct CP violation in B0 --> pi+pi- decays and model-independent constraints on the quark-mixing angle phi2

We report a new measurement of the time-dependent CP-violating parameters in B(0)-->pi(+)pi(-) decays with 535 x 10(6) BB pairs collected with the Belle detector at the KEKB asymmetric-energy e(+)e(-) collider operating at the Upsilon(4S) resonance. We find 1464+/-65 B(0)-->pi(+)pi(-) events and measure the CP-violating parameters S(pipi)=-0.61+/-0.10(stat)+/-0.04(syst) and A(pipi)=+0.55+/-0.08(stat)+/-0.05(syst). We observe large direct CP violation with a significance greater than 5 standard deviations for any S(pipi) value. Using isospin relations, we measure the Cabibbo-Kobayashi-Maskawa quark-mixing matrix angle phi(2)=(97+/-11) degrees for the solution consistent with the standard model and exclude the range 11 degrees <phi(2)<79 degrees at the 95% confidence level.

Observation of B decays to two kaons

Using 449x10(6) BB[over ] pairs collected with the Belle detector at the KEKB asymmetric-energy e(+)e(-) collider, we observe clear signals for B(+)-->K[over ](0)K(+) and B(0)-->K[over ](0)K(0) decays with 5.3sigma and 6.0sigma significance, respectively. We measure the branching fractions B(B(+)-->K[over ](0)K(+))=(1.22(-0.28-0.16)(+0.32+0.13))x10(-6) and B(B(0)-->K[over ](0)K(0))=(0.87(-0.20)(+0.25)+/-0.09)x10(-6), and partial-rate asymmetries A(CP)(B(+)-->K[over ](0)K(+))=0.13(-0.24)(+0.23)+/-0.02 and A(CP)(B(0)-->K[over ](0)K(0))=-0.58(-0.66)(+0.73)+/-0.04. From a simultaneous fit, we also obtain B(B(+)-->K(0)pi(+))=(22.8(-0.7)(+0.8)+/-1.3)x10(-6) and A(CP)(B(+)-->K(0)pi(+))=0.03+/-0.03+/-0.01. The first and second error in the branching fractions and the partial-rate asymmetries are statistical and systematic, respectively. No signal is observed for B(0)-->K(+)K(-) decays, and for this branching fraction, we set an upper limit of 4.1x10(-7) at the 90% confidence level.

Search for invisible decay of the Upsilon(1S)

We report results of a search for the invisible decay of the Upsilon(1S) via the Upsilon(3S)-->pi+ pi- Upsilon(1S) transition using a data sample of 2.9 fb-1 at the Upsilon(3S) resonance. The data were collected with the Belle detector at the KEKB asymmetric-energy e+ e- collider. No signal is found, and an upper limit for the branching fraction at the 90% confidence level is determined to be B(Upsilon(1S)-->invisible)<2.5 x 10(-3).

Observation of the decay B(0)-->J/psieta

We report the first observation of B(0)-->J/psieta decay. These results are obtained from a data sample that contains 449 x10(6) BB[over] pairs accumulated at the Upsilon(4S) resonance with the Belle detector at the KEKB asymmetric-energy e(+)e(-) collider. We observe a signal with a significance of 8.1 sigma and obtain a branching fraction of (9.5+/-1.7(stat)+/-0.8(syst)) x 10(-6).

Measurement of the near-threshold e+e- --> D(*)+/-D(*)-/+ cross section using initial-state radiation

We report a measurement of the exclusive e+e- -->D(*)+/-D*-/+ cross section as a function of center-of-mass energy near the D(*)+/-D*-/+ threshold with initial-state radiation. A partial reconstruction technique is used to increase the efficiency and to suppress background. The analysis is based on a data sample collected with the Belle detector with an integrated luminosity of 547.8 fb(-1).

Observation of time-dependent CP violation in B0 --> eta'K0 decays and improved measurements of CP asymmetries in B0 --> phiK0, KS0KS0KS0 and B0 --> J/psiK0 decays

We present improved measurements of CP-violation parameters in B(0) --> phiK(0), eta(')K(0), KS(0)KS(0)KS(0) decays based on a sample of 535 x 10(6) BB pairs collected at the Upsilon(4S) resonance with the Belle detector at the KEKB energy-asymmetric e(+)e(-) collider. We obtain sin2phi1(eff)=+0.64+/-0.10(stat)+/-0.04(syst) for B(0) --> eta(')K(0), +0.50+/-0.21(stat)+/-0.06(syst) for B(0) --> phiK(0), and +0.30+/-0.32(stat)+/-0.08(syst) for B(0) --> KS(0)KS(0)KS(0) decays. We have observed CP violation in the B(0) --> eta(')K(0) decay with a significance of 5.6 standard deviations. We also perform an improved measurement of CP asymmetries in B(0) --> J/psiK(0) decays and obtain sin2phi1=+0.642+/-0.031(stat)+/-0.017(syst).

Evidence for B-->eta'pi and improved measurements for B-->eta'K

We report evidence for the exclusive two-body charmless hadronic B meson decay B-->eta'pi, and improved measurements of B-->eta'K. The results are obtained from a data sample of 386x10(6) BB pairs collected at the Upsilon(4S) resonance, with the Belle detector at the KEKB asymmetric-energy e+e- collider. We measure B(B+-->eta'pi+)=[1.76(-0.62)(+0.67)(stat)(-0.14)(+0.15)(syst)]x10(-6) and B(B0-->eta'pi0)=[2.79(-0.96)(+1.02)(stat)(-0.34)(+0.25)(syst)]x10(-6). We also find the ratio of B(B+-->eta'K+)/B(B0-->eta'K0)=1.17+/-0.08(stat)+/-0.03(syst) and measure the direct CP asymmetries for the charged modes.

Evidence for large direct CP violation in B+/- --> rho(770)0K+/- from analysis of three-body charmless B+/- --> K+/-pi+/-pi+/- decays

We report results on a Dalitz analysis of three-body charmless B+/- --> K+/-pi+/-pi+/- decay including searches for direct CP violation. We report the first observation of the decay B+/- --> f2(1270)K+/- with a statistical significance above 6sigma. We also observe first evidence for large direct CP violation in the B+/- --> rho(770)0K+/- channel. The results are obtained with a data sample that contains 386 10(6) BB pairs collected at the Y(4s) resonance with the Belle detector at the KEKB asymmetric-energy e+e- collider.

Measurement of azimuthal asymmetries in inclusive production of hadron pairs in e(+)e(-) annihilation at Belle

The Collins effect connects transverse quark spin with a measurable azimuthal dependence in the yield of hadronic fragments around the quark's momentum vector. Using two different reconstruction methods, we find evidence of statistically significant azimuthal asymmetries for charged pion pairs in e(+)e(-) annihilation at a center-of-mass energy of 10.52 GeV, which can be attributed to a transverse polarization of the primordial quarks. The measurement was performed using a sample of 79 x 10(6) hadronic events collected with the Belle detector.

Observation of b-->dgamma and determination of /Vtd/Vts/

We report the observation of the flavor-changing neutral current process b-->dgamma using a sample of 386 x 106 B meson pairs accumulated by the Belle detector at the KEKB e+e- collider. We measure branching fractions for the exclusive modes B--->rho-gamma, B0rho0gamma, and B0omegagamma. Assuming that these three modes are related by isospin, we find B(B-->(rho,gamma)gamma)=[Formula: See Text] with a significance of 5.1sigma. This result is used to determine the ratio of Cabibbo-Kobayashi-Maskawa matrix elements /Vtd/Vts/ to be [Formula: See Text].

Improved evidence for direct CP violation in B0-->pi+pi- decays and model-independent constraints on phi2

We present a new measurement of the time-dependent CP-violating parameters in B(0)--> pi(+)pi(-) decays with 275 x 10(6) BB pairs collected with the Belle detector at the KEKB asymmetric-energy e(+)e(-) collider operating at the Gamma(4S) resonance. We find 666 +/- 43 B(0) --> pi(+)pi(-) events and measure the CP-violating parameters: S(pipi) = -0.67 +/- 0.16(stat) +/- 0.06(syst) and A(pipi) = +0.56 +/- 0.12(stat) +/- 0.06(syst). We find evidence for large direct CP violation with a significance greater than 4 standard deviations for any S(pipi) value. Using isospin relations, we obtain 95.4% confidence intervals for the Cabibbo-Kobayashi-Maskawa quark-mixing matrix angle phi(2) of 0 degree < phiv(2) < 19 degrees and 71 degrees < phi(2) < 180 degrees.

Studies of CP violation in B-->J/PsiK* decays

CP violation in B-->J/PsiK* decays is studied using an angular analysis in a data sample of 253 fb(-1) recorded with the Belle detector at the KEKB e(+)e(-) collider. The flavor separated measurements of the decay amplitudes indicate no evidence for direct CP violation. T-odd CP violation is studied using the asymmetries in triple product correlations, and the results are consistent with the standard model null predictions. The time-dependent angular analysis gives the following values of CP-violating parameters: sin(2phi(1) = 0.24 +/- 0.31 +/- 0.05 and cos(2phi(1)=0.56 +/- 0.79 +/- 0.11.

Observation of B+ --> plambdagamma

We report the first observation of the radiative hyperonic B decay B+ --> plambdagamma, using a 140 fb(-1) data sample recorded on upsilom(4S) resonance with the Belle detector at the KEKB asymmetric energy e+e- collider. The measured branching fraction is [symbol: see text](B+ --> plambdagamma) = (2.16(+0.58)(-0.53) +/- 0.20) x 10(-6). We examine its M(plambda) distribution and observe a peak near threshold. This feature is expected by the short-distance b --> sgamma transition. A search for B+ --> pepsilon0gamma yields no significant signal, and we set a 90% confidence-level upper limit on the branching fraction of [symbol: see textB+ --> pepsilon0gamma <4.6 x 10(-6).