Dear Dr. Minh,
I am running ModelFinder on a nt dataset (with only 1st and 2nd codon positions, 3rd codons removed) using IQtree version 1.6.3. The command line I used is as following:
.......
WARNING: Numerical underflow for lh-branch
WARNING: Numerical underflow for lh-branch
WARNING: Numerical underflow for lh-branch
WARNING: Numerical underflow for lh-branch
WARNING: Numerical underflow for lh-branch
WARNING: Numerical underflow for lh-branch
ERROR: rategammainvar.cpp:244: double RateGammaInvar::optimizeWithEM(double): Assertion `newPInvar < 1.0' failed.
ERROR:
ERROR: *** IQ-TREE CRASHES WITH SIGNAL ABORTED
ERROR: *** For bug report please send to developers:
ERROR: *** Log file: 2.1_5561genes.nt.1and2codons.IQtree.MFM10.log
ERROR: *** Alignment files (if possible)
/var/spool/slurmd/job2523971/slurm_script: line 20: 4113 Aborted iqtree -nt 30 -st DNA -s ChalCOS.5561genes.nt.1and2codons.fas -spp ChalCOS.5561genes.nt.1and2codons.partition.nex -m MF+MERGE -rcluster 10 -rcluster-max 10000 -pre 2.1_5561genes.nt.1and2codons.IQtree.MFM10 -cptime 60 -safe
I put the original fasta alignment and partition files in the dropbox folder (following link) for the trouble-shooting purpose. And in this folder, you can also find the .sh file I used for submitting the job to the biocluster and the output log file.