Numerical model for macroscopic quantum superpositions based on phase-covariant quantum cloning

Macroscopically populated quantum superpositions pose a question to what extent the macroscopic world obeys quantum mechanical laws. Recently, such superpositions for light, generated by an optimal quantum cloner, have been demonstrated. They are of fundamental and technological interest. We present numerical methods useful for modeling of these states. Their properties are governed by a Gaussian hypergeometric function, which cannot be reduced to either elementary or easily tractable functions. We discuss the method of efficient computation of this function for half-integer parameters and a moderate value of its argument. We show how to dynamically estimate a cutoff for infinite sums involving this function performed over its parameters. Our algorithm exceeds double precision and is parallelizable. Depending on the experimental parameters it chooses one of the several ways of summation to achieve the best efficiency. The methods presented here can be adjusted for analysis of similar experimental schemes. Program summary Program title: MQSVIS Catalogue identifier: AEMR_ v1_ 0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEMR_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 1643 No. of bytes in distributed program, including test data, etc.: 13212 Distribution format: tar.gz Programming language: C with OpenMP extensions (main numerical program), Python (helper scripts). Computer: Modern PC (tested on AMD and Intel processors), HP BL2x220. Operating system: Unix/Linux. Has the code been vectorized or parallelized?: Yes (OpenMP). RAM: 200 MB for single run for 1000×1000 tile Classification: 4.15, 18. External routines: OpenMP Nature of problem: Recently, macroscopically populated quantum superpositions for light, generated by an optimal quantum cloner, have been demonstrated. They are of fundamental and technological interest. Their properties are governed by Gaussian hypergeometric function 2F₁ of half-integer parameters, which cannot be reduced to either elementary or easily tractable functions. Computation of the photon number distribution, visibility, and mean number of photons, necessary for characterization of these states, requires evaluation of infinite sums involving this function performed over its parameters. Solution method: The MQSVIS program suite computes various quantum indicators, such as photon number distribution, visibility, mean number of photons, and variance, for macroscopic quantum superpositions of light. It takes losses (modeled with a beam splitter) and imperfect photodetection (modeled with a Weierstrass transform applied to the photon number distribution) into account. Cutoffs of the infinite hypergeometric sums are estimated dynamically, and precision is enhanced with the computation of expressions in logarithmic form. Depending on the experimental parameters, the program chooses one of several ways of summation to achieve the best efficiency. The program is parallelized using OpenMP standard, which ensures the best utilization of multicore processors, and splits the work into tiles computed with different nodes of a computer cluster. This allows computation of the required indicators for realistic values of the parameters. Running time: 1-2 h for a 1000×1000 tile, depending on the values of the parameters.