Strange behavior with l_mperp=T although l_noco=F

Summary

Two calculations, which only differ by the switch l_mperp=F/T, yield completely different results:

l_mperp=F starts with a reasonable charge distance after the generation of a starting density and converges to the correct magnetic moment of 1.9 bohr
l_mperp=T starts at a super high distance and converges to a non-magnetic state.

Strangely, the flag l_mperp is only supposed to affect the calculation if l_noco=T. However, l_noco=F in both cases.

Input and output file:

inp.xml_bug out.xml_bug out_bug

inp.xml_correct out.xml_correct out_correct

Compute environment

************************Summary***************************
Compiler   :  /usr/local/intel/compilers_and_libraries_2019.3.199/linux/mpi/intel64/bin/mpiifort
Compiler ID:  Intel
Flags      :   -mkl  -mkl -mkl -qopenmp -assume byterecl -r8
Libraries  :  -lxml2;-lmkl_scalapack_lp64;-lmkl_blacs_intelmpi_lp64;-mt_mpi

These Libraries are required:
  XML Library found       :     TRUE
  LAPACK Library found    :     TRUE
These Libraries are optional:
  ELPA (one node) Library found :     FALSE
  FFT from MKL found            :     TRUE
  LibXC Library found           :     FALSE
  HDF5 Library found            :     TRUE
  Wannier90 1.2 Library found   :     FALSE
  Wannier90-4   Library found   :     
  Wannier90-5   Library found   :     
  MAGMA Library found           :     FALSE
  MPI Library found             :     TRUE
    SCALAPACK Library found     :     TRUE
    ELPA Library found          :     FALSE
    ChASE Library found         :     FALSE
  Compile GPU version           :     FALSE


Compile serial version        :     FALSE
Compile parallel version      :     TRUE


Git describe  : snapshot-02.19-612-gbb20fec
Git hash      : bb20fec99f000d0a6e0016724afbb69c5d838714
Doxygen found : YES
************************-------**************************

This is BUG because:

Cobalt is magnetic

# Summary

Two calculations, which only differ by the switch `l_mperp=F/T`, yield completely different results:
- `l_mperp=F` starts with a reasonable charge distance after the generation of a starting density and converges to the correct magnetic moment of 1.9 bohr
- `l_mperp=T` starts at a super high distance and converges to a non-magnetic state.

Strangely, the flag `l_mperp` is only supposed to affect the calculation if `l_noco=T`. However, `l_noco=F` in both cases.

## Input and output file:
[inp.xml_bug](/uploads/0ee6a13d733c84026da65890307db1a1/inp.xml_bug)
[out.xml_bug](/uploads/dd2814a7bc59c403fda280c834eab446/out.xml_bug)
[out_bug](/uploads/dc8b55c520684ed1137a0d950185b712/out_bug)

[inp.xml_correct](/uploads/e129b308d85c8da16f26490e9fe66f58/inp.xml_correct)
[out.xml_correct](/uploads/73893a69f88cf9c4ae8e267c79b31db9/out.xml_correct)
[out_correct](/uploads/8bdfec297730857900a4c3719666c2b3/out_correct)

## Compute environment
```
************************Summary***************************
Compiler   :  /usr/local/intel/compilers_and_libraries_2019.3.199/linux/mpi/intel64/bin/mpiifort
Compiler ID:  Intel
Flags      :   -mkl  -mkl -mkl -qopenmp -assume byterecl -r8
Libraries  :  -lxml2;-lmkl_scalapack_lp64;-lmkl_blacs_intelmpi_lp64;-mt_mpi

These Libraries are required:
  XML Library found       :     TRUE
  LAPACK Library found    :     TRUE
These Libraries are optional:
  ELPA (one node) Library found :     FALSE
  FFT from MKL found            :     TRUE
  LibXC Library found           :     FALSE
  HDF5 Library found            :     TRUE
  Wannier90 1.2 Library found   :     FALSE
  Wannier90-4   Library found   :     
  Wannier90-5   Library found   :     
  MAGMA Library found           :     FALSE
  MPI Library found             :     TRUE
    SCALAPACK Library found     :     TRUE
    ELPA Library found          :     FALSE
    ChASE Library found         :     FALSE
  Compile GPU version           :     FALSE

Compile serial version        :     FALSE
Compile parallel version      :     TRUE

Git describe  : snapshot-02.19-612-gbb20fec
Git hash      : bb20fec99f000d0a6e0016724afbb69c5d838714
Doxygen found : YES
************************-------**************************
```

## This is BUG because:

Cobalt is magnetic

Discussion
Designs

The one place for your designs

To enable design management, you'll need to meet the requirements. If you need help, reach out to our support team for assistance.