Comments (3)
Hi Mike,
That error comes about in the removeRedundantPhases
function, which is needed when you trace the high-temperature evolution of two phases individually, and the two end up splitting from a common higher temperature phase via a phase transition. The higher temperature phase can get included twice, and this gets rid of one of the copies. Without getting rid of one copy you'd end up with micro transitions between nearly identical states. The problem occurs when you have two phases that join together in something that doesn't really look like a phase transition — dphi/dT remains finite — so the algorithm makes no distinction between the low-temperature phase and the high-temperature phase. So you end up with the temperature evolution of two different phases, both of which are the same at high temperature but different at low temperature.
The ideal solution would be to implement the algorithm in the not implemented error. I don't think this would be too hard to do, you just need to find the distance between the two phases as a function of temperature, mark where two diverge, and split them into low temperature and high temperature parts. Then you can get rid of one of the duplicate high-T pieces.
A much easier work around would just be to tighten your numerical tolerances. I don't think you can actually have two phases join each other without causing dphi/dT to diverge in any physically reasonable model (i.e., the second derivative of the field will go to zero in a second order transition), so if you tighten your tolerances you might be able to resolve the divergence. Try decreasing x_eps
and/or T_eps
in the generic_potential
class or your subclass thereof.
Let me know if you have any luck.
Cheers,
Carroll
from cosmotransitions.
Hi Carroll,
Thanks a lot for replying so quickly! Yes, ok, that's pretty much what I thought the error meant, from your exception description.
Thanks for sketching what the function should do. For now, as you suggested, I've found that I can improve things by varying x_eps (although in my case T_eps doesn't seem to have any effect).
I'll stick to this for now, but will let you know if I try to implement this function in the future.
Thanks again,
Mike
from cosmotransitions.
I am running cosmoTransitions for vev values 10^(3 to 7) GeV and kept the x_eps and T_eps as
self.x_eps = .00001
self.T_eps = .001
Still I am getting
NotImplementedError: Two phases have been found to coincide at one end but not the other. Ideally, this function would find where the two diverge, make a cut, and join them such there are no more phase redundancies.
Instead, just raise an exception.
Any suggestions
from cosmotransitions.
Related Issues (20)
- Eigenvalues of field-dependent mass automatically HOT 3
- getPhases() -- TraceMultiMin Error
- Issues with critical temperatures and minimum at T=0 HOT 5
- Queries on form of Boson_mass_sq
- Pass numerical potential to class SingleFieldInstanton to tunneling1D.py HOT 4
- Import issue in finiteT HOT 2
- RuntimeWarning: invalid value encountered in add
- Daisy improved corrections to the effective potential? HOT 1
- INPUT REHEATING TEMPERATURE
- raise TypeError( TypeError: object of type <class 'numpy.float64'> cannot be safely interpreted as an integer. HOT 1
- Compilation error HOT 1
- How to demtermine a thick-walled bubble or thin-walled bubble HOT 1
- A problem on diagonalization of mass matrix automatically HOT 6
- Cannot resolve bubble profile, keep getting "Integration Error: r > rmax"
- Euclidean action as a function of temperature (and derivative) HOT 1
- Error while calculating d(S3/T)/dT at T_nucleation HOT 3
- Low vev issue for 2nd-order phase transition
- Quick Question - Unit scales
- Issue with extend_to_minima in findAllTransitions()
- How to define path_pts (array_like) in full.tunneling() HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from cosmotransitions.